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Turbulence in Supersonic Flow 


LESLIE S. G. KOVASZNAY# 
The Johns Hopkins University 


SUMMARY 
First-order perturbation theory indicates that the Navier 
Stokes equations for a compressible, viscous, and heat-conductive 
gas can have three distinctly different types of disturbance fields, 
obeying three independent differential equations. These three 
“modes’’ vorticity mode, entropy 
mode, and sound-wave mode. 
the intensities of the fluctuations are small, but they interact at 


are: 
The modes are independent when 


of disturbance fields 


larger intensities when linearization is not permissible 

The principal tool for measuring the fluctuating fields is the 
hot-wire anemometer adapted for supersonic flow. In a super 
sonic flow the hot-wire anemometer responds to two variables: 
mass flow fluctuations and stagnation temperature fluctuations. 
Each of the above three modes contributes to both of these vari 
ables. By taking fluctuation measurements at different wire 
temperatures, the relative sensitivity to the flow parameters 
varies and a ‘‘fluctuation diagram’’ can be obtained. The char 
acter of the fluctuation diagram differs for each mode, and even 
coexisting modes can be separately determined by analyzing the 
fluctuation diagram under certain restrictive conditions 

Hot-wire measurements were obtained in various supersonic 
Pure temperature spottiness (entropy mode) and pure 
sound-wave fields were explored. Detailed measurements were 
taken in a turbulent supersonic boundary layer. With the aid of 
the fluctuation diagram, the intensity of temperature spottiness 
(entropy mode) and of turbulent velocity fluctuations (vorticity 
Energy spectrum and 


flows. 


mode) was measured across the layer. 
probability density measurements gave clues concerning the 
character of the fluctuations. 


INTRODUCTION 


p I ‘He PROBLEM OF TURBULENCE in a compressible 
medium brings up a number of specific questions. 
Will the turbulence be essentially similar to that in the 


Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953. 

* The author is much indebted to Drs. Francis H. Clauser and 
Mark V. Morkovin for their stimulation, vigorous interest, 
and helpful criticism; to Dr. F. R. Hama for the construction and 
operation of the high-frequency spectral analyzer; and to Miss 
Patricia C. Clarken for her active and efficient participation in 
all phases of the work. 
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incompressible case with only slight compressibility 
correction or will there be some entirely new phe- 
nomena? What kind of measurements of the fluctua- 
tions can be performed? What experimental methods 
can be used? How will the measurements be inter- 
preted? What quantity will take the role of the ‘‘tur- 
bulence level’ in the high-speed case? 

One can venture an intuitive guess that random fluc- 
tuations of more than one kind will be present in the 
compressible case. Sound waves are likely to accom- 
pany the vorticity fluctuations, and we may expect 
that the fluid particles, having passed through different 
shear regions, may have suffered different amounts of 
entropy changes. When the pressure becomes uniform, 
they may exhibit entropy spottiness appearing as tem- 
perature spottiness to the fixed observer. In the first 
part of the paper it will be proved that the three 
‘“‘modes” are the solutions for the perturbation fields 
superimposed on a uniform steady flow of a compres- 
sible, viscous, and heat-conductive gas. Since all 
three modes represent fluctuations both in momentum 
and energy when superimposed on a uniform super- 
sonic flow, the dynamic effects (e.g., the transition of a 
laminar boundary layer) are nonnegligible for the 
three modes. 

The extension of the hot-wire technique to super- 
sonic velocities has been reported earlier. The princi- 
pal objective of the present investigation is to find the 
proper method of disentangling the contributions to 
the hot-wire signal from the three different modes. If 
this effort is successful, the determination of the in- 
tensity and correlation of coexisting modes becomes 
possible. 

One of the most intriguing practical questions of 
supersonic turbulence measurements is the determina- 
tion of the fluctuation level in the free stream of a good 
supersonic wind tunnel and the exploration of the rela- 
tion of this level to the flow conditions upstream of the 
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nozzle. However, we have not chosen to investigate 
this case first because extremely low-level signals are 
involved which bring us dangerously close to the noise 
level of the electronic equipment. Intense fluctuations 
can be measured much more conveniently, and the 
character of ‘compressible turbulence’ can be more 
readily observed and understood. 

Following this reasoning, the first detailed measure- 
ments were obtained in a reasonably thick supersonic 
turbulent boundary layer. The scale of turbulence 
in this case was not prohibitively small and the inten- 
sity of fluctuations was sufficiently large to permit 
analysis of the signals and separation of contributions 
from the different modes. 

The present work is an integral part of our long-range 
program of turbulence investigation at supersonic 
velocities pursued for a number of years at Hopkins. 
During these years the financial support of two spon- 
sors made the investigations possible: U.S. Navy 
under contract NOrd-8036 and Squid Project under 
contract N6ori-105-Task Order III. Some of the ex- 
periments were conducted outside of our laboratory. 
The facilities of the Ballistic Research Laboratories, 
Aberdeen Proving Ground, were made available. The 
author also conducted some experiments at Ames 
Aeronautical Laboratory, NACA, and some of the 
results are reported here. The hot-wire equipment used 
in most of the tests was developed as a cooperative 
project of the National Bureau of Standards and The 
Johns Hopkins University. 


SYMBOLS 


position vector in coordinate system £4, £2, £3 


Fs - 


1» &2, &3 = Cartesian coordinates moving with the 
average velocity of the fluid (components 
of ) 

aS = Cartesian coordinates fixed with respect to 


the wind tunnel 
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T = relaxation time of medium 
¢ = velocity potential 
ly = characteristic length of the viscous com- 
pressible medium 
to = characteristic time of the medium 
r = wave length 
L = characteristic length of domain G 
k = wave number 
Wy = vorticity amplitude 
n = normal vector 
ty = decay time of vorticity mode 
d, = decay distance of vorticity mode 
P, = pressure disturbance amplitude 
Cc = complex propagation velocity 
ts = decay time of sound waves 
dy = decay distance of sound waves 
So = entropy disturbance amplitude 
ts = decay time for entropy waves 
d; = decay distance for entropy waves 
Ae(t) = voltage fluctuation across hot wire (signal 
compensated for thermal lag) 
m(t) mass flow fluctuation 
d(t) = stagnation temperature fluctuation 
To = stagnation temperature 
peta ( = velocity fluctuation observed by fixed probs 
- in x, y, z coordinate system 
AW(t) } . 
M = Mach Number 
a(M), B(M) = function of the Mach Number 
Aer = hot-wire sensitivity to stagnation tempera 
ture fluctuation 
Aem = hot-wire sensitivity to mass flow fluctuation 
yr = Aew/ Aer = sensitivity ratio 
F = (prime), root mean square (r.m.s.) value of a 
fluctuating quantity 
dx = Ae’/Aey = virtual r.m.s. stagnation temperature fluc 
tuation | 
x; = hot-wire resistance at equilibrium tempera 
ture 7, 
Rw = hot-wire resistance when operated at tem 
perature 7 jw 
, Rw —R. ? : 
aw’ = . = overheating ratio 
R. 
med ; 
Rur = — mass flow stagnation temperature correla 


tion coefficient 


= solenoidal velocity fluctuation (per cent ) 


temperature spottiness (entropy mode) in 


per cent of total temperature 








t = time 

p = density 

v = velocity vector in coordinate system &, , & 

UV, Ve, Vs = velocity components of v 

p = pressure 

B = dynamic viscosity 

r = absolute temperature 

G = specific heat at constant pressure 

K = heat conductivity 

oy = rate of dissipation 

¥ = ratio of specific heats 

v3 = Laplace operator 

G = space-time domain 

V = volume 

ta, te = beginning and end of time interval 
= bar indicating average value 

A = fluctuation 

r = nondimensional pressure fluctuation 

5 = nondimensional entropy fluctuation 

AE = entropy change 

vy = p/p = kinematic viscosity (average value) 

a= V u(p p) = speed of sound (average value) 

q = divu = source strength of the velocity field 

® = curlv = vorticity 


AU, = solenoidal velocity fluctuation 
Ty-1 
ATs =zAT—- —- ~— Ap = temperature spot 
P Y 
tiness 
Pn) = plane wave pressure distribution 
f = frequency 
o = mass flow fluctuation of sound waves that 
are almost Mach waves (per cent) 
Ls = integral scale of turbulence (longitudinal ) 
Uy, U2 } 


O71, 92, O3\ 


components of fluctuations 


(1) THREE MopeEs oF FLUCTUATIONS 


Governing Equations 


A simple analysis demonstrates that a viscous com- 


pressible heat-conducting medium can have three dif- 
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TURBULENCE 


ferent kinds of small perturbation fields obeying three 
independent differential equations. 

The following assumptions will (A) The 
medium is a perfect gas with constant specific heats, 
viscosity, and heat conduction; and (B) the Prandtl 
The governing 


be used: 


Number of the gas is equal to */,. 
equations of such a fluid contain four independent 
variables £), &, £3, ¢ and six dependent variables p, p, T, 
There are six equations: three equations 
(or a single vector equation) for the conservation of 
momentum, one equation for the conservation of 
energy, one for the conservation of mass, and one is 
the equation of state for the medium. 


Momentum: 
Ov p n 
 —e grad p — © grad v? + p (v X curl v) + 


u(V°v + '/3 grad div v) (1) 


Energy: 


i = «V°T + oP + v-grad p — pC,v-grad T + ® 
Of Of 
(2) 
Mass: 
(Op/Ot) + div (pv) = 0 (3) 
State: 
b/p = Uy — L/yIGT (4) 


9 


where V’ is the Laplace operator, 
0? 0” O° 
Of:> Of? ~—- OE; 


and ® is the viscous dissipation, 


® = »[(curl v)* + V2(v?) — */3(div v)? — 
2 div (v X curl v) — 2 v grad div v] 


The above set of partial differential equations are non- 
linear, and a simple inspection does not reveal any clues 
of the fluctuation fields. 


’ 


about different ‘‘modes’ 


Small Perturbation 


In the present analysis interest is focused on the 
small fluctuations in a flow field. Linearization of the 
above set of equations may yield an interesting limit- 
ing case. However, a similar analysis was carried out 
in reference 1 without heat-conduction effects. 

The following new assumptions will be introduced, 
in addition to (A) and (B): (C) The flow is described 


within a space time domain G consisting of a volume V 


in space and a time interval t, — tg. (D) The coordi- 
nate system is chosen in such a way that the average 
velocity in G is zero, (E) The fluctuations in pressure, 
density, and absolute temperature are small compared 
with their respective average values in G. (F) The 
velocities are small compared with the average speed of 
sound in G. (G) There are no solid boundaries within 


G or at its boundaries. 


IN SUPERSONIC 


FLOW 659 

The physical picture corresponding to the above as- 
sumptions is a compressible flow field in a finite domain 
of spave observed for a finite time interval viewed by 
an observer moving with the average velocity. For 
the sake of small perturbation theory, it is assumed 
that the departures from free-stream conditions are 
small. A typical case to which such a treatment may 
apply is the fluctuation field in a supersonic wind tun- 
nel. The flow is observed for a length of time that is 
less than the time of passage of a fluid particle through 
the test section. 

A trivial solution of the set of equations [Eqs. (1), 
(2), (3), (4)] is 

y == 0, p = const., T = const. 
p = yb/(y — 1)C,T = const. 

This ‘“‘flow field’ is a gas at rest in the above-defined 
coordinate system. From the conventional coordinate 
system x, y, 2, /, it also can represent a uniform parallel 
flow in a supersonic wind tunnel. 

We now assume small perturbations. Let us intro- 


duce mean values and fluctuations, 


T=T+ AT, p=pt+Ap, p=p+ Ap 


where 


oe 
T= + / / / / T (kx, >, fs, t) dk, do dé dt (5) 
mei ia a 


By definition 


Ap = Ap = 


The velocity fluctuation is the velocity itself by virtue 
of the special choice of coordinate system. 

Now it is convenient to define nondimensional tem- 
perature, pressure, and density fluctuations 


P = Ap/yp = Ap/pa’? (6) 
s = (AT/T) — (y — 1)P (7) 


where s represents the ‘“‘nonisentropic temperature 
fluctuation’’ and is really the nondimensional entropy 
fluctuation. 

Since the equation of state expressed with the en 
tropy change AF is 


1)/y,4E/Cp 


T2/T; = (pe pi)? ~ a 


linearization gives 
AT/T = (y — 1)P + (AE/C,) 
or 
s = AE/C, (8) 
The nondimensional density fluctuations become 


A A AT 
—— =P-s (9) 


p 6 Oop T 


Using the new variables the equations of motion can be 
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rewritten. Dropping the higher order terms we obtain (0qg/0t) — (4v/3)V7q = 0 (19 
a set of linear equations. (Os/Ot) — (4v/3)V*s = 0 (20 


Momentum: 
Ov . : 
> = —a’ grad P + vo(V*’u + 3/3 grad div v) (10) 
Energy: 
Os 4 4 — 1) 
_ sy a ty ” v2P (11) 
ot 3 3 
Mass: 
div v = —(OP/Ot) + (Os/Ot) (12) 


The linearized equation of state has been used to elimi- 
nate the density fluctuations. 

Having a set of linear equations, superposition holds; 
therefore, we are seeking different types of solutions 
We can obtain 
two different forms of the equation for the conservation 


to satisfy the above set of equations. 
of momentum involving the divergence and the curl 
of the velocity field, 


(13) 
(14) 


div v = q 
curl v = @ 


Taking the curl and the divergence of Eq. (10), we find 


(15) 
(16) 


(Ow/Ot) — vV*a = 0 
(0qg/Ot) — (4¥/3)V*g = —a°V*P 


The vorticity field obeys a separate equation not in- 
volving the other dependent variables, indicating that 
the vorticity field has no first-order interaction with 
the pressure and entropy (temperature) fields. 

Now we can show that the pressure also obeys a sep- 
arate equation. From Eq. (12) we can express the 
divergence of the velocity field, 


q = —(OP/dt) + (Os/Odt) (17) 


Substituting g in Eq. (16) and substituting the terms 
involving s from Eq. (11), we find 


4yv OV*P 7 


0 18) 
a” Of? 3a*_ Of 


The pressure obeys this equation, which is almost the 
wave equation. The extra term accounts for sound 
absorption at high frequencies. 
obeys the above equation will be called, in general, a 


A pressure field that 


sound wave. 

If we compare Eqs. (11) and (16), we find that both 
equations are of the heat equation type and that a 
quantity proportional to the Laplacean of the pressure 
acts as “‘heat source.” 

The general solution of both Eqs. (11) and (16) con- 
sists of two parts, one part depending on the pressure 
and another being the solution of the ‘shomogeneous”’ 
equations, 


(20) is not inde- 
pendent, since it must obey the conservation of mass, 
Eq. (12). The pressure-dependent solutions of Eqs, 
(11) and (16) can be found in the following way: We 
can express V’P from Eq. (18), which can be written 


as 
(: . 2 cep 10 = 2) ° 
~ = - - —_ ad 9 
” Or a@a\dt 3. (21 


where 
r = 4(y — 1)v/3a? 


The general solution of Eqs. (19) and 


7 can be regarded as the relaxation time of the medium, 
Substituting Eq. (21) into Eqs. (11) and (16), we obtain 
the special pressure-dependent solutions 


[1 + 7(0/Ot)]g = —OP/Ot (22) 
[1 + 7r(0/dt)]s = 7r(OP/OL) (23 


These equations indicate that the source field g and 
entropy field s follow the pressure waves not exactly 
in phase, but with a small time lag 7. This time lag is 
quite analogous to the thermal time lag of a hot-wire 
anemometer. However, the time lag 7 is extremely 
small for all practical cases. For atmospheric air it is 
of the order of 10~' sec. For many concrete examples, 
the operator 7(0/Of) is negligible compared to unity. 
Expressing the potential flow field with a velocity 


potential ¢ instead of g, we find, from Eq. (16), 


0 \ 0¢ fe) 
(1+, ) "=a )P (24 
or/ dr | Ot 


If the flow field is steady, 0¢/0/ = 0, and we find from 
Eq. (10) that P = 0. 
field. 
“sound waves” 


TF 
Y- 


This is a steady harmonic flow 
We shall take out this special case from the 
and add it to the ‘‘vorticity mode,’ 
since it does not contradict Eq. (15). 

Another special case occurs when V*P = 0. Then 
0°P/dt? = 0. From Eq. (16) we find g=0; therefore 
the field is harmonic. Eq. (10), we 
find that 


In this case, from 


P = —(1/a’)(0¢/Ot) (25 

We are now ready to define the three independent 
“modes.” 

(1) Vortictty Mode.—The principal variable is the 


The 
velocity field consists of two parts: One part is induced 


vorticity w. The governing equation is Eq. (15). 
by the vorticity found within the volume, and the other 
part is a steady harmonic field that is completely defined 
by specifying the time-averaged velocity at the bound- 
aries, 
l [ (E — &’) XK w(f’)dV(E’) 
i. 
for | 7 
, 


v(%) = 
: E— &') 


grad ¢; (26) 
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with O¢; Of = O and V’*¢g, = 0, P = 0, s = 0. The 
steady harmonic field fits naturally into this particular 
mode, since the vorticity field outside of the volume V 
induces exactly such a harmonic field within the volume 
as long as the diffusion does not progress so far that the 
vorticity outside of the volume will propagate into the 
volume. 

(2) Sound-Wave Mode. 
pressure, and it obeys Eq. (18). 
areq and s [see Eqs. (22) and (23)]. 


The principal variable is 
The other variables 
The special case 


1S 


oP/at = 0, V7P=0, g=0 


The flow is then harmonic, but since the case 0¢g/0/ = 0 
is taken out, the only solution is 


g = git, Vg’ = 0, Og.’/t = 0 


where g», the potential, varies proportional to time and 
y' is a harmonic function independent of time. 

The small entropy variation in the normally “‘isen- 
tropic’ sound field is caused by heat-conductive sound 
absorption, but it is a negligible quantity for most prac- 
tical cases. 

(3) Entropy Mode.—The compatible solutions of 
Eqs. (19) and (20) yield this mode. The principal vari- 
able is the entropy fluctuation s obeying Eq. (20). The 
decay of temperature (entropy) spots by heat conduc- 


tion changes the density and causes a motion. The 
motion is a potential flow 
v = (4v/3) grad s (27) 


or v = grad g3, where ¢g3; = 4vas/3. The motion as- 
sociated with the decay of entropy spottiness is rather 
small in all practical cases. 

We see from the results that two of the modes obey 
the heat equation, although they have different diffusion 
coefficients (in the ratio of 3:4), and the third mode 
obeys the wave equation with an extra absorption term. 
Since this is a linearized theory, the modes do not inter- 
act within the domain G. Interaction may take place, 
however, outside of the domain either at solid bound- 
aries or in regions where the fluctuations are not small 
and therefore nonlinear interaction is not negligible. 
If a point is outside of the domain G, it may mean one 
Either it is outside of the volume V or 
If in its past history 


of two things: 
prior to time f, or after time fp. 
the fluid went through a grid, the high-intensity turbu- 
lence produced strong sound waves and the viscous dis- 
sipation produced entropy spots. The former have 
propagated away according to the wave equation, but 
the latter have traveled with the fluid, diffusing slowly 
by heat conduction. When our observation starts (at 
t,), the modes are already sufficiently low in intensity 
so that nonlinear interaction between any two modes is 
negligible, although nonlinearity may have had a de- 
ciding role in the creation of these modes. 
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On the other hand, if these noninteracting modes 
pass through a region of a strong gradient, such as a 
shock wave, they may interact with the strong field and 
conversion from one mode to another may take place.** 


Second-Order Theory 


If the different modes are not sufficiently weak to 
warrant linearization and, on the other hand, are not 
strong enough to destroy completely the above picture 
of modes, a theory of weak interaction can be developed, 
keeping only the first- and second-order terms in the 
dependent variables. 

If the second-order terms are small compared with 
the they may be considered as 
“sources” 
Formally, we obtain the following equations: 


first-order terms, 
in the governing equations for the modes. 


(Ow/dt) — vV*2w = Q (a, P, 5) (28) 
Os 4 
Soles V2s = S(o, P, s) (29) 
ot 3 
1 O°P 4yv OV?P 
grees GIP a ee = II(w, P, s) (30) 
a’ Of? 3a* Ol 


where Q is the rate of vorticity creation by the action 
of all three modes, S is the rate of entropy spot creation 
by the interaction of all three modes, and II is the rate 
of sound generation by the interaction of all three 
modes. 

Since we assumed only second-order interaction, each 
term breaks down into six terms, each involving only 


one pair of modes. For example, 


IT,(@, @) + II,(P, P) + II3(s, S) + 
IT ,o(@, P) + IIo3(P, S) + IIs; (s, @®) 


Ii(e, P, s) = 


where II, is the source strength of sound waves gener- 
ated by the vorticity field, Il, is the nonlinear term in 
sound propagation (e.g., steepening of sound waves 
into shock waves), II; is the source strength of sound 
waves generated by the decay of temperature spotti 
ness, IT,. is the scattering of sound waves by the moving 
medium (vorticity field), IIp3 is the scattering of sound 
waves on the entropy spots (variation in speed of 
sound), and IIs; is the source strength of sound waves 
generated by the vorticity field hurling around tem 
perature spots. 

Similarly, S and Q each can be broken down into six 
interaction terms. The term Q,(@, @) is exactly the 
nonlinear term in incompressible turbulence theory 
or the 


accounting for the “‘stretching of vortices” 
“energy transfer term’’ in the spectral representation. 
There are 18 interaction terms in the second-order 
theory, and, depending on the specific problems, some 
can be more significant than others. 

Already a few attempts have been reported in the 
literature to estimate some of these interaction terms. 
Lighthill‘ and later Proudman” treated the sound 


generated by the vorticity (II, term), and Chandrasek- 
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har estimated the scattering of sound on the turbulent 
(vorticity mode) field.” Other attempts were made 
with a more statistical approach such as Moyal’s,' 


giving an explicit procedure to find the double correla- 
tion in coexisting isotropic turbulence (vorticity mode) 
and isotropic random sound waves. 


Further Simplification of the Modes 


Up to this point the size of the space-time domain G 
was completely unspecified. There are three charac- 
teristic dimensions entering into the problem. The 
first is the characteristic length of the medium and a 
related characteristic time 


(31) 


(32) 


lo = val 
ly = 


” 
v/a- 


These two quantities are naturally related to the 
mean free path of the gas; however, they are defined 
here on a continuum basis. 

The second important length parameter is the charac- 
teristic length of the disturbances. The simplest way 
to visualize this is to assume sinusoidal plane wave 
disturbances with wave length X. 

The third length is the dimension of the domain G. 


ble 


This can be denoted by L. 

The characteristic time interval can be defined as 
the divided by the speed of 
sound. 

If the time interval of observation is short and, in 


characteristic length 


the case of sound waves, the extent of the domain in 
space is correspondingly small and at the same time 
the characteristic length of the disturbance pattern \ 
is large compared with the characteristic length of the 
medium /), we find that the diffusion terms become 


negligible and we obtain, for the three modes:  vor- 
ticity mode, 
Ow/ Of = 0 (33) 
sound-wave mode, 
(1 /a?)(0?P/ot?) — V*?P = 0 (34) 
entropy mode, 
Os/dt = 0 (35) 


Eqs. (33), (34), and (35) simply state that the vorticity 
and entropy are ‘frozen patterns’ and that the sound 
waves obey the pure wave equation. 

Regarding the turbulence pattern as a ‘‘frozen pat- 
tern’ is not a new concept. It is a standard assump- 
hot-wire turbulence 
It is usually re- 


tion for evaluation of measure- 
ments in a low-speed wind tunnel. 
ferred to as For the low-speed 
If the 
turbulence level is low, it is simply assumed that the 


pattern does not change much during its passage in the 


“Taylor's hypothesis.”’ 
(incompressible) case, it was treated by Lin.’ 


neighborhood of the hot-wire probe; therefore, a hot- 


wire record in time at a fixed location x can be regarded 
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in first approximation as a space record (along the §, 
axis) of the turbulence pattern at an instant. 

The present result also indicates that Taylor's hy. 
pothesis may be applicable for two modes (vorticity 
and entropy) and that unattenuated wave propagation 
must be considered for the sound waves. 

A simple numerical example will indicate the orders 
of magnitude involved in this last simplification. — Per; 
odic patterns in space are possible solutions for all 
three modes: longitudinal waves for the sound-waye 
mode, transversal waves for the vorticity mode, and 
the Without the 
simplifications expressed by Eqs. (33), (34), and (35 


scalar waves for entropy mode. 
the waves decay in time at different rates for the dif 
ferent modes according to Eqs. (15), (18), and (20). 
For the numerical example we shall assume conditions 
prevailing in the supersonic wind tunnel that has been 


used for most of our experiments. Stagnation conditions 


were near atmospheric; pp) = 80 em. Hg; 7) = 300°K. 
Mach Speed of sound in working 
section a@ = 
0.5 cm.” per sec.; the mean velocity with respect 


Number = 1.75. 
28,500 cm. per sec.; kinematic viscosity 


y= 
to the wind tunnel (” = 49,800 cm. per sec. 
We take a periodic plane wave pattern with a wave 
length \ = 0.1 em. or wave number k = 63 cm. | 
(1) Vorticity Mode. 
nential decay 


Transversal waves with expo 


(m-é t/t » 
7 (50 


tk 
® = Gye 
(@o:n) = O- since div wow = 0 
n> = 1 


Substituting into Eq. (15), we obtain the decay time 


For our example, 4, = 5.0 X 107‘ sec. The vorticity 
pattern will travel in the wind tunnel during the decay 


Ut, = 25.2 cm. 


time a distance d,; = 
(2) Sound Waves. 
with a velocity and exponentially decaying in time 


Longitudinal waves propagating 


Psi: ™ (38 


The real part of Cis C, and gives the actual propagating 


velocity. The imaginary part gives the decay time 


Substituting into Eq. (18) we obtain 


2k . 2k 1 
ee ee 
od o oR 


2Ryv/ 31 1s 


Since in all practical cases the quantity 
small compared with unity, the propagation velocit) 
differs little from the nominal speed of sound a. The 
exponential decay is characterized by the decay time 


(40)| 


to = 3/yk?v 
For our example, % = 5.3 X 10-4 sec. The wave 
travels with the speed of sound during this time a dis: 
tance d: = ats. For our example, d; = 15 cm. This 
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distance is about 150 wave lengths. The relaxation 


time 7 1S 


T 4(y — 1)v/3a? 


For our example, r = 3.3 XK 10~" sec. The phase lag 
of the velocity field behind the pressure field is less 
than !/3 of 1 deg. The smallness of 7 also indicates 
that the entropy fluctuations caused by the decay of 
sound waves [Eq. (23)] are negligible compared to the 
isentropic fluctuation of temperature. 

(3) Entropy Mode.—A periodic pattern with expo- 
nential decay in time is a possible solution 


s = se" si (41) 
Substituting into Eq. (20) gives the decay time 
ls = 3 dv? (42) 


For our example, f; = 3.8 X 10~‘ sec. The flow pat- 
tern would travel a distance d; in the wind tunnel, d; = 
Ut. For our example, dj = 18.8 cm. The velocity 
field induced by the decay of the entropy pattern is 
periodic (longitudinal waves). The maximum velocity 


from Eq. (27) is 
[¥ |mar. = (4vk/3)50 (43) 


For our example, 


Vi maz./@ = 1.48 K 10-350 


This estimate indicates that the potential flow, induced 
by the entropy mode, is rather small. 

The numerical examples indicate that the viscous 
decay of the vorticity fields and the entropy field and 
the viscous attenuation of the sound field may be legiti- 
mately neglected if the wave lengths involved are not 
smaller than 0.1 em. and if the extent of the domain 
considered is not larger than a few centimeters. Natu 
rally there is no logical obstacle to reduce the size of the 
domain to the scale of the minimum size disturbance 
involved, and in this manner the argument for extreme 
simplification is strengthened. 

The analysis was carried out in a coordinate system 
moving with the mean flow. Our measuring instru- 
ments, however, are fixed with respect to the wind 
tunnel. We must later reinterpret these findings from 
the coordinate system of the measuring probe. 

(II) DETECTION OF THE FLUCTUATION MODES BY THE 
Hot-WIRE ANEMOMETER 


The Response of the Hot-Wire Anemometer 


The random fluctuations occurring in a supersonic 
wind tunnel can be conveniently detected and measured 
by the hot-wire anemometer. The principal advantage 
of the hot wire is not its inherently fast response but 
much more the fact that its response can be greatly 


improved by electronic methods.*!!'” 
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Naturally the hot-wire probe is fixed with respect to 
the tunnel walls and gives a signal proportional to fluc- 
tuations of flow variables at a point fixed in the wind- 
tunnel coordinate system (x, y, z).. On the other hand, 
the fluctuation modes introduced in the previous sec- 
tion were defined in a coordinate system moving with 
the average velocity (£1, &, &); therefore interpretation 
of measurements must take into account this differ- 
ence. 

The response of a hot wire placed normal to a super 
in some detail in 


sonic air stream has been treated 


references 8-12. The principal findings are the follow- 
ing: 

(A) The unheated wire assumes an equilibrium tem- 
perature 7, that is rather close to the stagnation or total 
temperature 7). 

(B) For moderate amounts of heating the heat loss 
of the wire is proportional to the temperature difference 
Tw-— 7,. It is also a linear function of the square root 
of the mass flow and is independent of the Mach Num 
ber as long as the flow is supersonic. 

(C) The hot wire heated by a constant electric current 
responds to a linear combination of mass flow fluctua- 
tions and stagnation temperature fluctuations. 

(D) The thermal lag of the hot wire placed in a super- 
sonic air stream does not essentially differ from the 
thermal lag observed at low speed, and the same type 
of electronic compensation can be used for both. 

The voltage fluctuations Ae(/) developed across the 
hot-wire probe (assuming proper thermal lag compensa- 
tion in the amplifier) can be expressed as a function of 
the two flow parameters: mass flow fluctuations m(?) 
and stagnation temperature fluctuations 3(/). 


A(pl’) Ap Al 
m = 100 = = 100 + 100 -= (44) 
pl Dp | 
8 = AT,/T, (45) 
Since 
Ty = T + (02/2C,) (46) 
d= a(AT/T) + B(AU (") (47) 
where 
a = 1/{1 + [(y — 1)/2]4P?} (48) 
B= aly — 1M? (49) 


The hot wire unfortunately does not respond directly 
to any one of the fluctuation modes derived in the first 
section but to two quantities that have no direct inter- 
pretation. From Eqs. (44) and (47) we see that the 
response has contributions from density fluctuations, 
velocity fluctuations (component in the direction of the 
main flow), and temperature fluctuations. The sensi- 
tivity to pressure does not enter directly. 

The output of the hot-wire anemometer is a function 


of time 
(50) 


Ae(t) = Aerd(t) Aeywm(t) 


where the sensitivity constants Sey, Aer characterize 
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Fic. 1. Fluctuation diagram for 1 per cent mass flow fluc- 
tuations and 1 per cent stagnation temperature fluctuations with 
varying degrees of correlation. 


the operation of the wire. Their values vary with the 
operating conditions of the wire, especially with the 
wire temperature. When the wire is barely heated, 
Aey — 0 and the wire responds mainly to stagnation 
temperature fluctuations. When the wire is well heated 
(a tungsten wire can be heated up to about 300°C. with- 
out danger of oxidation), the two sensitivity constants 
are about equal indicating that we cannot achieve 
mass flow fluctuation sensitivity alone. 

With a single wire we are unable to separate the dif- 
ferent fluctuations; however, if our interest is the 
measurement of the-mean square fluctuations, we can 
obtain those by taking repeated readings with the same 
wire at different operating conditions (i.e., different 
wire temperatures). 

The mean square fluctuation of the hot-wire output 


Ae? = Aer?3? + Aen?m? — 2AerAeymd (51) 


has three contributions: mean square stagnation tem- 
perature fluctuations 3*, mean square mass flow fluc- 
tuations m2, and the correlation m# (cross-correlation ) 
between mass flow and stagnation temperature fluc- 
tuations. Taking three readings at three different 
wire temperatures, we obtain three different sets of 
corresponding values of Ae?, Aey, and Aer, and we can 
solve three linear equations for the three unknown 
quantities 3?, m?, and md. 

Naturally we often find that the mean square voltage 
has experimental error and that we desire a more 
“graphic’”’ description of the situation. This desire 
has lead to the introduction of the ‘fluctuation dia- 
gram.” 

We define a ‘“‘virtual stagnation temperature fluctua- 
tion” t« 

Ae? 2 


v2 = -— = §? — 2rmd + r*m? (52) 
Aer" 


where r = Ae,,/ Aer is the sensitivity ratio. 
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Plotting J« versus r, we obtain the “fluctuation dia- 
gram.’ The sensitivity ratio is almost synonymous 
with wire temperature, and for tungsten wires the sen- 
sitivity ratio is almost equal to ay’, the overheating ratio 
of the wire. 


Gy’ = (Rw aa x.) R, 


(53) 


where R, is the resistance of the wire at equilibrium (un- 
heated) temperature and Ry is the resistance at actual 
temperature 7. For an unheated wire, ay’ = 0. 

Fig. 1 shows a fluctuation diagram for 1 per cent 
mass flow fluctuation and | per cent r.m.s. stagnation 
temperature fluctuation with varying degrees of posi- 
tive and negative correlation. 


m’ = Vn, 8 = V8, Rur = md/m'd' 


Although the fluctuation diagram enables us to deter- 
mine m’, 3’, and Ryr, it does not reveal how the mass 
flow fluctuations and stagnation temperature fluctua- 
tions should be attributed to the different fluctuation 
modes; however, a closer investigation reveals that all 
three fluctuation modes (vorticity, entropy, sound) 
have distinctly different fluctuation diagrams. 


Fluctuation Diagrams of the Different Modes 


If only one of the three fluctuation modes occurs in 
the supersonic flow, we can work out explicitly its 
fluctuation diagram. For our present purpose we 
shall use only the simplified forms characterized by 
Eqs. (33), (34), and (35). 

Vorticity Mode (Turbulent Velocity Fluctuations). 
The flow field is solenoidal and there are no pressure, 
density, or temperature fluctuations. The hot wire 
will respond to only AU’, the x component of the velocity 
fluctuation field. Introducing the 


velocity fluctuations in the direction of the main flow 


nondimensional 


u(t) = 100[AU (t)/U] (54 


The mass flow and stagnation temperature fluctuations 
become [Eqs. (44) and (47) ] 


m= u (55) 
ov = Bu 
The voltage output of the wire 


Ae/ Aer = (8B — r)u 
or 


ov. = IB—rin’ (58) 


The fluctuation diagram consists of two straight lines 
intersecting at r = 8 and having slopes of absolute 
magnitude u’. The response to velocity fluctuations 
becomes zero when 6 = r. 
ating point of the hot wire, since it is sensitive only to 
the and sound 
waves). 

Fig. 2 gives the fluctuation diagram for different 
Mach Numbers for solenoidal velocity fluctuations 


This is an important oper- 


two other modes (entropy mode 
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that have an r.m.s. value of 1 per cent of the mean 
velocity. 

Entropy Spottiness.—In the simplified form the en- 
tropy mode is another “frozen pattern’’ just as is the 
solenoidal velocity pattern (vorticity mode). We also 
neglect the potential velocity pattern set up by the 
decay of the temperature spots. There are no pressure 
fluctuations, but there are density fluctuations. The 
stagnation temperature 

The temperature fluc- 


mass flow fluctuations and 
fluctuations both are nonzero. 


tuation is defined as 


0 = 100(ATs, 7) = 100as (59) 
From Eqs. (44) and (47), 

v= 0 (60) 
m = O0/a = 100s (61) 

The voltage output of the hot wire becomes 
Ae/Aer = [1 + (r/a)JO (62) 

or 

dx = [1 + (r/a)]O’ (63) 


The fluctuation diagram shown in Fig. 3 is a straight 
line intercepting the y axis at r = —a. Since r cannot 
be made negative, no operating condition of the wire 
exists where the sensitivity for entropy spottiness 
vanishes. For wind-tunnel measurements it is more 
convenient to use, in place of s, the nondimensional 
temperature fluctuation 0 expressed in per cent of the 
total mean temperature than expressed relative to the 
free-stream temperature, especially because the total 
temperature changes relatively little in a wind tunnel 
even when the Mach Number varies through a wide 
range. 

This mode of fluctuation is somewhat 
Both the vor- 


ticity and entropy modes are ‘frozen patterns’’ altered 


Sound Waves. 
different from the previous two modes. 


only by slow diffusion, but the sound waves are propa- 
gating rapidly with respect to the fluid. The propaga- 
tion velocity is the speed of sound with respect to the 
fluid, but the apparent propagation velocity with re- 
spect to the tunnel depends on the orientation of the 
waves. Ina supersonic tunnel, sound waves can propa- 
gate both upstream and downstream and can even be 
stationary (Mach waves). 

The governing differential equation in the simplified 
form is Eq. (34). Introducing the velocity potential 


A 


Ap/p = P, v = grad ¢ 
P = (1/a?)(O¢/Ot) (64) 
AT/T = (y — 1)P (65) 


A plane wave solution propagating in the direction n 
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ness (entropy mode). 


P = P,(n-& — at) (66) 


where P,(n) is an arbitrary function of 7. The velocity 


field is 
0, = anyPo(mk, + nok + naks — at) (67a) 
V2 = ANeP (mb + Neko + sks — at) (67b) 
V3 = an3P (me + Mok, + Msk3 — at) (67¢) 


The mass flow and stagnation temperature fluctuations 
are functions of the orientation of the wave [Eqs. (44) 
and (47)] 
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space for plane sound waves. 


n 
m = 100 ( ~ = ') P,(n-& — at) (68) 
ve = aly = ee +- nM)P)(n-é — at) (69) 
When the sound wave becomes a Mach wave or 
n, = —1/M (70) 
the stagnation temperature fluctuation vanishes. This 


is a well-known fact obtained as a limiting case from 
an oblique stationary shock wave through which the 
enthalpy is unchanged. 

The voltage fluctuation across a hot wire is 


100 | Aera(y — 1)(1 + mM) — 


Ae : + i) Pi(t) (71) 
Mu M 0 ( 


As we see, there is no typical fluctuation diagram for 
However, taking the 


Ae(t) = 


sound waves of all orientations. 
finite frequency resolution of the hot-wire anemometer 
into account restricts the problem. 

The sound waves move with respect to the medium, 
and the medium is translated by the supersonic mean 
velocity. If we treat plane sinusoidal sound waves of 
wave length A, the sound-wave field can be expressed 
by the vector wave number 


k = (2nr/\)n (n? = 1) (72) 


P = P, exp t[ Rye + Roky a k3é3 _ akt | (73) 
Introducing the coordinate system fixed with respect to 
the tunnel (x, y, z), we find that the observed frequency 
f of the sinusoidal plane waves depends on the vector 
wave number. 


P = P, exp i[kix + key + kgs — 2nxft] (74) 
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where 


f = (ka + k,U)/20 (75) 
The observed frequency becomes zero when the wave 
becomes a Mach wave. The situation can be clearly 
seen in the wave number space given in Fig. 4. 
tor k represents a train of sinusoidal plane waves char- 
Dis- 
turbances that appear as fluctuations of a given fre- 


A vec- 
acterized by their wave length and orientation. 


quency fp to the fixed observer have wave number vec- 
tors forming two hyperboloids with asymptotic cones 
B’ and B” parallel to the cone of zero observed fre- 
quency A whose apex is at the origin. This particular 
cone is the Mach cone in wave number space formed by 
all wave number vectors representing a plane sinusoidal 
Mach wave train. 

The sound waves that are detectable by the hot-wire 
anemometer, which has a flat frequency response up to 


fo, must have wave number vectors between the cones 


of this constant observed frequency. If the sound- 
wave field is a homogeneous random field, it can be 
characterized by its three-dimensional power spectrum 
If the field is approximately isotropic, the energy dis- 
tribution is spherically symmetric in the wave number 
space. It is conceivable that the energy contribution 
from the small wave numbers (long wave length) is 
negligible compared with the contribution from the 
higher wave numbers. In such a case the most signifi- 
cant contribution to the hot-wire signal comes from the 
This annular volume in wave 
We 
can state, therefore, that the principal contribution to 


shaded area in Fig. 4. 
number space is close to the cone of Mach waves. 


the hot-wire signal in this case comes from waves that 
are almost Mach waves, provided that the contribution 
from long waves is negligible. Physically, this fact 
can be understood, since the waves that are close to 
being Mach waves move with a small velocity with re- 
spect to the hot wire; therefore, they can be easily 
resolved even when their wave length is short. 

The fact that intense sound waves steepen up into 
shock waves can be visualized in wave number space 
as a strong energy transfer toward the higher wave 
numbers. The stronger contribution in the high wave 
numbers, on the other hand, strengthens our previous 
arguinent. 

Shadow photographs from projectiles often show a 
sound-wave field emanating from the turbulent bound 
ary layer on the surface of the projectile, and these 
waves, although not exactly Mach waves, usually have 
angles close to the Mach waves. 

If the detectable sound waves can be assumed to be 
almost Mach waves, there are no stagnation tempera- 
ture fluctuations and the mass flow fluctuations alone 
Hence, the fluctua- 
tion diagram is a straight line from the origin. The 
strength of these sound waves can be characterized 
most conveniently by the mass flow fluctuations 


contribute to the value of Ae(f). 
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m =o = [1 — (1/M*)|P (76) 

or 
P = [M?/(M — 1)]o (77) 
Ae/Aer = —ro (78) 

or 
de = ro’ (79) 


The fluctuation diagram for sound waves that are the 
type described above is shown in Fig. 5. 

Naturally there are situations in which the sound 
waves of longer wave length with orientation different 
from Mach waves contribute significantly to the meas- 
ured r.m.s. fluctuation. In such a case the fluctuation 
diagram becomes different from the one shown in Fig. 
5, but the present approximation seems to suffice for 
sound generated by boundary layers, especially by 
transition bursts. 

Fig. 6 gives a summary of the fluctuation diagrams 
of all three modes at Mach Number = 1.75. 
Separation of Coexisting Modes 


In an actual flow in the wind tunnel we may expect 
to find all three modes with different intensities. They 
may also be either statistically independent or some- 
what correlated. Three fluctuating quantities at a 
point can be described statistically by the joint proba- 
bility distribution of these three variables. Restrict 
ing our interest only to second moments of this joint 
probability distribution, we find that there are six 
quantities to specify: the r.m.s. value of the three 
fluctuations and the three cross-correlation coeffi- 
cients. 

Since we saw that the fluctuation diagram specifies 
only three quantities, in general, it is not possible to 
determine all six quantities without some additional 
information or assumptions. 

There are, however, important cases when the de- 
termination of only three quantities is necessary. If 
the level of one of the three modes is negligible com- 
pared with the other two, there are only two r.m.s. 


























values and one cross-correlation to be determined. 
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Fic. 5. Fluctuation diagram for sound waves that are almost 


Mach waves having 1 per cent pressure fluctuations (Ap/p = 
0.01), 
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Fic. 6. Comparison of fluctuation diagram for three modes 
Such is usually the case when in a shear flow region both 
the vorticity mode and entropy mode are so intense 
compared with the sound waves that the latter can be 
neglected and the fluctuation diagram can be fully in- 
terpreted. Another important case occurs when all 
three modes can be assumed to be entirely uncorrelated. 
In this case only the three r.m.s. values need be deter- 
mined. Both of these special cases can be explicitly 
solved. 

The three fluctuating quantities that are observed 
at a point as a function of time are defined on a non- 
dimensional basis, using characteristic mean quantities 
referring to the supersonic flow as contrasted to the 
quantities w, P, s which were defined disregarding the 
mean flow. 

Entropy Mode. 
in per cent of the mean stagnation temperature. 


Temperature spottiness is expressed 


O(t) = 100[A7s(t)/To] = 100as 


Vorticitty Mode. -Solenoidal velocity fluctuation com 

ponent in the mean flow direction is expressed as per 
cent of the supersonic mean velocity. 
AU ,(t) 

u(t) = 100 Fa 


Detectable sound waves that 
Mass flow fluctuation is 


Sound-Wave Mode. 
are close to Mach waves. 


expressed as per cent of mean mass flow. 


( 2 
o(f) = {1 -—- 
M?*/ yp 


The fluctuation diagram of all three fluctuations oc 
curring together in the general case can be obtained 
by combining Eqs. (57), (62), and (78). 


0.* = k ts 0 + (8 — r)u — ro (SOa) 
a 
or 
J,* = (* = ‘0° + (Bp - aa + r°o*? + 9° = x 
a ai 
(8 — r)Ou — 2 ~ * " "80 —2(B—r)ruao (S0b) 
c 
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Fic. 7. Fluctuation diagram for 3 per cent temperature 
spottiness; 7.5 per cent anticorrelated and 5 per cent noncor 
related turbulent velocity fluctuations with sound field neglected; 
M = 1.58; a = 0.667; B = 0.667; 8’ = 3; wu’ =9; Ryo = —0.83 
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Fic. 8. Fluctuation diagram for uncorrelated modes at \J = 
1.75; temperature spottiness 0.1 per cent; turbulent velocity 


fluctuations 0.2 per cent; sound waves (detectable) 0.1 per cent 
of mass flow fluctuations. (Dotted lines show separate contribu 
tions. ) 


Fluctuating quantities can be decomposed into per- 
fectly correlated and uncorrelated parts with respect 
to another fluctuation. If one fluctuation is chosen 
as the primary variable, the second fluctuation can be 
decomposed as the sum of two parts, one perfectly cor 
related to the first and the other uncorrelated to the 
first. A third fluctuation can be decomposed into 
three parts, one part perfectly correlated to the first 
fluctuation, the second part perfectly correlated to the 
second part of the second fluctuation, and a third part 
uncorrelated to both 
Using these six fluctuation quantities we obtain de- 


first and second fluctuations. 
scription equivalent to the three r.m.s. values and three 
correlations. Using the 
first fluctuation 0, the velocity fluctuations can be de- 


temperature spottiness as 


composed into two parts: 


u= WM + Uo 


where “; = velocity fluctuation perfectly correlated 


to the temperature spottiness and “: = velocity fluctua- 


tion uncorrelated to the temperature spottiness; and 
og = 6, + oF + oz 
where o; = sound waves correlated to 0, o. = sound 


waves perfectly correlated to #, and 0; = sound waves 
uncorrelated to 0 and u. 
The relations defining these fluctuations are 


(Sla) 
(SIb) 


0(1)(@u/ 0") 
u(t) 


U;(1) —_ 


u(t) = — u(t) 
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a, = O(t)(Oc/0?2) (Sle) 
a, = U2(t)(uzo/u2?) (81d) 
03 = a(t) _— oi (t) “re o2(t) (Sle) 


Uy? + up? = u? (S1f) 
a1” + a2? + a3? = o? (Slg) 


The fluctuation diagram in terms of these six new 
variables is 


k +r 
a 


Special Case (I)—Sound Waves Are Negligible 


9 


0’ + (8B — r)u’ — ra | + 


x?(r) = 


[(8 — r)u.’ — roo’|? + [ro3’]? (82) 


2 6’ + (8B - al + [(8 — r)u'}? 
(83) 
The fluctuation diagram for the following values: 
1.58, 6’ = 3, uw; = 7.5 
Rie = —O0.83 


M = 


UW. = 09, u = 9, 


is shown in Fig. 7. These values are realistic and have 


been encountered in a turbulent supersonic boundary 














layer. 
Three important quantities can be read from the 
diagram 
I (B) the value 3x at r = 6B 
(d3«/dr), =, = the slope at r = B 
r=1f% = the intercept of the slope 
dx (To) = the value of d« at r = 7% | 


The three quantities can be obtained as 


8’ = ads(B)/(a 
(d3%/dr), ~3 — (0'/a) 
B«(70)/(B — ro) (S6 


+ B) (S4 
Uy’ = (5 | 


Uy’ == 


The value of 3«(7) is somewhat uncertain, since it Is 
often only inferred from the curve when 7) is nega 


tive. 


Special Case (II)—Uncorrelated Modes 


When the modes originate at different regions, it is 
possible that the three modes are uncorrelated. Such 
may be the case for the residual fluctuations encoun 
tered in a good wind tunnel. The fluctuation diagram 
is convex from below 


2 a =a r 2 7 = 
3%-(7) = 6? + (8B — r)*-u*> + rar (Si 


a 


By taking three different readings, the three unknowns 
can be found. Fig. 8 shows a fluctuation diagram ol 


this type. 





In actual hot-wire measurements these two special 
cases were utilized in the evaluation, since the situation 
found in the turbulent supersonic boundary layer 1s 
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close to Special Case (I) and the fluctuations found in 
the free stream are not unlike Special Case (II). 


(III) EXPERIMENTS 


Equipment 
The principal obstacles to the introduction of the 
hot-wire technique into the supersonic field were of a 
nature. Wires broke in the supersonic 
and no electronic equipment was available 


technical 
stream, 
which would give sufficient time resolution. 

The first measurements in supersonic flow were 
taken by the author with two purposes in mind. One 
was the verification of the previously expounded gen- 
eral ideas about fluctuation modes; the other was to 
perform exploratory measurements in turbulent flows 
in order to estimate orders of magnitude of fluctuations 
and the relative importance of different phenomena 
encountered. 

The wires used in the measurements at Johns Hop- 
kins were tungsten wires of 0.00015 in. diameter. They 
were electroplated with copper leaving a bare portion 
0.070 in. long and mounted by soft solder between 
needles set 0.110 in. apart. The copper-coated portion 
provided insensitive supports in the region behind 
the conical shock waves produced by the needle tips. 
It was found that many wires were lost because of 
fatigue failure of the copper sleeve, especially when the 
wires were operated in a highly turbulent shear region 
such as the turbulent boundary layer close to the wall. 
A drop of rubber cement at the soldered points provided 
satisfactory damping. Fig. 9 shows the arrangement. 
Such wires had long records of use, and one was oper- 
ated for 14 hours in the supersonic tunnel, accomplish- 
ing a ‘“‘mileage”’ 14,000 miles. We 
consider this record quite satisfactory. 

The electronic equipment is described briefly in 
The useful 


of approximately 


reference 9 and in detail in reference 11. 
frequency range of the combined hot wire and compen- 
sating amplifier is 2-70,000 cycles per sec. The noise 
level is correspondingly high. With the 
compensation level of 0.3-0.5 X 10~* sec., the mini- 
mum detectable fluctuation is about 0.2mV root mean 
square. 

The wire was overheated above equilibrium tempera- 
ture 7, by an amount varying from 15° to 300°C., and 
the sensitivity ratio r varied between 0.05-1.00. At 
Mach Number 1.75, the fluctuation diagram should be 
extended up to prose r = 0.75, but at lower Mach 
layer) a lower limit was 


time-lag 


Numbers (within boundary 
sufficient. 


Measurement of Pure Temperature Spottiness 


One of the earlier experiments was performed in 1951 
in the Bomb Tunnel of the Ballistic Research Labora- 
fluctuation 


tories, Aberdeen Proving Ground. The 


level was determined in the working section of the tun- 
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16.9. Hot-wire mounting for supersonic flow. 
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Fic. 10. Fluctuation diagram obtained in 1951 in the Aber- 
deen Bomb Tunnel; M = 1.73. Pure temperature spottiness, 


0’ = 0.3, explains results. 


nel at Mach Number 1.73. The above-described ideas 
about fluctuation ‘‘modes”’ 
and the measured data were filed away. 
it was suspected that the fluctuations were principally 
temperature spots, and their cause was determined 
later. The arrangement of the tunnel was such that 
the working section was supplied from four independent 
aftercoolers dis- 


were not yet crystallized, 
However, 


compressors, each having separate 
charging four air streams of different temperature into 
the supply section. It was conjectured then that, with 
the drifting of the relative temperature levels in the 
four units, the temperature differences could vary, and 
therefore the intensity of temperature spottiness would 
vary from run to run. At the conclusion of the meas- 
urements this condition was remedied, but later the en- 
tire tunnel was abandoned in favor of a newer design. 
It was rather interesting, however, to replot the 1951 
data in the light of present findings. Fig. 10 shows the 
fluctuation diagram with the straight line intercepting 


at r = —a corresponding to pure temperature spotti- 
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Fic. 11 Hot-wire detection of random shock-wave oscillations 
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Fic. 12 
wave; M = 1.75; 
cycles per sec 


Hot-wire signal obtained from oscillating weak shock 
A(pl')/pU = 0.04; calibrating sine wave 200 


ness (compare with Fig. 3). It is now believed that 
this one mode was so strong that it overcame the resid- 
ual fluctuations due to other modes present in the 


working section. 


Detection of Weak Waves 


Mach waves or weak shock waves that appear sta- 
tionary in a wind tunnel actually do oscillate either in a 
regular or in a random manner. If a hot wire is placed 
at the average position of such a wave, we can regard 
the flow field in two ways. From a coordinate system 
fixed to the wind tunnel we have a slightly moving, 
almost stationary, wave, and the hot wire is upstream 
for a portion of the time and downstream otherwise. 
From a coordinate system moving with the mean flow, 
these waves are almost plane sound waves with orien- 
tation close to the Mach angle. The fluctuation dia- 
gram should correspond to Fig. 5. 

A hot-wire probe was placed into a weak shock wave 
produced by the leading edge of a flat plate. Fig. 11 


OCTOBER, 1953 
shows the disposition. The hot-wire signal is shown jn 
Fig. 12. It appears as a random square wave, indj- 
cating that the hot wire oscillates in a random fashion 
between two distinct states. From mean flow (d.c.) 
hot-wire measurements taken before and after the 
shock wave, it was established that the mean mass flow 
change across the wave is about 4 per cent correspond- 
ing to a flow angle deflection of !/, deg. 

The fluctuation diagram is given in Fig. 13. The 
diagram indicates pure sound waves as in Fig. 6. The 
mass flow fluctuation indicated is 1.7 per cent r.m.s, 
If the mean mass flow change is 4 per cent across the 
wave, a pure square wave fluctuation would correspond 
to +2 per cent variation or to a 2 per cent r.m.s. value. 
On the otker hand, a sine wave type fluctuation with a 
! per cent peak-to-peak amplitude would give only a 
1.41 per cent r.m.s. fluctuation. Since the oscillo- 
graphic records indicated trapezoidal fluctuations, the 
1.7 per cent found from the fluctuation diagram is well 
bracketed by the 1.41 and 2 per cent values. The 
fluctuation diagram indicates that there are other fluc- 
tuations present (nonzero intercept) besides the domi- 
nating sound waves. 


Detection of Turbulent Regions by the Hot Wire 


Naturally the detection of transition in a supersonic 
boundary layer is of paramount interest. Generally 
speaking, it is not difficult to distinguish between the 
appearances of a hot-wire signal on the oscilloscope 
when the wire is in a true turbulent region or in a sound 
field emanating from a turbulent region. 

The author carried out some preliminary experiments 
in the 6- by 6-in. heat-transfer tunnel at Ames Aero 
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Fic. 13. Fluctuation diagram for random oscillating weak 
shock wave (arrangement in Fig. 11). 
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nautical Laboratory of the National Advisory Com- 
mittee for Aeronautics during the summer of 1952. 
An instrumented flat plate was explored with a hot- 
wire anemometer (the wire was a 0.0003 in. diameter 
tungsten about 0.200 in. long welded to steel needles). 
The arrangement is shown in Fig. 14. The suction slot 
was intended to remove the entire turbulent boundary 
layer from the tunnel wall before the flat plate. It 
was found that at low Reynolds Numbers (low pres- 
sure) the hot-wire still indicated intermittently turbu- 
lent bursts, and a strong vibration of the leading edge 
shock wave was also observed by the hot-wire probe. 


Detailed exploration showed the following state of 
affairs: The corrugated outer boundary of the turbu- 
lent boundary layer had crests that spilled above the 
leading edge, causing turbulent bursts in the otherwise 
laminar boundary layer. When the Reynolds Number 
was low, these bursts drifted downstream in the laminar 
boundary layer without much growth, making the 
boundary layer turbulent only a small fraction of the 
time. This picture also explained the generally laminar 
character of the boundary layer as detected from an 
averaging type of instrumentation (recovery factor 
measurement). However, at somewhat higher Reyn- 
olds Numbers, the boundary layer was found in a 
transition, with the turbulent bursts 
This 


slow state of 
gradually overtaking the entire boundary layer. 
explained some of the anomalous results obtained in the 
same boundary layer by other experimental means. 


Even an experiment of such a preliminary nature 
proved that hot-wire measurements of a qualitative 
character can greatly increase the understanding of 
peculiar phenomena, since it gives a detailed and in- 
stantaneous picture rather than merely long time-aver 


aged effects. 


Turbulence Measurements in the Supersonic Boundary 
Layer 
After the above-described preliminary experiments, 
the exploration of the supersonic turbulent boundary 


Department of Aero- 


layer was undertaken. The 
nautics at The Johns Hopkins University has a closed 
circuit continuously operating supersonic wind tunnel 
with a 7- by 1Il-in. working section. The tunnel is 
operated at a Mach Number of 1.75 and at atmospheric 
stagnation conditions (py) = 80 cm. Hg; JT) = 300°K.). 
The most conveniently available boundary layer lo- 
cated on the ceiling of the working section was chosen 
for the experiment. Its thickness was approximately 
0.500 in. Since no previous measurements were per- 
formed on turbulent fluctuations in a supersonic bound- 
ary layer, it was felt that the thickness was more at a 
premium than the controlled history of the layer. 


The probe mechanism permitted definition of a dis- 
tance about +0.002 in., which we considered satisfac- 
tory since only the exploration of the supersonic por- 
tion of the layer was planned. 





HOT -WIRE PROBE 





SUCTION 


Fic. 14. The turbulent boundary layer on the tunnel wall is 
not removed completely. Intermittent bursts are caused by 
peaks of turbulent boundary layer contaminating laminar bound- 


ary layer on flat plate (Ames, NACA, 1952) 
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Mach Number distribution across turbulent boundary 


Fic. 15 
layer 
Preliminary traverses indicated the following 
facts: The hot-wire signal is extremely characteristic 


when the outer edge of the layer is reached. The in 
termittent outer edge causes a strongly one-sided 
asymmetrical signal that gives the appearance of 
“growing grass’ as the probe plunges into the layer. 
The signal is rather high just as in low-speed turbulent 
boundary layers, but the fluctuations are rapid and the 
frequencies observed reach well into the 30-40 ke. 
range. 

The mean total temperature does not seem to vary 
across the layer within the accuracy of detection (about 
2°C.). Mean mass flow measurements, on the other 
hand, can easily be performed. With the assumption 
of constant total temperature (enthalpy) and constant 
pressure across the layer, the mean velocity and Mach 
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Fic. 16. Hot-wire signal in turbulent boundary laver; 1/ 
1.75. Vertical scale indicates approximately 10 per cent 
mass flow fluctuations 
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Number distribution can easily be determined. Fig. 
15 shows the mean Mach Number distribution across 
the layer. 

The experimental scatter of different sets of measure- 
ments is also indicated. The hot wire naturally is not 
the best instrument for mean flow measurements, since 
its advantage lies primarily in its ability to resolve 
rapid fluctuations, but it is still available for measuring 
mean properties with the degree of accuracy shown in 
Fig. 15. 

The character of the hot-wire signal at different dis- 
tances from the wall is shown in Fig. 16. The amplifi- 
cation is approximately the same in all records, and 
the marks on the amplitude scale indicate +10 per 
cent mass flow departure from the mean mass flow. 
The strong ‘‘single-sided’’ character of the signal close 
to the outer edge is visible. After the rough explora- 
tion the real question was to establish what modes of 


fluctuations were present. 
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Fluctuation diagrams were taken at different loca. 
tions within the boundary layer and just outside of jt. 
As a result of these measurements it was decided that 
the sound waves may be neglected within the turbulent 
boundary layer compared with the larger turbulent 
velocity fluctuations and temperature  spottiness, 
Since the mean flow has no discontinuity at the outer 
edge of the boundary layer, it is believed that the sound 
generated within the boundary layer propagates freely 
through the outer boundary and may only suffer some 
scattering, but essentially no appreciable reflection 
back into the layer. If this contention is true, the 
intensity of sound waves just outside of the boundary 
layer cannot be substantially different from its intensity 
just inside the layer unless there is some mechanism 
that would “‘trap’’ the sound waves within the bound- 
ary layer. The sound-wave level, expressed as go’ 
(1.m.s. mass flow fluctuation), certainly did not exceed 
0.2 per cent just outside of the boundary layer and most 
likely was only 0.1 per cent. On the other hand, the 
turbulent velocity fluctuation in the boundary layer 
was of the order of 2-3 per cent. Such being the case, 
the author believes that the effect of sound waves may 
be neglected in evaluating the fluctuation diagram from 
the signals within the boundary layer. This statement 
does not imply that the sound waves are negligible as a 
It only neglects their contribu- 
Sound 


physical phenomenon. 
tion to the overall mass flow fluctuation level. 
level is usually given in a decibel scale. Fluctuations 
in turbulence are given as per cent of mean flow proper- 
ties. The sound levels that are negligible for our pur 
poses may still be enormous in the ordinary sound in- 
tensity scales. This simplification, however, seems to 
be legitimate only in highly turbulent sheer regions 
where turbulence and temperature spottiness are so 


predominant. 


Using this assumption the fluctuation diagram could 
be interpreted unequivocally. It was found that both 
turbulent velocity fluctuations and temperature spot 
tiness are rather intense. It was also found that there 
is a high negative correlation (Rp, ~ —0.7) between 
the velocity fluctuation and temperature spottiness. 
This fact can be interpreted to mean that, whenever a 
lump of fluid is encountered which is warmer than the 
average, it comes from a layer closer to the wall and has 
also a momentum defect giving a negative velocity 
fluctuation compared with the mean value. Because 
of this high negative correlation between velocity and 
temperature, the character of the oscillograph traces 
did not vary much with changes in the overheating of 
the wire (or changing r), at least in the outer layers of 
the boundary layer. However, close to the wall, es- 
pecially at y = 0.050 in., the signal exhibited a different 
character by changing the wire’s sensitivity to different 
fluctuation modes. 


Fig. 17 shows the difference in traces at y = 0.050 
in. with sensitivity ratios r = 0.7 (marked ‘‘TEMP. 
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SPOTS”) and r = 0.05 (marked “MIXED"). The 
first case gives sensitivity to temperature spottiness 
alone; in the latter case there is about equal sensitivity 
to temperature spottiness (0) and to turbulent fluctua- 
There is a noticeable ‘‘spikiness’’ downward 
indicating cold intermittent This spikiness, 
however, disappears when the sensitivity is not pre- 
Since the oscillographic records 


tions (7). 
spots. 


dominantly for 90. 
give short samples from which to make statistical infer 
ences, it was felt that obtaining the probability density 
function of the signals would be more instructive. 
Probability densities were determined by taking long 
(1 min.) exposures on rather insensitive plates from the 
face of the cathode-ray scope. Reference signals gave 
the standard deviation, and the photographic process 
was standardized and calibrated. 

Photodensity measurements gave the probability 
All of the 


and a 
The 


density curves that are shown in Fig. IS. 


curves have the same standard deviation, 
Gaussian distribution is shown for comparison. 


signal is positive (wire hotter) toward the right. 


Two locations in the boundary layer (y = 0.050 in. 
and y = 0.350 in.) are compared; one is in the outside 
region where the intermittency of the outer boundary 
gives single-sided pulses toward the positive direction 
(the wire becomes hotter whenever the probe is within 
the hotter and slower fluid of the boundary layer as 


on 
\ / \ 
(d) AI 
/ 3 
/ 
/ 
/ 
/ | 
i 
/ 
Te al A. — = i —— A. 
| 2 374 3 24 0 1 2 3 498 


Probability density of hot-wire signals at two locations in boundary layers with two different sensitivity ratios 


The other location is 
With 


the wire sensitive to temperature spottiness, the single 


compared with the free stream 
at y = 0.050 in., where the situation is different. 


sided character is reversed, indicating that cool spots 
may cause a pulse structure. However, when the wire 
is sensitive to both modes, the probability density re- 


turns closer to Gaussian. 


The author must confess that this peculiar phenom- 
enon occurring close to the wall (vy = 0.050 in.) and 
also suggestively close to the sonic region of the layer 
(M = 1) is not satisfactorily explained as yet, and it is 
not even clear whether it is a real phenomenon of the 
turbulent supersonic boundary layer or merely an in 
strumentation phenomenon connected with operating 
a probe in the shear layer close to the Mach number 
unity. 


The distribution of fluctuation intensities is given in 
The r.m.s. temperature spottiness 0 is given 
If the tempera 


Fig. 19. 
in per cent of the total temperature. 
ture spottiness were given as per cent of the temperature 
difference across the layer (about 40 per cent of the 
mean total temperature), the two curves (u’ and 0’) 
would be almost identical, suggesting the not too sur 
prising fact that the momentum and heat transfer 
across the shear layer are similar; therefore the fluc 
tuation levels are also of the same order of magni 


tude. 








674 JOURNAL OF THE AERONAUTICAL SCIENCES 


% fluctuation 


} 








spectrum measured at 





== Ne _ 
lO inches 


Fic. 19. Temperature spottiness and turbulent velocity fluctua 


tions in turbulent boundary layer. 


Naturally the hot wire resolves only the large eddies 
in the layer. A measurement of the energy spectrum 
was performed in the range 500-70,000 cycles per sec. 


The results were plotted in Fig. 20 in the nondimen- 
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Fic. 20. Comparison of the energy spectrum of turbulent fluc 
tuations in supersonic boundary layer, y/6 = 0.5. 
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in reference 14. The intercept of the spectrum at zero 


wave number (frequency) determines the integral 


scale L,. When the nondimensional plot was com 
pared with the data obtained for low-speed boundary 
layer, essentially one parameter is free and can be ad 
justed. We found that the data from the supersonic 
layer at y = 0.250 in. matched surprisingly well with 
the choice L, = 0.110 in. This value is about a quarter 
of the boundary-layer thickness (0.500 in.) as contrasted 
to the findings in low-speed boundary layers'' where 
L, was found to be about 40 per cent of the boun lary 
layer thickness. This disagreement of a factor of 2 js 
not too surprising, since a low-speed boundary layer at 
moderate Reynolds Number has been compared with a 
supersonic boundary layer of much higher Reynolds 
Number. 


CONCLUSIONS 


(A) Random fluctuations in a supersonic flow can be 


decomposed into three distinct modes: vorticity, 


sound, entropy. The modes obey separate linear differ 
ential equations if the fields are weak and have moder 
ate interaction if the fields are stronger. 

(B) If the spatial scale of the disturbances is not too 
small, the effect of viscosity and heat transfer can be 
neglected for the short-time history of the disturbances, 
and the velocity, pressure, and temperature fields can 
be easily split into the three modes. The velocity 
field is split into a solenoidal and an irrotational field, 
the latter 
The pressure fluctuations 


The former constitutes the vorticity mode; 
belongs to the sound waves. 
belong to the sound-wave mode. The density and 
temperature fluctuations are split into isentropic fluc 
tuations that belong to the sound-wave mode and into 
a nonisentropic part that constitutes the entropy mode. 
density, or 


The vorticity mode has no 


temperature fluctuations. 


pressure, 


(C) The hot-wire anemometer can be operated suc- 
cessfully in a supersonic flow. The key to separation 
of the modes is the “fluctuation diagram’’ representing 
the r.m.s. output of a hot wire at different operating 
wire temperatures. 

(D) A variety of preliminary hot-wire measurements 
in supersonic wind tunnels indicate that all three 
modes are significant but that their relative importance 
may vary from one flow region to another. Sound 
waves may have more relative importance in the ‘‘free 
stream’ and the other two in shear layers (boundary 
layers, wakes). The 
spottiness) was always significant when the vorticity 
mode (turbulent velocity fluctuations) was intense. 


entropy mode _ (temperature 


(E) The fluctuation field in a supersonic turbulent 
boundary layer appears to be similar to a low-speed 
boundary layer. The only important difference is the 
intense temperature spottiness caused by the intense 
turbulent heat transfer across the layer. 


(Concluded on page 682) 

















It ha 
speed 
experin 
laver 
vestiga 
nomen: 


sudden 


A ti 
section 
series | 
end of 
larged 
changi 
The e3 
also vi 

The 
pressu 
ary-la) 
dimen 
range 
ind a 
tative! 
of inte 


irea ti 


Re 
"7 
the | 
fulfill 
Philc 
Ts 


ZETO 
gral 


‘Oln- 


lary 
ad 

onic 
vith 
irter 
Sted 
here 
ary 

2 is 
r at 
tha 
olds 
































The Effect of Boundary Laver on Sonic 
Flow Through an Abrupt Cross- 
Sectional Area Change 


ROBERT S. WICKi 
California Institute of Technology 


SUMMARY 


It has long been known that the pressure at the base of high 
Since 


experimental work on projectiles has indicated that the boundary 


speed projectiles is lower than the ambient pressure 


laver definitely affects the base pressure, it was decided to in 
vestigate the effect of the boundary layer on the similar phe 
nomena occurring in internal sonic or supersonic flow through a 


sudden enlargement of cross section. 


A two-dimensional flow channel consisting of a convergent 
section, a straight section, and an enlarged section connected in 
series Was constructed. The flow reached sonic velocity at the 
end of the throat section and expanded suddenly into the en 
larged section. The boundary layer in the throat was varied by 
changing the mass rate, pressure, and the length of the throat 
The expansion ratio (throat area to expanded section area) was 


also varied 


The results of these experiments indicated definitely that the 
pressure in the corner of the expansion was related to the bound 
ary-layer type and thickness upstream of the expansion. A 
dimensional analysis was used to correlate the results over a 
range of area expansion ratios and boundary-layer thicknesses, 
concept was devised which explained quali 


ind a “‘jet-pump”’ 


tatively the effects observed. Means of comparing the results 
of internal flow over corners and through sudden expansions in 


irea to external flow over corners were also considered 


NOMENCLATURE 


1,, A upstream and downstream internal cross-sectional 
areas 

B angle between are and expansion line (Fig. 2) 

H,, H height of upstream and downstream section of two 


dimensional test apparatus (Fig. 4) 


hi, length of upstream and downstream section in test 
apparatus (Fig. 4 

V Mach Number 

‘ag pressure 

R radius from corner to expansion line (Fig. 2) 

Re, Reynolds Number at entrance to expansion based on 


L, as length 


HW width of two-dimensional test apparatus 

A; rectangular coordinates (Fig. 2) 

a Mach angle 

6 angle isentropic flow turns through to reach Mach 
Number 17 

o direction of characteristics (Fig. 2) 

Y = ratio of specific heats 


Received May 25, 1953. 

* The experimental work described in this paper was done at 
the University of Illinois as part of a thesis submitted in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy 

t Senior Research Engineer, Jet Propulsion Laboratory. 


Subscripts 


I conditions upstream of expansion 

2 conditions downstream of expansion 

| conditions at exit of apparatus (atmospheric 

b conditions in corner between upstream and downstream 
0 stagnation conditions 


INTRODUCTION 


denen THERE IS A SUDDEN INCREASE in flow cross 
sectional area in either external or internal flow, 
a jet boundary will usually appear at the corner which 
separates the main flow from the wake in the corner. 
The 
expansion, whether for internal flow or for flow over a 


pressure in the corner downstream of such an 
projectile, is referred to as the base pressure. In 1776, 
Jean-Charles Borda! investigated the flow of water 
through a sudden increase of duct cross section, and the 
resultant Borda-Carnot 
give a fair estimate of the head loss that occurs in 
Nusselt* ap 


formula is still considered to 


sudden enlargements of cross section. 
pears to be one of the first to conduct experimental 
tests with high-velocity gas flow (subsonic flow of air 
and steam) through sudden increases in flow cross 
He later concluded from his more extensive 
that the 


section. 
experiments with sonic and supersonic flow 
base pressure will be equal to the entrance pressure if 
the entrance velocity is subsonic; but, if the entrance 
flow is supersonic, the base pressure could be equal to, 
The 


base pressure ratios (defined as the ratio of the pressure 


less than, or greater than the entrance pressure. 


at the base to the pressure upstream of the corner) 
that Nusselt obtained for supersonic flow in the entrance 
to the expansion ranged from 1.138 to 0.992 for a 4 
to 1 area increase and 1.294 to 0.695 for a 6.37 to 1 
area increase. In his experiments, Nusselt used both 
a convergent nozzle (J; = 0.970 to 0.981) with no 
straight (constant area) section and a Laval nozzle for 
the supersonic (.1/; < 2) cases. No attempt was made 
by Nusselt to determine the factors governing the base 
pressure, although a one-dimensional analysis based on 
the continuity, energy, and momentum equations was 
made. This, as will be shown later, is not usually 
satisfactory. 

Workers in the field of ballistics have long been 
concerned with the problem of sudden expansions of 
external compressible flow over the rear of projectiles 
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Fic. 1. Comparison of internal and external sudden ex- 


pansions. 


and its interrelationship with the base pressure be- 
cause the base drag, which is a considerable portion 
of the total drag, depends directly on the base pressure. 
It has long been known that the pressure at the base 
of high-speed projectiles is lower than the ambient 
pressure, and the manner in which most ballistics test 
data have been presented would lead one to the con- 
clusion that the base. pressure ratio is only a function 
of the flight Mach Number. 
tests, Chapman‘ and Kurzweg’ definitely showed that 
viscosity (i.e., boundary layer) effects have a definite 
effect upon the base pressure. Chapman also suggested 
the use of a criterion based upon a modified boundary- 
layer thickness to correlate data from recent tests on 
Since recent experimental 


In a series of postwar 


projectiles and missiles. 
work on projectiles has indicated that the boundary 
layer definitely affects the base pressure, it was de- 
cided to investigate the effects of the boundary layer 
on the similar phenomena occurring in internal sonic 
or supersonic flow through a sudden enlargement in 
cross section (cf. Fig. 1). The experimental study of 
this type of phenomena with an internal flow apparatus 
has a number of distinct advantages over usual ballis- 
tics test procedures. The size of the air supply needed 
is greatly reduced by eliminating the need for tunnels 
with large enough cross sections so that wall inter- 
ference, etc., will not disturb flow over the model. 
“Stings” and other support mechanisms can also be 
eliminated. The most important advantage of an 
internal flow apparatus is that complete static pressure 
and surface temperature measurements can be made 
not only along the entrance section to the expansion 
(analogous to the body of the projectile) but also in 
the wake region. These measurements are particularly 


valuable if one wants to test theoretical predictions 
adequately. 


COMPARISON OF INTERNAL AND EXTERNAL FLow 


There are a number of factors that prevent the 
direct comparison of base pressure phenomena present 
in internal two-dimensional flow with similar phenom- 
ena present in external axially symmetrical or two- 
dimensional flow. The first factor is that there are 
wall shear stresses acting on the wake and the re- 
attached boundary layer in the internal flow apparatus, 
and there are no equivalent shear stresses in the ex- 
ternal flow case. The second, and possibly more im- 
portant, factor is that internal flow upstream of the 
sudden expansion is invariably accompanied by severe 
pressure gradients if the approach section has parallel 
walls, while the pressure gradient along the cylindrical 
portion of a projectile is usually very small. The third 
factor is that the jet boundary will be intersected by 
‘Mach lines’? (expansion waves) originating at the 
opposite corner, which is not the case in external two- 
dimensional flow (cf. Fig. 1). The point of inter- 
section of the first expansion wave with the jet bound- 
ary represents the maximum expansion ratio possible 
for a given initial Mach Number and base pressure 
ratio for which the jet is not intersected by Mach 
lines from the opposite corner. This point can be 
calculated. 

Fig. 2 shows the expansion lines (characteristics of 
the first set) originating at one corner being inter- 
sected by the first expansion line (characteristic of 

















R 


Fic. 2. Interaction of expansion waves (characteristics ) 
from one corner with those from the opposite corner in 
internal two-dimensional flow. 
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EFFECT OF BOUNDARY 


the second set) from the opposite corner. The problem 
js to determine this intersecting characteristic. After 
this characteristic intersects the last expansion line 
from the corner, it is a straight line until it meets the 
jet boundary. The flow between the corner and this 
frst intersecting characteristic can be described by the 
Prandtl-Meyer solution. The system of coordinates 
for the present solution is rotated through the same 
entrance angle 6, sonic flow would turn through in 
order to reach the initial Mach Number 1/; under 
consideration. ¢; is the direction of the first expansion 
line, and its length is R;. Now consider the conditions 
at the point R. The angle between the characteristic 
of the first set (and also of the second set) and the 
streamline is a, the Mach angle, and is equal to arc 
sin 1/ AZ. 

Calling the angle between an r = constant line and 
the first characteristic of the second set B, the fol- 


lowing relation is obtained: 


—dR T 1 — tan’ a 
= tan | 2a — = (1) 
Rd@ 2 2 tan a 


The relationship between / and the angle ¢ is deter- 
mined by the Prandtl-Meyer solution, 


tan B = 


M? = (y + 1 — 2 cos? €)/(y — 1) cos’ « 


where € = Vy — 1)/(y + 1). 
With algebraic simplification Eq. (1) becomes 


dR l 
= — tan ede — ctn «de (2) 
R 3 Y¥+ 
and, when integrated between ¢€; to €« and R; to Rs, this 
becomes (for y = 1.4) 
R,/R. = (sin &/sin «) “(cos &/cos ¢) (3) 


R, is related to [/, in the following way: 
Il, = R, cos (of: — 6) 


@ and thus € depend on the base pressure ratio, and 
I, can be found from the intersection of the expansion 
line with the jet. The expansion line is straight after 
passing through the expansion lines from the opposite 
corner, and its slope in our system of coordinates will 
be (6. — 0, + az). 6 is the angle the flow would have 
turned through if it had started from sonic velocity, 
and a, is the Mach angle corresponding to the state 
determined by the base pressure ratio. 

If the entrance Mach Number is unity, no area of 
the jet will be unaffected by the expansion waves from 
the opposite corner. This feature alone makes the 
case of sonic flow through a sudden expansion unique. 
The factors listed indicate that the problem of internal 
flow is not strictly analogous to the corresponding 
problem in external flow, but there are enough over- 
all similarities to make the study of internal flow valu- 


able. 
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ANALYTICAL TREATMENT 


The behavior of compressible fluid flow depends 
that is, whether it is 


greatly upon the type of flow 
The corresponding flow 


subsonic, sonic, or supersonic. 
equations will be elliptic, parabolic, or hyperbolic, 
respectively. The type of equation is important since 
it determines how the boundary conditions enter into 
the problem. In the case of the parabolic and hyper- 
bolic type, the problem is an initial value type, and 
disturbances cannot travel upstream. In the case of 
the elliptic type, the problem is a boundary value type, 
and downstream conditions will affect the flow up- 
stream. 

In an internal flow problem such as the flow of a 
compressible fluid past a sudden increase in cross- 
sectional area, there could be more than one type of 
flow present. For instance, the flow in the smaller 
(entrance) section could be sonic or supersonic and the 
flow in the larger section subsonic or supersonic. In 
the case where the flow in the larger section is subsonic, 
the flow is determined upstream as far as the entrance 
if flow at the entrance is sonic or supersonic and past 
the entrance if the flow there is subsonic. A more 
interesting case arises when the flow in the larger sec- 
tion is supersonic and the flow in the entrance section 
is sonic or supersonic. In this case, the pressure in the 
corner is not determined by the downstream conditions 
but rather by the upstream conditions. The purpose 
of this investigation is to consider how changes in the 
upstream conditions affect the pressure in the corner. 

In any case, whether the flow in the entrance section 
is sonic or supersonic, there will be a region, because of 
viscous effects, in the neighborhood of the upstream 
surface where the flow is subsonic, and, therefore, the 
flow conditions (such as boundary-layer type, momen- 
tum thickness) in this region will then be interrelated 
with the conditions in the corner. These facts cannot 
be considered in a one-dimensional analysis, since the 
one-dimensional analysis (which requires that the 
momentum theorem be evaluated at the boundaries) 
does not go into the details of the internal mechanism 
of the flow. A two-dimensional theoretical treatment 
would make it possible to elaborate further on the 
details of the flow, but it still would not supply ade 
quate means of solving the problem of sudden expan- 
sions even when the flow is completely hyperbolic in 
nature, since it is necessary to know the base pressure 
before the two-dimensional flow pattern can be con- 
structed. 

A one-dimensional analysis will however be useful in 
establishing generally the range of physically possible 
solutions. If the one-dimensional energy, continuity, 
and momentum relations are applied between the en- 
trance to the expansion and some distance down 
stream where one-dimensional conditions might be 
approximately fulfilled, the following relations can be 


derived. 
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Fic. 3. Possible solution for internal expansion—2 to 1 
area ratio (result of one-dimensional analysis). 




















0 





Momentum: 
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(4) 
P, yr ey \ 
Energy and continuity: 
Pi\A, — Mz jl + [(y — 1)/2] Mo - 
PA, Mi \1 + [(y — 1)/2]0"4 ? 
Energy and continuity: 
Po = A\M, ( + [(y — 1), st + 1)/2(y — 1) / 
Po ry AsM, \1 + [(y — 1)/2]M;? (6) 


Combining Eqs. (4) and (5), P2/P; may be eliminated, 
and values of 4, and MM: which allow positive base 
pressure ratios can be determined for the particular 
expansion ratio A2/A, in question. A further physical 
limitation is that the stagnation pressure ratio given 


by Eq. (6) should be less than 1. The two foregoing 
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Schematic diagram of the two-dimensional test section 
g 





























Fic. 4. 


THE AERONAUTICAL SCIENCES 


OCTOBER, 1953 

physical limitations give both an upper and lower limit 
for the Mach Number 1/2 in the expanded section 
when the entrance Mach Number M, is specified. A 
further limitation results from the specification that g 
normal shock is the strongest that can occur at the 
entrance to the enlargement. The latter restriction 
gives a lower limit for A/. which lies within the limits 
defined by entropy limitations and thus defines the 
least value for M, for a given value of WW. The limita- 
tions are a function of the area ratio, and Fig. 3 is 4 
representation of the possible solutions for a 2 to | 
expansion ratio. 


EXPERIMENTAL PROCEDURE 


A two-dimensional flow channel consisting of a con- 
vergent section, a straight section, and an enlarged 
section connected in series (cf. Fig. 4) was constructed, 
The width of the flow channel was '/2 in., and the height 
(77,) at the throat was 0.381 in. These two values re- 
mained constant throughout the tests. 
amount of airflow available was '/2 Ib. per sec. at 110 
Ibs. per sq.in. gage stagnation pressure. The expansion 


The maximum 


ratio //,/H/, (throat area to expanded section area) was 
varied from 1.76 to 3.54 by varying the height //, of 
the expanded section. The flow reached sonic velocity 
at the end of the variable-length throat section (L, 
and expanded suddenly into the enlarged section. The 
boundary layer in the throat was varied by changing 
the mass rate, total pressure, and the length of the 
throat, allowing Reynolds Numbers at the throat based 
on the approach length to be varied from 1 to 15 times 
1. 
The 
was determined by wall pressure measurements, as 


behavior of the flow in the throat section 
well as by measured temperature recovery factors, and 
the base pressure ratio by measuring the pressure in 
the corner and at the entrance to the expansion. The 
test procedure was to vary the stagnation pressure 
over the available range for the 
tions of throat section lengths and enlarged section 


various combina- 
heights. 


EXPERIMENTAL RESULTS 


As shown in Fig. 4, the apparatus actually consists 
of two sudden enlargements in series. The conditions 
at the first (the sonic throat) are explicitly controlled, 
while those at the second (the flow from the expanded 
section into the surrounding atmosphere) are only im 
plicitly controlled. The expansion ratio for the second 
expansion is essentially infinite, and its base pressure 
i.e., P;. If the pressure 
in the expanded section were lower than the atmos- 


must be nearly atmospheric 


pheric pressure P;, the boundary layer in this section 
would be expected to thicken considerably and reverse 
flow might be expected. These conditions are ap- 
proached as the overall pressure ratio P,/ Pp is increased 


(i.e., the reservoir pressure Py) is reduced). Fig. 4, a 














plot | 
sure 
(Ly ] 
eter, 
ratio 
const 
the € 
ratio 
as th 
pitch 
chan, 
pand 
wall : 
the fi 
limit 
ther 
deviz 
more 
from 
appa 
only 
vest; 
resull 
abov' 
P,/P 
funct 
ratio. 
highe 
woul 
discu 
was | 
base 
expal 
curve 
the 1 
Reyn 
Fis 
lengt 
even 
was 
Plots 
direc’ 
value 
ratio 
of th 
ratio 


Th 
by of 
defin: 
thick 


T Fe 
equal 
sure 
throat 
critica 


limit 


ection 
d. A 
that a 
it the 
‘iction 
limits 
‘S the 
imita- 


” 


con- 
arged 
icted. 
eight 
eS re- 
mum 
t 110 
nsion 
) was 
Ts of 
aeity 
(L, 
The 
ging 

the 
ased 
imes 


tion 
» as 
and 
e in 
The 
sure 
ina- 
tion 


ists 
ons 
led, 
led 
m- 
ynd 
ure 
ire 
OS- 


ion 


ip- 


ed 











Seen 





EFFECT OF BOUNDARY 
plot of base pressure ratio P,/P,+ vs. the overall pres- 
sure ratio P;/Po for constant approach section length 
(L,/H;, = 11.81) with the expansion ratio as a param- 
eter, shows the effect of raising the overall pressure 
ratio Ps/Po. The base pressure ratio P,/P, remains 
constant until a definite value of Ps/P ) (dependent on 
the expansion ratio) is reached and the base pressure 
ratio in one corner starts to rise. It was noted that, 
as the reservoir pressure was lowered still further, the 
pitch of the sound emitted by the jet stream suddenly 
changed, and it was apparent that the jet in the ex- 
panded section was completely separated from one 
wall of the expanded section. This would indicate that 
the flow was symmetrical about the axis until an upper 
limit of P;/P» was reached; and then, as Py was fur- 
ther decreased, the base pressures in the two corners 
deviated as the jet in the expanded section became 
more unsymmetrical and finally separated completely 
from one wall. Because the mechanical design of the 
apparatus permitted measurement of the pressure in 
only one corner, the same type of phenomena was in- 
vestigated by means of the hydraulic analogy, and these 
results verified the existence of the unsymmetrical jet 
above definite values of P;/Po. The upper limit of 
P,/P) for the existence of a symmetrical jet was a 
function of the area ratio; the larger the expansion 
ratio, the lower the critical value of P,/P (i.e., the 
higher the critical value of Py below which the jet 
would become unsymmetrical). The remainder of the 
discussion will be restricted to the case where the jet 
was symmetrical or, in other words, cases where the 
base pressure was not affected by conditions in the 
Thus, only the flat portion of the 


expanded section. 
This restriction limits 


curves in Fig. 5 is considered. 
the range of variation of the flow variables such as 
Reynolds Number, etc. 

Fig. 5 shows that at a given area ratio and approach 
length (L,) the base pressure ratio remained constant 
even as the Reynolds Number in the approach section 
was varied by changing the reservoir pressure Pp. 
Plots obtained similar to Fig. 5 were shifted in the 
direction of higher base pressure ratios for smaller 
values of throat length. In Fig. 6 the base pressure 
ratio (symmetrical case) is plotted against the reciprocal 
of the expansion ratio //,///, with the approach length 
ratio L, //, as a parameter. 


ANALYSIS OF RESULTS 


Theoretical work by Crocco® and experimental work 
by others*® indicate that projectile base pressures are 
definitely related to the upstream boundary-layer type, 
thickness, and the location of the transition point. It 


t For the tests reported, the entrance Mach Number was 
equal to unity, and P; can be replaced by P*, the critical pres- 
Since isentropic flow is assumed in the core flow in the 
the 


sure 
throat section, P,/Po is equal to P,/P, times a constant 


critical pressure ratio. 
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Fic. 6. Effect of area ratio on base pressure ratio. 


would seem reasonable that the data obtained for the 
base pressures in internal flow could be correlated in a 
similar manner. Referring to Fig. 5, it is obvious that 
the Reynolds Number varies considerably over the 
flat or constant base pressure ratio portion of the 
curves. In fact, the Reynolds Number varies from 
4.6 XK 10° at the right-hand end of the curve for the 
1.76 expansion ratio to 10.6 & 10° at the left-hand end. 
The theories of Crocco and the available data on base 
pressures indicate that, if the boundary layer upstream 
of the expansion is not turbulent, the base pressure 
changes rapidly with Reynolds Number. The present 
data on base pressures for internal expansions would 
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= 1.76. 


indicate that the transition of the boundary layer from 
laminar to turbulent has occurred far upstream of the 
expansion. In fact, factors 
measured in the parallel throat section ranged from 
0.88 to 0.89. 

Since all that the 
boundary layer upstream of the expansion was turbu 
lent, Chapman's modified boundary-layer thickness 
ratio seemed to be a reasonable parameter for corre 
lating the data. Chapman reasoned that the thickness 
of the boundary layer just upstream of the corner de 


temperature recovery 


available evidence indicated 


termines the base pressure for projectiles and used a 
parameter that was the slenderness ratio (length to 
diameter of projectile) divided by the Reynolds Num- 
ber (based on length) to the one-fifth power (for turbu- 
lent boundary layers) to correlate his test results. If 
the present test results are put in this form using 
L,/f, as the slenderness ratio and L, as the length for 
evaluating the Reynolds Number at the sonic throat, 
a result similar to Fig. 7 is obtained for each expansion 
ratio (//://7,).. Chapman drew a curve similar to the 
dotted line in Fig. 7 through his data. This dotted 
line indicates that, asa first approximation, this method 
of correlation is applicable to internal flow results. 
However, it would seem that there is some justification 
for drawing separate curves through the data for each 
value of L,///,; as shown by the solid lines in Fig. 7. 
Similar curves were obtained for three other expansion 
ratios, but these did not coincide with Fig. 7 unless a 
factor was included which took into account the ex 
pansion ratio, as well as the approach length ratio. 
The following expression was obtained which expresses 
all the results as a function of the expansion ratio 
(/1,//I,) and the approach length ratio (L,///;) within 
the accuracy of the data (ef. Fig. 6): 


P, 
= {| 0.734 
P, ( ” 


0.55 


Dy ) 7 () 
exp 1.37 
L Hi, bE 


The following restrictions apply to Eq. (7): 


1.3 (the effect of the side walls on the 


(1) W/H, = 


two-dimensional aspects of the problem is not vet 


known). 
(2) 


18.37. 
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(3) P,s/Po < 0.05 + 0.1U07,/1,)T. 
(4) Transition point far upstream of the sudden 
expansion. 

Examination of Eq. (7) and Fig. 7 indicates that the 
boundary-layer thickness ratio L,/H, (Rez,)* is not 
sufficient for correlating the data, since the ratios 
L,/H], and /1,/ 11, are also necessary parameters. 
results, however, do not mean that the base pressure 


These 


ratio in internal flow is independent of boundary-layer 
thickness. The turbulent boundary-layer thickness js 
only directly proportional to the one-fifth power of th 
Reynolds Number in the case of incompressible flow 
with no pressure gradients. If there were no heat 
transfer, the fifth-power law might serve as an approxi 
mation even if the flow were high speed and compressi 
ble in the absence of pressure gradients. 
ditions are undoubtedly approximately satisfied along 
The essential difference 


(These con 


the cylinder of a projectile.) 
between the internal and external flow is that in in 
ternal flow the growth of the boundary-layer thickness 
has a much more pronounced effect on the pressure 
gradient. In fact, the boundary-layer growth and the 
pressure gradient in the constant-area approach section 
are closely interrelated by continuity considerations 
One must conclude that this interrelationship, together 
with the related shear stresses in the upstream section, 
acts in such a manner that the base pressure ratio is 
essentially constant in the range tested. 


A simplified ‘“‘jet-pump”’ type of theory can be de 
vised to explain qualitatively how the base pressure 
ratio depends on expansion ratio for a particular up 
stream geometrical and flow condition. The concept 
presented differs considerably from that presented by 
Hoerner.‘ Hoerner considers the boundary layer as 
an insulating layer that reduces the effectiveness of the 
jet as a pump, while the writer considers it as a source 
This is the essential difference 
The 
base corner is thought of as a sump with two supplies 
The first is the boundary layer flowing around 


of fluid for the corner. 
between his concept and the proposed concept. 


of mass. 
the corner and will be essentially constant if upstream 
conditions are constant. The second source is the back 
flow in the boundary layer along the wall of the ex 
panded section. This backflow occurs because of the 
pressure difference across the shock wave originating 
where the jet strikes the wall. This is schematically 
shown in Fig. 8. The amount of this flow will depend 
on the strength of the shock. 


the sump by the “ejector action’”’ of the jet and will, 


Mass 1s removed from 


among other factors, be proportional to the exposed 
area. Consider the sonic expansion pictured in Fig. 9, 
and let the line A-B represent the wall for the case of a 
finite expansion ratio. The expansion ratio is now in- 
creased to that represented by A’-B’. Since the slope 


of any sonic jet is greatest and is changing most rapidly 


t This restriction limits the values to the flat portion of the 


hase pressure ratio vs. P;/P») curves (cf. Fig. 5) 
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EFFECT OF BOUNDARY 


just past the entrance to the expansion, the resulting 
shock at B’ will be considerably weaker than that at B 
and less backflow will occur. The jet surface will in 
crease simultaneously, and the net result will be that 
more fluid will be pumped out in comparison to the 
Thus, the pressure in the sump will decrease 


influx. 
If the original 


yntil a point of equilibrium is reached. 
area ratio were that depicted by C-—D, an increase in 
area ratio to C’-D’ would also have caused a weaker 
shock at D’ to occur, but, since the jet slope is smaller 
in this region, the change in backflow would have been 
less than for the expansion change to A’—B’ from A~B. 
This indicates that, for the change of expansion from 
C-D to C’-D’, the decrease in base pressure would be 
less than for the change from A-B to A’—-B’. The two 
results of this explanation are: (a) As area ratio in 
creases for given upstream conditions, the base pres 
sure will decrease; and (b) the change in base pres 
sure for a given change in area ratio is much greater 
at small area ratios than at large area ratios. The 
first of these conclusions is confirmed by the experi 
mental results presented in Eq. (7) and Fig. 6. If 
Eq. (7) is differentiated with respect to //,///,, the 


Wl 
xp 1.37 (") 


(S) 


following result is obtained: 
(073 + 


which confirms the second result. 
‘jet-pump"’ concept can qualitatively predict the effect 
of the expansion ratio on the base pressure ratio for 


d(P,/P;) 0.55 ) 1.3 
d(Ho/H;) L,/H,) (e/M,)? 


Thus, a simplified 


internal sonic flow. 


CONCLUSIONS 

The study of sonic flow past a sudden increase in 
cross-sectional area reveals that: 

(1) The mechanism of internal and external flow is 
principally the same, and base pressure phenomena in 
external flow can be studied relatively easily by experi- 
ments with internal flow. 

(2) The essential problem is the determination of 
the base pressure when the flow in the enlarged section 
is supersonic. 

(3) The boundary layer upstream of the sudden 
expansion influences the base pressure ratio for a given 
expansion ratio. 

(4) An approximate correlation of data for any 
given expansion ratio can be established in the same 
manner as in external flow by plotting the dimension- 
less base pressure coefficient vs. a function of a modi 
fied boundary-layer thickness at the entrance. 


(5) A more accurate analysis of data indicates such 


a correlation involving the modified boundary-layer 


thickness does not account for the systematic deviation 


Irom a single curve (cf. Fig. 7). These deviations are 
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“Jet-pump” action in two-dimensional sudden enlarge 
ments 


Pia Ss 


Siol 
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jet with walls for finite two 


Interaction of sonic 
dimensional sudden enlargements 


apparently due to changes in the boundary-layer con 
figuration caused by the presence of pressure gradients 
in the entrance section. 

(6) A modified ‘‘jet-pump”’ 
which explains qualitatively the effect of area ratio on 


theory can be devised 


the base pressure ratio. 

As a result of this study it is thought that future 
work along the following lines would be profitable: 

(1) Extend two-dimensional studies to include 
supersonic flow at entrance. 

(2) Extend both sonic and supersonic two-dimen 
sional studies to different width ratios (W///,) in order 
to study the effect of secondary flow in side wall 
boundary layer. 

(3) Extend studies to include the effect of length 
ratio (L./II,) of the downstream section on jet sym 
metry. 

(4) Extend the concept of internal sudden expan 
sions to internal axially symmetric flow. This should 
have profound effects upon the symmetry of the re 
sultant jet at high P,/P, ratios, since the two corners 
now separated will be interconnected in the axially 
symmetric case. 

(5) Extend experimental work to cases where a 
laminar boundary layer is present upstream of the 
expansion. 


(Concluded on next page) 
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Subsonic Compressibility Corrections for 
Propellers and Helicopter Rotors 


E. V. LAITONE* ann LAWRENCE TALBOT? 
University of California at Berkeley 


SUMMARY 


Through the application of linearized theory a compressibility 
correction is obtained which gives the thrust or torque coefficient 
in subsonic flow, for a propeller or a helicopter rotor, directly 


from the incompressible flow values. 


The blade element theory is applied to subsonic flow by intro 
ducing the linearized compressibility correction, while the first 
order compressibility effect on the axial inflow is obtained from 
the simple momentum theory. The subsonic compressibility 
correction for the lightly loaded propeller is found to be expressed 
in terms of the forward flight Mach Number and the Mach 
Number based on the rotational speed of the propeller tip 
These relations are directly applicable to a helicopter rotor in 
vertical flight and reduce to a simpler form depending only upon 
the rotational Mach Number at the rotor tip. It is shown that 
for the helicopter rotor a satisfactory approximation, for tip 
Mach Numbers less than 0.8, may be obtained by using the 
simple Prandtl-Glauert with the Mach Number 


given by the rotational speed at the three-quarter radius. 


correction 


The derived compressibility corrections are found to be in 


agreement with available experimental data. The limitations 
of the relations are analyzed by comparison with a numerical 
solution calculated from a more exact procedure. 

A critical survey is made of other work on subsonic propeller 
theory. It is shown that in some cases fundamental mistakes 
have been made which result in considerable error in predicting 
the effect of compressibility 


. 
A correction for propellers or helicopter rotors. 
The problem seems to be too difficult even when 
the theory sub- 
The previous investigations have been 


INTRODUCTION 


PRESENT THERE IS no general compressibility 


considering vortex for linearized 
sonic flow. 
mainly on the basis of the momentum theory alone 
(see references 1 and 2), and practically nothing has 
been done on the blade element theory in subsonic flow. 
This is understandable since the propeller character- 
istics depend upon the blade shape and inflow distribu- 
tion, resulting in a complicated problem having no 
simple solution even in the incompressible case. 
However, the effect of compressibility may be closely 
approximated if we are only interested in the ratio of 
the propeller characteristics in subsonic flow as com- 
pared with those computed in detail for the same 
propeller by assuming incompressible flow. In this 
case we may easily obtain a general expression for this 
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ratio, at least fer lightly loaded propellers, since this 
ratio is not greatly affected by variations in the blade 
shape or inflow distribution. This is true because we 
actually compare the effect of compressibility on each 
section of the blade, element by element, and integrate 
the effect of the change in radius in order to get the 
principal effect of compressibility on the ratio itself. 
The influence of compressibility on the inflow dis- 
tribution is only a higher order effect that may be 
neglected, since the axial inflow itself only produces 
a second-order effect that may be neglected for lightly 
loaded propellers. 

The general compressibility correction, in the form of 
a ratio of the values in subsonic flow to the corre- 
sponding values from the incompressible solution, will 
be obtained by assuming the constant inflow given by 
the momentum theory for compressible flow. Then it 
will be shown that even the constant inflow produces 
only a second-order effect that may be neglected for 
lightly loaded propellers. This procedure is_prac- 
tically the same as introduced by Pistolesi (see pages 
245 and 248 of reference 3), who assumed constant 
values of the inflow velocities in order to estimate the 
effect of small changes in the solidity of propellers in 
the ordinary incompressible flow solution. His results 
showed that the actual varying inflow velocity distribu- 
tion had a negligible effect in predicting the change or 
ratio of the propeller characteristics, even with 50 
per cent changes in the solidity. 

In essentially the same way, assuming constant 
inflow, the ratio of the subsonic propeller thrust and 
torque, in relation to the incompressible value, will be 
evaluated by using the two-dimensional Prandtl- 
Glauert rule. The ratios found will be restricted by 
the limitations of linearized subsonic flow theory in 
not being able to predict the critical Mach Number 
where the lift of the blade element drops off from its 
linear relation. This effect can only be determined by 
experiment, although as a rule, for thin blades, it can 
be assumed to occur when the propeller tip helical 


Mach Number exceeds 0.8. 


List OF SYMBOLS 
Airflow Symbols 


V vp/p = local velocity of sound, ft. per sec 


a = 

p = local static pressure, Ibs. per sq.ft 
p = local density, Ib. sec.?/ft.* 

Y = 1.4 = specific heat ratio of air 
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P,= Po 
U y¥———y 
Uo(1 + 83) 








(3) 
Fic. 1. 

q = potto”/2 = free-stream dynamic pressure, lbs. per 
sq.ft. 

uy = axial advance velocity of propeller relative to free 
stream, ft. per sec. 

Wi = resultant local velocity relative to the propeller 
blade element, ft. per sec. 

w = propeller angular velocity, rad. per sec. 

M = u/a = local axial-flow Mach Number 

Ur = Rw/ay = rotational Mach Number of propeller 
tip 

r = uo/Rw = My/Mpr = propeller tip speed ratio 

MI = m/Ro, = A[(1 + 6)/(1 — 6’)] = actual rate of pro 
peller tip advance relative to air 

1+6 = u/u = axial interference velocity factor due to 
induced axial inflow 

| — 6’ = w/w = rotational interference factor due to slip 
stream rotation 

Subscript Symbols (See Fig. 1) 

0 = free-stream conditions far ahead of the propeller 

1 = local conditions induced at the propeller 

2 = local conditions immediately behind propeller 

3 = slipstream conditions an infinite distance downstream 


where static pressure is again free stream 


Propeller Symbols 


! = rR* = circular disc area swept by the propeller 
blades, sq.ft. 

B = number of blades 

ri = blade-section chord, ft. 

R = blade radius, ft. 

r = radial distance to blade section, ft 

0 = blade-section pitch angle, measured from plane of 
rotation to zero lift chord line, rad. 

a = (0 — ) = angle of attack of blade section measured 
from zero lift chord line, rad. 

o = tan! (m/ra,;) = tan~! (R/r)d\, = inflow angle in- 
duced at blade section, rad. 

n = Tu/P = Tu/wQ = propulsive efficiency 

z = propeller thrust, lbs 

Q = propeller-shaft torque, ft.Ibs. 

r = w? = propeller-shaft power, ft.lbs. per sec 

my = pirR*u, = assumed mass rate of flow handled by 


propeller disc as shown in Fig. 1 
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m = section lift-curve slope for two-dimensional incom 
pressible flow 

Cp) = drag (Ibs.)/goc = section profile drag coefficient 

tr = lift (Ibs.)/qoc = (dCx/da)a = section lift coefficiens 

Cr, Cg = standard NACA thrust and torque coefficients 


ANALYSIS OF LIGHTLY LOADED PROPELLERS 


From Fig. 1 the general expressions for the thrust 
and torque of the blade element are given by 


dT = B(1/2)p,W,?(C, cos ¢ — Cp sin ¢)c d? (1 
dQ = B(1/2)p,W,2(C_, sin @¢ + Cp cos $)rc dr (2 


where 
Wi? = uy? + (rw)? (3 
7 dC, m 
Cc. = a (0 — ¢) (4 
da V1 — (W/a)? 
op = tan™! (1;/rw)) (5 
In Eq. (4) the lift coefficient has been approximated 
by the linearized (Prandtl-Glauert) compressibility 


correction for two-dimensional subsonic flow, where 
m is the incompressible flow lift-curve slope (27 by 
linearized theory) and the angle of attack a must 
remain small. 

Now since Cp << C, (for Mach Numbers less than 
critical) and ¢ is small (except near the hub where the 
compressibility effects fortunately are negligible), we 
are justified in neglecting the term Cp sin ¢ in Eq. (1); 


therefore, 
ee be , 
T= S pi lds? t 
“ 0 
r m(@ — ) cos $ ; 
(ra)*| cdr (6 
V1 — (u%/a;)? — (re;/a;)? 


Even for uniform inflow Eq. (6) results in logarithmic 
expressions that become infinite at the lower limit for 
all distributions of 6 except for the particular case 


c(@ — ¢) cos ¢ = constant (7 


Since we are interested in a general expression for 
the ratio of the compressible thrust to the incom- 
pressible thrust for the same propeller, any distribution 
of c and 6 will be satisfactory as long as it is a reason- 
able one, and we can show that small variations in ( 
and @ have only a higher order effect upon the resulting 
ratio. Now, Eq. (7) is an appropriate distribution of 
c and @, since it corresponds to a constant chord pro- 
peller twisted so as to maintain a constant angle of 


attack as long as ¢ remains small-——that is, 
c = constant 
(9 — d) = ar = constant (S 
cos@ = 1 — ¢@7/2!+...=1 


If we define the propulsive efficiency following Glauert 


(see pages 204 and 232, reference 3) as 
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COMPRESSIBILITY 


uy dT a ( | ) ‘4 | tan @ | 
= 7 = Nes = ne 
w dQ 1 + 6, tan (@+ yo). 
(9) 
that is 
m = 1/(1 + 61) (10) 
m=1- 6’ (11) 
n3 = tan ¢/tan (¢ + Yo) (12) 
where 
tan yo = Cp/CL (13) 


then the constant angle of attack defined in Eq. (5S) 
corresponds to a maximum value of nz; or the minimum 
loss of thrust due to frictional drag on the propeller 
blades. The only approximation involved is that @ 
remains small, an assumption not justified near the 
hub. However, in this hub region the compressibility 
effects are negligible at forward speeds; 
therefore, the ratio of the compressible to the incom- 
pressible thrust is actually determined by the outboard 
sections of the propeller where Eq. (8) is valid. This 
is also the best choice for assuring the validity of 


ordinary 


Eq. (4), since the angle of attack remains a small 
constant along the entire blade. 

Substituting Eq. (8) into Eq. (6) and integrating for 
the case of uniform inflow, we obtain 


l YR 
, = = p;(Rw,)*Re m — x 


' + M2 . ( Mr, ) 
sin 
Mr,’ V1— Mp,’ 


V1 — Me? - Ma) 


Mr,’ 2 
(14) 
where 

Ra Rw(1 — 6’) Mr — 6’) - 
Mn, = — = (15) 

ay o(a;/ ao) Q;/ap 

uy uy (1 + 4) Mo(1 + 6) 
M, = = = (16) 

ay ao a,/ apg a;/ ag 


For the corresponding incompressible flow case we 
have pp = p,; anda— ™, so Eqs. (6) and (8) reduce to 


_ ] 2 ; QR ‘i ( uy < = 
I ie = 9 po( Rw)? Rc m 2 l+s3 Re (17) 


Taking the ratio of the compressible thrust to the 
incompressible value for the same propeller, we find 
that, when the induced inflow velocities are the same, 
then the thrust increases with either Mach Number 


by the ratio 
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1 + 3(M,2/Me,’) 
(18) 


Later, by using the momentum theory, we shall 


prove that, to the first order for a lightly loaded pro- 


peller, 
Pi1/ po = 1 — Mo? [( ig A ) 4qo| 
M, ( y-1 ‘)(" =) 
= l M,? +- 
Mo . 2 1qo ; (19) 
M — 1 T/A 
on Me ( ) 
Mr 2 4do 
and 
6) = (T/A) $90; 6’ 0 (20) 
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Consequently, for lightly loaded propellers, it is seen 
that Eq. (18) may be written to the first order as 


Cr r 3 
GG Tea % " 
l+M? ( Mr ) — Vi Mi? - MM 
Mr, V1 — Me? Mr? 
1 + 3(Mo?/ Mp?) 
(21) 
where 
Mr = Rw/ay (22) 


is the actual rotational Mach Number of the propeller 
tip, based on the free-stream speed of sound, and 


Mo = Uo/do (23) 


is the usual forward flight or free-stream Mach Number. 
Fig. 2 is a plot of Eq. (21) for constant values of J/p, 
while Fig. 3 shows the ratio, obtained for constant 
values of Mo, of the thrust at varying values of \/p to 
that when Mr = 0. In Fig. 2 the dotted lines shown 
for comparison are those given by the ordinary Prandtl- 
Glauert rule using the helical Mach Number occurring 
at the three-quarter radius so that 
= 1/V1—M,/? 


Cr Cine. . 


where 


M,,/2 = Mo? + [(3/4)Me}* 


A similar analysis of the torque expressed by Eq. (2) leads to the equivalent of Eq. (6) 
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It is evident from Fig. 2 that the Prandtl-Glauert 
rule cannot be applied directly to a propeller by the 
simple procedure of using the helical Mach Number at 
some fixed per cent of the blade radius, especially fo; 
larger values of either Mach Number. However, when 
Mr = 0, Eq. (21) exactly reduces to the simple Prandtl. 
Glauert rule as would be expected, 


(f), 
Cr VW 0 


inc. “"R 
(1 + Mo?)(1 — Mo?) 
2M, 





we (1 — M,?)' : : | 
V1— MM; 
(2la 
Also shown in Fig. 3 are the values obtained from | 
combining Eqs. (18), (19), and (20) for various values 
of the dise loading. 


It is obvious that Eq. (21) isa 
satisfactory approximation, even for relatively large 
dise loadings as long as AJ is small. 


It is interesting to note that the thrust remains finite 
at the limiting value of unity for the helical Mach 
Number at the tip (shown in Figs. 2 and 3), because 
the infinite lift predicted by the Prandtl-Glauert rule 
used in Eq. (6) only occurs in the infinitesimal area at 
the tip itself. 


Actually Eqs. (4) and (6) become invalid 
when the so-called critical Mach Number is exceeded, 
since the lift decreases and the drag abruptly increases. 
Consequently, only Eqs. (1) and (2) are applicable 
when the critical helical Mach Number (approximately 
0.8 for thin blades) is reached. 


namely, 


} 
m(@ — $) sing a | 





B {® 
OQ = - | pil,” + (ra,)?] | Cp cos@ + cr dr (26) | 
yA 0 (") (=) 
Ss \ ay ay 
Then if we take ¢ small, as in Eq. (8), and write, in addition, 
: u 1+ 6 
sing = @= P= = a ( ) (27 | 
Yow) ] = 6 
we can integrate Eq. (26) for ag = (@ — ¢) a constant 
) I I 11+ M,? Mr, V1— Mr? -— M?? 
ES ‘ - = Co ( = a) + Mar | — sin -7 -— = . (28) 
(1 2) pi(Rw,)?R (Rc) 4 2 a Mr,’ Vi1i- M,? Mr? 
where 
(e) LL + (m,) Jr 
rdr 
. i Ry, 
Do = z (29 
[(1/4) + (1/2)A,?7] 
Similarly, for the same propeller we have the incompressible flow solution, 
Cin C (, et l r ) rs At (1 4 3.2 0 
ae = i . ma 3A1") (30) 
(1 /2) po( Rw;)?,R?c : 4 2 ‘ 1 3 —_ a 


Cons 
govel 


Tl 
belo 
Cp 
trac’ 
fron 
pute 
(the 
add 
obte 
pres 
ider 
Figs 
S 
tor¢ 
Ma 
cier 
Eq. 
inc! 
par 
mai 
Ma 
the 


uy 
fic 


al 


slauert 
by th 
iber a 
tlly fo 
» Wher 
randt| 


- M,; 
(2la 


| from 
values 
) is a 
large 


finite 
Mach 
cause 
. Tule 
ea at 
valid 
ded, 
ASES, 
sable 


itely 


26) 


Wi 


8) 


e 
t 
T 
l 


























COMPRESSIBILITY CORRECTIONS 


Consequently, the ratio of the compressible torque to 


that is, for a lightly loaded propeller 


= 1 a A] 1 
Cp, 4 + 5 A?) + mar 5 


govern Eq. (21) 


C  Q 


C ™ Fine 7 a l 
~~ on c», ( ‘ 


inc, 


This ratio depends upon Cp,; however, if we remain 
below the critical Mach Number we may assume that 
Cp, is practically constant. Therefore, we may sub- 
tract the torque due to blade friction (the Cp, term) 
from the incompressible solution. Then we can com- 
pute the effect of compressibility on the remaining part 
(the induced drag torque or C,@ term). Finally, we 
add the incompressible torque due to blade friction and 
obtain the total torque in subsonic flow. The com- 
pressibility correction for the induced drag torque is 
identical to Eq. (21) and is therefore also given by 
Figs. 2 and 3. 

Since the thrust coefficient and the induced drag 
torque coefficient have the same ratio of increase with 
Mach Number, it is evident that the propulsive effi- 
ciency for a lightly loaded propeller, as defined by 
Eq. (9), will remain constant as the Mach Number 
increases to the critical value. This occurs because the 
part of the torque coefficient due to blade friction re- 
mains a small constant. However, when the critical 
Mach Number is reached, Cp, abruptly increases and 


the propulsive efficiency decreases. 


APPLICATION TO A HELICOPTER ROTOR IN VERTICAL 
CLIMB 


The application of these results to a helicopter rotor 
in vertical climb is especially simple since the disc 
loading is so light that not only is the inflow small but 
Mo and u,? may be considered zero. Consequently, it 
is easy to analyze the effect of any arbitrary blade 
twist in this case. 

Since M/) = 0 for a helicopter rotor, Eqs. (21) and 
(31) for a constant angle of attack rotor become 


Cr _ 3 (a a. ae ue’) (32) 


Cr. 2\ Mp Mr? 
Co, N ~ ‘Me V1- ~) 
— ma _ 
c. 4 *o\ Mr? Mr? i 
= ~ 00) 
CO: ne (Cp, 4) — Map| r 3) 
These correspond to the curves Mo = 0 in Figs. 2 and 


3 and are replotted in Fig. 4. Of course, they may be 
obtained directly by integrating Eqs. (6) and (26) with 
u,> = 0, pr = po, and a; = dy. Now with these simpli- 
fications we can finally integrate Eqs. (6) and (26) for 


any arbitrary blade twist. For example, the so-called 
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the incompressible value, for the same conditions that 


is 
+ M,° . 1 Mr V1 ving Mr’ — M,? 
Sil ; 
Me ~~ V1i—M;? Me? - 
; r Me? (31) 
A” } + mar (1 +3 : 
3 Mr? 
ideal blade twist defined by 
0 = Or(R/r) (34) 


which is the only one actually giving a uniform inflow 
distribution (see page 68, reference 4), gives 


(0 — ) 


Then, with cos ¢ = 


Cr/Cr., = 2 [a — V1 — M;?), Me | (36) 


= (Or — \)(R/r) = ar(R/r) (35) 


1, Eq. (6) yields 


Similarly, Eq. (26) gives 


Co, (’ - Yi- us’) 
. — Marr = 
Ce 4 Mr’? : 


rm (Cp,/4) — mar(d/2) (37) 


Fig. 4 shows that Eq. (36) agrees almost exactly with 
Eq. (24), the simple Prandtl-Glauert rule using the 
Mach Number occurring at three-quarters of the blade 
radius, at least for Mach Numbers less than 0.8. 
Eq. (36) also agrees fairly closely with Eq. (32) indi- 
cating that the effect of the blade angle distribution is 
It is obvious that, if we substitute 


(0 — o) = [0 — A(R/r)] 


slight. 


into Eq. (6), then the compressibility correction for any 
blade angle distribution, even the flat blade with no 
twist, will fall between the two values defined by 
Eqs. (32) and (36). Since Eq. (36) is for the ideal 
blade twist, which is the only one that actually pro- 
duces the uniform inflow distribution assumed in the 
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derivation, therefore it must be the most nearly correct. 
This indicates that Eq. (21) may slightly overestimate 
the effect of compressibility due to the impossibility of 
actually obtaining a uniform inflow distribution when 
the angle of attack is kept constant. However, this 
effect is slight, as may be shown by the fact that the 
Taylor series expansions of Eqs. (32) and (36) are 
nearly identical for the first three terms. 

Since the compressibility correction is limited by the 
critical Mach Number of 0.8, or less, it is evident from 
Fig. 4 that a satisfactory approximation for any heli- 
copter rotor in vertical climb is given by the simple 
Prandtl-Glauert rule using three-quarters of the tip 
Mach Number. A comparison with experimental data 
is also shown in Fig. 4. The dotted line is calculated 
from the tables given in reference 5 for a two-bladed 
tapered propeller having a 6 per cent profile thickness, 
a 2-ft. diameter, and a solidity of 0.038. The tests 
were made in a wind tunnel at forward speeds of 90 ft. 
per sec. or less, corresponding closely to My) = 0. The 
experimental data prove that the theoretical results 
are valid as long as the critical tip Mach Number of 
0.8 is not exceeded. 


THE MOMENTUM THEORY FOR SUBSONIC FLOW 


The subsonic inflow velocity distribution may be 
calculated by combining the blade-element theory 
with the general momentum theory. This is done by 
expressing the blade-element thrust increment as in 
Eq. (1) in terms of the unknown inflow velocity incre- 
ments at the propeller, station (1). This is then 
equated to the thrust increment determined by the 
velocity increments far downstream at station (;). 
This solution is relatively straightforward for incom- 
pressible flow, since the velocity increments are easily 
related to one another. For example, for light loadings 
the velocity increments at the propeller are one-half 
the magnitude of those far downstream. However, in 
subsonic flow this simple relation is further compll- 
cated by the density variation and the increase in en 
tropy due to any finite discontinuity or finite blade 
drag. 

To simplify this complicated problem we can use the 
one-dimensional axial-momentum theory, which yields 
only the average inflow, considered constant across the 
stream tube, and essentially replaces the propeller by 
the actuator disc. Then for the assumed 
stream-tube flow shown in Fig. 1 we 
Po = ps (see pages 330-334 of reference 6), 


so-called 
have, with 


) 


Uy) = MyUob3 (38) 
(39) 


T = m(u3 — 

m, = Apolo => A pil, = A spol» = A3p3lt3 
where A = rR? is the propeller disc area. Since the 
flow is isentropic from station (0) to station (1), 


ay" »s y¥-1!1 thie 
( ) a -|1 +7 Me (1 - a) 
ao Po < Uo” 


1/(y¥-1) 
Mo?(28, + 3 | 


Y 


(40) 
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Now the only known area in the assumed stream 
tube shown in Fig. | is A; = A = 7R*; therefore We 
can evaluate the mass rate of flow from Eqs. (39) and 


(40) as 


m, = pi,A = polyA (1 + 61) * 
1/(y—1) 


nel 7 
: Mob + 6:4 | (41 


1 — 
In addition, the total work done by the propeller 
may be evaluated from the one-dimensional energy 





equation by assuming no heat transfer from the stream 
tube to the surrounding air so that 
P dy” Uo” a3” U3” 
+ 


- eS (42 


mM, y— 1 . y— 1 2 


where P is the total power actually supplied to the 


a3" oo @ ? 
y—1 
U3" —— Uy" YPo | l 
my, -" i ( = (43) 
~ y= i Nee Po 


So far, all the relations derived are exact for the 
assumed conditions of an adiabatic process with uni- | 
form inflow as depicted in Fig. 1. Now a difficult 
problem arises in determining the relation between 4, 
and 6;. In reference 1 Krzywoblocki based his numeri- 
cal calculations on the assumption that p, = p), arguing 
that the pressure jump through the propeller disc 
happens too quickly for the density to change. How- 
ever, this is an impossible assumption for a perfect gas, 


propeller so that 


Tuo — Uy” 


U3" 
my, ; 
n é 


P = wQ = - 


since any pressure discontinuity in an adiabatic flow 
must be accompanied by a corresponding density dis- 
continuity. On the other hand, in reference 2 Vogeley 
based his calculations on the assumption that p; = py, 
the so-called isentropic propeller with no blade friction. 
That is, since the entropy change is given by 


S3 — Sp = C, In (2 )(2) ‘| = 
Wo] \ os) § 


¥7C, In (”) >0O (44) 
P3 


p3, the process becomes isentropic and 


S3 —_ Ry = 


therefore, if po = 
Eq. (43) reduces to exactly the same expressions ob- 
tained for incompressible flow, 


Ty l - 4 7 03 
P, = wQ; = = 5m | us” — uw") = Tui{1 + = 


m1 2 
(45) 


However, this is incompatible with the assumed thin 
discontinuity disc having A, = Aj, since the energy 
equation and the impulse-momentum equation can be 
combined to prove there must be an entropy increase 
with p; < po if a finite discontinuity exists between /» 
and p;. Consequently, if one assumes p3 = po, as in 
the calculations of reference 2, then one cannot tacitly 
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assume A; = A» and the existence of a thin actuator 
disc. However, a valid first-order approximation can 
be obtained by using the exact incompressible flow 
relations as the first step of an iterative process. Since 
Eq. (45) is exact for incompressible flow, it is seen that, 
if we also use the exact incompressible flow relation 


Fy = Tu => Tuo(1 + 6) (46) 


then we immediately obtain explicit relation between 


§, and 6; by combining Eqs. (45) and (46) so that 
l i 1 — 5 
6) =~ és = = ” (47) 
2 2m Uo m 


As shown in reference 3, Eqs. (45), (46), and (47) are 
valid for lightly loaded propellers in in- 


rigorously 
Therefore, we are now in a position 


compressible flow. 
to justify the use of all these equations, at least as a 
rigorous first-order approximation for lightly loaded 
propellers at low subsonic Mach Numbers. 

Since the relations are true only for lightly loaded 
propellers (for example, the rotation of the slipstream 
has been neglected, see page 194 of reference 3), we 
may further simplify the results by neglecting the 
higher powers of all the 6 increments. 


Then Eqs. (16), (38), (40), (41), (46), and (47) reduce 
to 
y, i — 
5, = CaN - ( 14) (1 — | - > (48) 
ee es m1 
r = 4qA5b, {1 + 6,(1 = M,*) | (49) 
PP, = = Tu(1 + 61) = Tu 1 (50) 
mm, = polio A [ l + 61 ( 1 — M,’)] (51 ) 
— | 
ee 1 — Mo75;); a (1 ot Mc) (52) 
pr ay 2 
M, : Se I ® “9 
M, — l + 6; ( +- 2 Me) (9.3) 


These provide the expressions used in Eqs. (19) and 


(20) for lightly loaded propellers. It is seen that 


6, ~ (1/A)/4qo (49) 


as for incompressible flow; consequently, the Mach 
Number effects on the inflow velocity are negligible 
when the disc loading is less than one-tenth the free- 
stream forward flight dynamic pressure (qo), see Fig. 3. 

Since the rotational velocity increment 6’ is always 
much less than 6,, except near the hub where the as- 
sumption of uniform inflow is invalid (see Fig. 46, 
page 3), it is obvious that 6’ = 0 is 
an extremely 
loaded propellers. 
found to be negligible; 
is found to be of the second order. 
the rotational increment 6’ in the general momentum 
theory shows that the main effect of slipstream rota- 
tion is in increasing 6, relative to 6; by the following 


242, of reference 
satisfactory approximation for lightly 
The effect of compressibility is 
actually the main effect of 6’ 
Consideration of 
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relation (see reference S): 


6; = (1 2)53[1 + 6,(M)? Mr’*) | (54) 
This does indicate, however, that Eq. (46) and the 
all that are really 


subsequent approximations are 
justified from the uniform axial-flow momentum theory, 
since the neglected slipstream rotation produces a 
second-order effect that must be considered, especially 
at large AJ, if a higher order effect is sought for the 
larger values of the disc loading. 


Fig. 5 compares Eqs. (21) and (24) with the values 
My is relatively small, 


given in reference 7. 
the difference between Eqs. (21) and (24) is slight. 


Since 


However, Eq. (21) seems to be in better agree- 
ment. 
A comparison of Eqs. (18) and (21) with all the 


experimental data in reference 7 shows that Eq. (21) 
agreement with the experimental 


is usually in better 
data are for 


data. This is not surprising, 
two-bladed propellers and the effect of a finite number 
of blades, which was entirely neglected in deriving 
Eq. (19), produces an influence that is opposite to that 
given by Eq. (19). This may be seen by using the 
linearized subsonic vortex theory to the first- 
order effect of an infinite number of blades for lightly 
power series (see reference 


since the 


find 


loaded propellers as a 
8), 
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I 6 
63 (1 = ? : ) (55) 
2 2V1 — M.? 


Another factor that indicates that Eq. (21) is a 
better approximation, even for finite loads, is the fact 
that Eq. (18) already slightly overestimates the com- 
pressibility because both a and 6, were assumed con- 
stant in its derivation. This may be seen by combining 
Eq. (1) with the differential form of Eq. (47), 

dT = pyt4,(26,u9)2ar dr (56) 
The numerical integration shows that the simultaneous 
assumption of constant a and 6; slightly overestimates 
the compressibility effect. 

Since the influence of a finite load, as shown in Fig. 3, 
is to increase the compressibility correction and Eq. 
(18) already overestimates the effect, therefore Eq. (21) 
is recommended for general use. The actual range of 
dise loadings that can be analyzed by using Eqs. (18) 
and (19) can be determined only by comparison with 
more extensive experimental data, since Eqs. (19), 
(54), and (55) provide second-order corrections that 
tend to cancel. Fig. 2 shows that the simple Prandtl- 
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Glauert correction using the helical Mach Number at 
the three-quarter radius, Eq. (24), provides a satis. 
factory approximation even for propellers as long as 
the Mach Number Mach 
Number remain less than 0.7. 


forward and the helical 
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A Steady-Flow Aeroelastic Study by 
Electrical Analogy 


ROBERT H. 


SCANLAN* 


Centre National des Recherches Scientifiques, France 


SUMMARY 

The present paper describes the use of the electrical tank ‘‘lift 
ing surface calculator” of Malavard and Duquenne in France for 
an aeroelastic study. The new items brought out here are: (1) 
the applicability of this electrical analogy to the calculation of 
certain ‘‘elementary’’ aerodynamic effects, in incompressible 
flow, which are of considerable utility; and (2) the specialization 
of structural influence coefficient concepts, permitted by the 
availability of the new means of aerodynamic calculation, which 
allow the calculator to avoid oversimplified notions such as that 
of an elastic axis and of the rigidity of the lifting surface in the 
wind or flight direction. 

With the analogy in mind, suggestions are also made for sub- 
sonic compressible flow aeroelastic lift distribution studies and 
for accounting for large dihedral effects on incompressible flow 


distributions 
INTRODUCTION 


oo AND DUQUENNE! HAVE DEVELOPED an 
effective electrical analogy for the determination 
of steady-state lift distributions on thin lifting surfaces 
of arbitrary plan form. In this analogy the electrical 
potential in a three-dimensional field is identified with 
the perturbation potential of the flow about the surface. 

Pines®® has used a matrix form of the lift distribution 
on more or less conventional wings to solve their steady- 
state aeroelastic problems. 

In the present study, the analogy of Malavard and 
Duquenne is fitted to the needs of an aeroelastic study 
of the type considered by Pines. Structural influence 
coefficients are so chosen as to be applicable to a large 
class of lifting surfaces, to avoid questions of an elastic 
axis, and to account for curvatures in the lifting surface 
in the direction of the wind. The electrical analogy 
proves to be sufficiently flexible to permit the deter- 
mination of new elementary aerodynamic effects di- 
rectly suited to the structural coefficients chosen. The 
resulting combination of aerodynamic and structural 
effects takes on a simpler form than in the original work 
of Pines. An example is worked out for a delta wing 
of aspect ratio 2. 
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tional des 


Other possibilities of the electrical analogy such as 
new uses of the elementary effects, study of the effect of 
effects at advanced Mach 


dihedral, and aeroelastic 


Numbers are mentioned. 


ELEcTRIC TANK ANALOGY OF MALAVARD AND 
DUQUENNE 


Malavard and Duquenne' have developed an elec- 
trical analogy for potential flow problems concerning 
thin lifting surfaces. They have exploited this analogy 
by means of their apparatus called the “‘lifting surface 
calculator.’’ It is worth while here to give a brief 
résumé of this apparatus and the electrical analogy 
involved. 

It is usual in electrical fluid-flow analogies to identify 
the hydrodynamic potential ¢ with the electrical poten- 
tial |’. Here this is done in a_three-dimensional 
(“deep”) tank of electrolyte of dimensions approxi- 
mately 5 by 3 by 4 ft. More precisely, the potential 
flow about a thin lifting surface as in Fig. 1 is prescribed 
by a potential 


g = Ux + (a, y, 2) (1) 


where ® is the perturbation potential about the surface 
in a free stream of velocity l’ far from the surface. It 
is found more practicable for the particular analogy at 
hand to consider only the perturbation potential &, and 
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it is accordingly this ® that is identified with the elec- 
trical potential |’ of the field within the tank. The 
tank is made of slate or other insulating material. The 
electrolyte found by Malavard to be most practical 
from a variety of points of view is ordinary tap water. 
Also, in practice, alternating current of low voltage and 
of frequency 1,000 cycles per sec. or over is used. 

The pressure difference 7 over the thin lifting surface 
of Fig. 1 (which is confounded with its projection in the 
xy-plane) is 


x = —pl(0®/Odx) (2) 


The condition that the airflow be tangent to the sur- 
face of Fig. 1 defined by z = f(x, y) is 


(v 4 ~*) of ‘ Od of 7 OP 
Ox / Ox Oy Oy Os (3) 
which simplifies to 
O*’/d0zs = L(Of/Odx) (4) 


for small slopes Of/ Ox, Of Oy and perturbation velocities 
OP/Ox, O'/Oy. The circulation I'(y) about a single 
section of the surface parallel to the x-axis is related to 
the potential & by 


r(y) = 28(x,, y) (5) 


where x, is the x-coordinate of the trailing edge at that 
section. 

The application of Ohm's Law to a current 7 passing 
into an electrolyte of potential distribution |’ through 
an electrode of area dS is 


1= —oa(dV/dn) dS (6) 


where o is the conductivity of the liquid and » is the 
outward normal to its surface. 

The analogy of Malavard and Duquenne consists of 
The tank is filled with 
The free surface of 


the following arrangement. 
water to the full depth possible. 
the water represents the plane of the surfaces of dis- 
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continuity which are to represent the wing, its wake. 
and the surrounding area. The wing, S, is represented 
by a flat surface, in contact with the free surface of the 
electrolyte, and covered by a series of closely spaced 
flush electrodes separated by insulating spaces (Fig. 2). 
Each of these small electrode surfaces, of area dS, has 
a controlled current passed into the electrolyte through 
it. The amount of this current, if dV/dn is identified 
in the analogy with 0®/0z, is given by Eq. (6) as 


i = —ol/(Of/dx) dS (7 


In the wake, the surface ¥ of Figs. 1 and 2 is in contact 
with the free surface of the water, and it consists of 
long thin electrode strips separated by narrow insulating 
bands. The surrounding area A in the same plane is 
placed at zero potential by means of a large metallic 
surface plate, at zero potential, just in contact with the 
water surface, and covering the rest of the free surface 
of the water in the tank. 
(7), are introduced and maintained, the trailing-edge 


As the correct currents, Eq. 


Kutta-Joukowsky condition is satisfied by placing each 
of the wake electrode bands at the same potential as 
the wing trailing-edge electrode just ahead of it. (To 
aid in accomplishing this simply in practice, each 
trailing-edge electrode is sometimes made to lie half 
in the wing and half in the wake. See Fig. 5.) 

The scale relation between © and V is given by 


®P = [U(of Ox) (2 a dS)|V (S 


or, if the “reduced circulation” y = 2/U/ is intro- 
duced (where / is the semispan of the wing being 
studied), 

, FR[(Of/Ox)dS} V 

= 2 " (9 
l Al 

where Al’ is a means of measuring the current 7 through 
use of a known resistance R: Ri = AV. 

The fact that the whole manipulation of this analogy 

is concerned with potential values (in particular, per 


turbation potential values) bears emphasis. For ex 
ample, from Eq. (2) lift values per unit span 
/ rdx = —pU® (10 


are available directly from the potentials to be read at 
the lifting surface electrodes. The studies normally 
made require only a semispan wing situated with its 
x-axis horizontal along one flat vertical wall of the tank 
and the y-axis in the free-water surface. Conditions 
of symmetry or antisymmetry of loading are repre 
sented, respectively, by making this wall insulating or 
conducting. 

A further exposé of this deep tank work of Malavard 
and Duquenne would be repetitious, but suffice it to 
point out that for wings of sweptback and other modern 
plan forms they have found excellent agreement with 
the lift distributions predicted by the analytical meth- 


ods of Falkner,’ Weissinger,* and others. The appara 
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tus described is in current use and is not limited to 
wings of regular plan form; the possibility of studying 
the lift distribution over wings of decidedly bizarre 
plan form is not the least of the advantages of this 
analogy. With modifications in technique even the ef- 
fects of tip tanks and end plates could be introduced 
into the studies. 

In the case of an untwisted rigid wing at an angle of 
attack a, the value a = —Of/ Ox is constant throughout 
the lifting surface, and, hence, by Eq. (7), the current 
iat each electrode is proportional only to the respective 
area dS served. Duquenne? has profited by this cir- 
cumstance to replace, in this case, the surface S (Fig. 2) 
by a simple hole in the plate A cut to the plan form of 
the wing. The current to replace the uniformly dis 
tributed electrode currents 7 is then furnished by a plate 
covering the entire bottom of the deep tank and raised 
to a potential above zero. The perturbation potential 
field of the wing is then read by means of a traveling 
probe dipped slightly into the water. The location of 
the probe in the horizontal plane is easily followed on 
an accompanying drafting board by a tracing point at- 
tached to an extension of the probe carriage. 

Both the wing plan form represented by a great num- 
ber of electrodes and that replaced by a cutout may be 
deformed geometrically as necessitated by potential 
flow studies at elevated Mach Numbers less than 1, so 
that the analogical system of Malavard and Duquenne 
may be said to be useful in the entire subsonic range of 


steady flow. 


ELEMENTARY EFFECTS 


As is commonly the case in structural studies, so also 
in aerodynamic studies on a new design it is desirable 
first to concentrate attention upon ‘‘elementary”’ ef- 
fects. For example, for a variety of loading conditions 
on a rigid aerodynamic surface it is convenient to have 
at hand a set of “unit’’ studies combinable at will. 
This is indispensable in aeroelastic studies, where ele- 
mentary structural effects (often given by a matrix of 
influence coefficients) interact with corresponding aero 
dynamic effects. It was then natural to require for the 
present study a flexible means for finding elementary 
aerodynamic effects applicable to a wide variety of lilt 
ing surfaces. This requirement is admirably fulfilled, 
in the author's opinion, by the lifting surface calculator 
of Malavard and Duquenne. 

One common structural assumption in wing studies 
(even of sweptback and delta wings) is that sections 
parallel to the wind are rigid. It was decided here to 
depart from this assumption, which is decidedly weak 
in certain cases, notably that of a flexible delta wing. 
The question of what type of structural influence coef 
ficients to use is immediately posed. Here it was de 


to demand certain “angles of attack per unit 


The actual wing under study (see ‘‘Example”’ 


cided 
load.”’ 
and Fig. 3) was divided up into a number (here 14) of 
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segments, assumed small enough so that only the forces 
(without residual moments) acting upon them needed 
to be taken into account at the centroids of the seg 
ments. The different segments in the chord direction 
being permitted to have different slopes, the following 
structural influence coefficients 6;; were sought: ‘‘the 
increase in angle of attack 6;; in the wind direction at 
the segment 7 due to a unit normal load at segment 7.’ 
The resulting square matrix of structural influence 
coefficients is here designated 0. It is to be noted that 
with these influence coefficients the subsequent analysis 
is rid of any “elastic axis.”’ 

The aerodynamic “influence coefficients’ required to 
correspond to the general structural ones defined above 
are purely ‘“‘mathematical’’ ones in the sense that, in 
dividually, they are isolated “unit effects’ scarcely ob 
tainable by an aerodynamic test on a physical model. 
However, the electrical analogy method provides them 
with ease, and here it has an advantage also over truly 
mathematical elementary effects, since these would be 
extremely difficult to calculate. Here the aerodynamic 
influence coeflicient (with reference, for example, to 
Fig. 3, where 7, 7 l, , 14) is “the lift on segment 7 
due to a unit change in the angle of attack of segment 


j.’ The resulting square matrix of aerodynamic in 


fluence coefficients is given by 
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[C/2)pl°P |G = g?’G (11) 


where the matrix G is composed of nondimensional ele- 
ments (Ay;,//) yi; 
7 due to a unit change in the angle of attack of segment 


“the reduced circulation on segment 


j multiplied by the ratio of the width Ay, of the segment 
7 to the semispan /.” 
air.) 


(Above, p is the density of the 
The technique of obtaining these aerodynamic 
coefficients will be explained in some detail in connec 
tion with the example presented. 


STEADY-FLOW AEROELASTIC THEORY 


Pines*®® has used elementary efiects in calculating 
wing loading distributions for aeroelastic studies on 
more or less conventional wings. The general form oi 
the work of Pines is well adapted to the present study. 
In fact, the present study takes an even simpler form 
for the reasons that the elastic axis is dispensed with and 
the elementary aerodynamic effects are drawn from 
simple analogical tests rather than from theory. 
the distribution of 
throughout a surface is considered. If the segments of 
the surface (as in Fig. 3) are numbered consecutively, 


Here aerodynamic loading 


the lifts p; at these segments may be represented by a 
column matrix {~} having as many items Pp; as there are 
segments 7. The aerodynamic lift forces are simply 
expressed by 


{Pp} = gl’G{a} (12) 


where 
fa} = far} + {ag} (13) 
is the angle of attack distribution of the segments, 
separated into the parts ag (due to rigid surface) and 
ag (due to elasticity of the surface). The elastic part 
depends on both the inertia loads p; and aerodynamic 
loads p: 
far} = O[{ pr} of tpt] (14) 


{nW}, W being the weight of the section 
Eqs. (12) to (14) com- 


where {~;} = 
subjected to  g’s acceleration. 
bined yield 


{p} = g?’Gl[ far} + O({pr} + {P})] 
= g’G\ar} + g@?GO{ pi} + g?GO}p} (15) 
or 
{Pp} = tp} + gA{p} (16) 


where /?GO = A. 

The directly calculable column {/}, is usually a good 
approximation to {p}. To include the 
part of Eq. (16), an iteration method then is easily set 
up by replacing {p} on the right by {p},, then by {p/p}, 
+ gA{p}i, ete. The result is the series 


{Pp} = (phi t+ gAlp}it @A*{p}it 
... FGA {p}hit+... 


aeroelastic 


(17) 


Because of the inherent nature of the matrix iteration 
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process, for some integer N large, A‘ {|p}, = \p}., 
where |/}, is a characteristic column associated with a 
characteristic value such that, for any integer ,, 


A'tpj. = Nt p}.. Thus, at the first point m such that 
A” |p}, = {p}., the powers of the matrix A in series 
(17) can be replaced by powers of X, 

(P} = thhi + GAt bhi + G@A* tp} + 


+ q™ipjc(l + gd + gr? + ) (18 


The series in parentheses sums to 1/(1 — gd) for gd <1, 
so that the divergence speed of the surface is given by 
q = I/X. 

Theory of this type is extendable to effectiveness of 
ailerons,° or other flaps and is, of course, simply appli- 
cable as it stands to the question of divergence. As 
presented it is suited only to the incompressible range. 
A word on the compressible range is given later in con- 
nection with a suggested iteration involving repeated 
applications of the electrical analogy. 


EXAMPLE 


A flat-plate delta wing of aspect ratio 2 is chosen for 
the present study (Fig. 3). Structural deflection in- 
fluence coefficients were obtained (in a somewhat crude 
fashion) from a 5-mm. thick Plexiglas model of semispan 
1.05 ft., and these were transformed by (also somewhat 
crude) calculation directly into the matrix 0 of slopes 
per unit loads. Greater care should be taken in ob- 
taining these coefficients in the case where results, 
rather than method, are important. It may be men- 
tioned in passing that another analogical method men- 
tioned by Malavard’ is useful in the approximate cal- 
culation of flat-plate deflection influence coefficients 
if it is desired to calculate rather than measure these 
coefficients. Williams* has suggested a semiempirical 
method for such calculation which requires the knowl- 
edge of the central flow line from a point of loading 


(source) to a symmetrically located point (sink) (on 
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the opposite half-wing) within a channel having the 
plan form of the wing. The flow lines (and their 
orthogonal family of equipotential lines) are rapidly 
given by an application of the “Teledeltos’’ conducting 
paper technique of Malavard.’ These flow lines were 
obtained for the present wing, for a single tip load 
only, by means sketched in Fig. 4, where thin copper 
electrodes are shown attached to the conducting paper 
and are placed at potentials of 100 and zero, respec- 
A few representative lines are included in the 
The conducting 


tively. 
sketch with their potentials marked. 
paper furnishes a plane field of uniform resistance, the 
electrical potential of which is measurable by a simple 
Wheatstone bridge. A _ pencil-like probe is pressed 
against the paper at various points to pick up the po- 
tential there and to measure it by balancing the bridge. 
The lines of equal potential are quickly plotted out. 
They are the analog of a fluid or stress-flow field. It is 
interesting that the Teledeltos technique rectifies di- 
rectly a weakness pointed out by Williams of the flow- 
line method—namely, the difficulty of making calcula- 
tions for a number of different load locations. 

The deflection influence coefficients of the Plexiglas 
plate were first calculated. The approximate slope 
influence coefficients calculated from them were next 
obtained (matrix 90). Because of the extremely small 
deflections at segments 1 to 4, their values and their 
effects elsewhere were taken as zero in this study. 
(Note that several tables, such as the above-mentioned 
influence coefficients and other lengthy numerical 
items are intentionally omitted from this paper in the 
interests of conserving space.) 

For the aerodynamics of the delta wing by the anal- 
ogy method, the model shown in Fig. 5 was constructed, 
having 110 electrodes composed simply of a conducting 
paint. Upon the figure are designated the relative 
values of the areas served by electrodes of typical con- 
figuration. The electrodes were grouped in such a 
way as to confine a distinct group of six or more elec- 
trodes within each of the 14 segments into which the 
lifting surface was divided. 

Fourteen separate ‘‘elementary”’ electrical analog 
tests were performed, in each of which currents propor- 
tional to the areas served by the electrodes of one single 
segment were made to flow, other segments not being 
fed electrically. Potentials were read at all the elec- 
trodes. Table | is a typical one (for segment 6) of the 
14 elementary distributions of potential, expressed in 
the reduced form of y per unit angle of attack (of seg- 
ment 6). 

As a means of verification of the techniques used in 
the elementary studies, the entire wing was fed elec- 
trically by currents proportional at each point to the 
area served. The result for the entire wing (represent- 
ing the potential distribution for a unit angle of attack) 
was compared with the sum of the elementary results. 
As a second means of verification, these two results 
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were compared with those obtained by performing the 
experiment with a cutout wing and current sent from 
the tank bottom. The results of the three different 
methods are presented in Table 2. It is seen that the 
sum of the elementary effects is high with respect to 
the results of the two other methods. Unfortunately, 
these results do not represent the accuracy obtainable 
by the methods, for certain manipulative errors were 
later found to explain the trend of these discrepancies. 
There may possibly be occasions where elementary 
effects defined by too-small groups of electrodes will 
falsify the results because of inadequate representation 
at a finite set of points of an essentially continuous 
phenomenon. Such occasions can always be remedied 
by changes of scale so as to increase the number of 
electrodes, etc. To improve elementary results, which 
are uniformly too high by a small percentage, it is also 
permissible to reduce all of them by the same factor 
that brings their sum into agreement with the usually 
excellent results of the cutout wing technique. Here 
it will simply be assumed that the elementary effects 
of the example are not in error. 


Eq. (2) shows that, to obtain the pressure difference 
distribution on the surface, each potential distribution 
must be differentiated with respect to x. On the other 
hand, Eq. (10) shows that lift per unit span is given 
directly by the potential; in fact the difference between 
potential values on either side of an elementary segment 
(in the x-direction) gives exactly the lift on that seg- 
ment when multiplied by the segment’s width. (The 
appropriate scale factor must be included when dealing 
with the values of electrical potential.) 


Proceeding thus, again in a somewhat crude manner 
and without graphing the potential distributions, the 
elementary lift distributions were obtained for unit 
angles of attack of the various segments. These were 
set up in a table of aerodynamic influence coefficients 
which represented the matrix G of reduced circulations 
yi, scaled by the width ratios Ay,// appropriate to their 
segments. The distribution of lifts per segment when 
the entire wing was studied as a unit is presented in the 
right-hand column of Table 3. The left-hand column 
represents G{l}. The disparity of the columns con- 
firms the crudity of the segment lift calculation—for 
which refinements are obvious. It was further assumed 
that all these lifts act at the area centroid of their respec- 
tive seginents. This is again a source of error permit- 
ting obvious refinements, such as the choice of a more 
probable (aerodynamic) centroid of loading under a 
representative loading condition. Also, if an increase 
in complexity of the analysis were permitted, introduc- 
tion of moment (as well as force) influence coefficients 
could account for more refinements. The subdivision 
of the wing into many small segments has the effect of 
cutting down the seriousness of the errors of geometry 
retained here. With automatic computing-machine 
aids many more than 14 segments could be employed. 
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The matrix GO was next calculated. ime of Ll (oy Of Oy O/ . 
For the calculation of this example the lift per unit ~ Ox Ox Ox = Oy ~Oy (20) 


angle of attack is given by 


Ip) = 9G + AGI) + PARGHLY + 

GARG} 1} + @AIGI{IL[1/C — gr)] 
since the iteration process converged on the fourth 
multiplication. Above, and for the example only, the 


following notation is used: 


Items of interest in the iteration of A are given in 


Table 4. 
‘rad. per Gm. = 1.79 X 10-* rad. per Ib. 


A = 398 X 10 


For the critical divergence velocity lL’, at sea level (p = 
po) 


(/2)pU PA = 1 
U, = 445 m.p.h. 


THE EFFECT OF LIFTING SURFACE DIHEDRAL 


It has been mentioned above that account can be 
taken of symmetry or antisymmetry of the wing sys- 
tem. It is, in the symmetric case, possible to take 
account of wing dihedral and the effect of spanwise 
flow on the overall results. This is a nonlinear effect 
as is seen from the form (3), where 
As an example, consider one of 
From Eqs. 


‘“‘small’’ terms are 
no longer neglected. 
the elementary effects mentioned above. 
(3) and (9) there is obtained 
2a0RdS V 

IAV 


7/0 = (19) 


where @ is the ‘‘generalized angle of attack,”’ 


of the segment in question. It is the distribution y/@ 
which is obtained directly from the electrical analogy. 


If Eq. (20) is written 


i an of (7 6) Of ; O(y/0) OF a) 
~ Ox = 2\ dx or dy dy 
it follows immediately that 
-_ Of /Ox 
- 1 (O(y/0) Of — A-y/8) "7 (21) 
2 ( Ox Ox Oy oy 


Thus @ is seen as a generalization of the angle of attack 


TABLE 3 TABLE 4 TABLE 5 
Effective Per Cent 
Drop in the 
Actual Angle of 
Attack of 20° 
Caused by a 
Dihedral of 20° 
at the Elementary 
Segment 
(Symmetry of 
Segment with 
By 7 . Result Opposite Half- 
102 Table |10? X Gilt of Iter Wing Is 
G}1} 2 19.7 ation Assumed ) 
l 2.4 2.1 0.124 0.0354 9.5 
2 3.3 3.2 0.169 0.0850 5.4 
3 6.4 6.2 0.326 0.1773 6.5 
} 1.3 10.8 0.563 0.1785 8.5 
5 2.4 2.4 0.122 0.0291 9.0 
6 3.5 2.2 0.180 0.0405 2.5 
7 6.4 6.4 0.325 0.1573 3.3 
8 7.4 16.8 0.885 0.4948 6.1 
9 2.5 2.4 0.129 0.0638 5.9 
10 4.7 5.0 0.240 0.1381 3.0 
11 19.3 18.6 0.981 0.7079 6.4 
12 2.6 2.8 0.133 0.1523 8.2 
13 19.7 19.2 1.000 1.0000 me 
14 15.0 14.5 0.762 0.9137 17.9 
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Of/Ox for the case where the elementary segment is 
canted both in the wind direction (Of/Ox) and in the 
spanwise (dihedral) direction (Of/Oy). Table 5 pre- 
sents some results for elementary effects showing 
the per cent reduction in effective angle of attack: 
i [(Of/Ox) — @]/(Of/Ox)} X 100 per cent for an assumed 
angle of attack of 20° and a dihedral angle of 20° of the 
elementary segment. Such studies might be of in- 
terest in certain V-empennage or star-fin arrange- 
ments where dihedral is large. 

If overall, instead of elementary, effects are dealt 
with, it is not possible to have y directly from a single 
experiment, for, by Eq. (20), y is given by a differential 
equation the coefficients of which vary over the wing 
area. An iteration for y can be envisaged simply as a 
preliminary distribution of y (unmodified), the slopes 
of which are to be used in Eq. (20) to reestimate 8, 
followed by a second experiment using new currents 


i= —aV6dS (22) 


and yielding a new approximation to y, etc. 


CONCLUSION 


The centering of attention upon elementary effects, 
as in the present study, appears feasible from the anal- 
ogy point of view and effective in a typical aeroelastic 
problem. In addition, it points to the possibility, aside 
from questions of aeroelasticity, of taking into account 
dihedral; of investigating the efficient location of ailer- 
ons, elevons, or flaps; of the orienting of the hinge 
lines of these elements so as to benefit from regions of 
decided spanwise flow; of the studying of various rigid 
wing twist distributions as sums of elementary dis- 
tributions with the aim of achieving a given surface 
lift distribution (with or without elastic changes). 

The possibility of compressible flow studies has been 
mentioned. For the case of aeroelastic effects, matrix 
iterations based on linear effects are no longer possible; 
but it remains possible to make such investigations on 
the ‘‘analog-iteration’’ basis, somewhat as indicated 


for the case of dihedral. The procedure is, for example, 
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to study a given “rigid wing’ loading at a given Mach 
Number. With the loading found, a new wing defor. 
mation is calculated from the structural data. With 
recourse again to the analogy, a new loading at the 
given Mach Number is found. The process may be 
repeated until convergence at the given Mach Number 
and likewise may be reproduced at several Mach Num 
bers. Actual speeds of reversal of control effectiveness 
are then obtained from extrapolation of the results, ete 
It is generally more useful to have at hand a continuous 
plot of ‘effectiveness’ data than to know simply the 
critical velocities of reversal or divergence. The itera- 
tive use of the electrical analogy in compressible flow 
cases, as suggested here, resembles some of the old 
incompressible flow methods in use prior to the advent 
of recent aeroelastic theory: such methods habitually 
calculated loading, then deflection, then new loading, 
etc. However, the analogy improves the situation by 
doing the difficult calculation—the aerodynamic load 
distribution—rapidly. 
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The Second-Order Lifting Pressure and 
Damping in Roll of Sweptback Rolling 
Airfoils at Supersonic Speeds 


JOHN C. MARTIN* ann NATHAN GERBER? 


Ballistics Research Laboratories, Aberdeen Proving Ground 


SUMMARY 


The lifting pressure distribution and the damping in roll are 
determined for rolling sweptback airfoils with arbitrary cross 


sections. The time variation in the partial differential equation 
is eliminated by utilizing a set of coordinate axes that are attached 
to and rolling with the airfoil. The second-order lifting pressure 
is determined by an iterative method, and the damping in roll is 


calculated from the pressure distribution. The results are shown 


to agree well with available experimental data. 


INTRODUCTION 


; pe LIFTING PRESSURE DISTRIBUTION and the damp- 
ing in roll are calculated on the basis of a second- 
order theory similar to the one introduced by Buse- 
mann! and extended by Van Dyke.** The airfoils 
considered here are infinite swept wings with arbitrary 
symmetrical cross sections. It is the authors’ opinion 
that from the study of these relatively simple airfoils 
effects of thickness on other airfoils can be estimated. 

Recent work by Milton D. Van Dyke? * indicates 
that second-order solutions of the partial differential 
equation of steady supersonic flow can be obtained 
by the use of iterative methods. The partial differential 
equation considered here is not the equation of steady 
supersonic flow, but it is similar, and we shall assume 
that the second-order solution can be obtained by 
iterative methods. 

It will also be assumed that the characteristics are 
the same for the first- and second-order solutions (this 
assumption was made in references 2 and 3). For 
steady plane flow the second-order solution’ found by 
an iterative method based on the preceding assumption 
yields the correct second-order pressure of the Busemann 
second-order theory. 

Unfortunately, no such justification of the preceding 
assumption is known to the authors for three-dimen- 


sional flows. 


THE PARTIAL DIFFERENTIAL EQUATION 


The partial differential equation to be used is a spe- 
cial case of the three-dimensional time-dependent equa- 


Received March 17, 1953. 


* Aerodynamicist, Airflow Branch, Exterior Ballistic Labo 
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+ Mathematician, Airflow Branch, Exterior Ballistic Labo- 


ratory 


tion for the potential function of a nonviscous com 
pressible fluid. This equation can be written as [refer- 


ence 4, Eq. (8) ] 


where c is the velocity of sound in the undisturbed 
stream ahead of the airfoil, y is the adiabatic exponent, 
and V is the velocity of the free stream. 

A perturbation will be used which is defined by 


2 = V(x’ + &) (2) 


Note that the perturbation potential, &, has been nor- 
malized through division by the free-stream velocity. 
Eq. (2) indicates that the coordinate axes have been 
chosen so that the free-stream velocity is parallel to 
the x’-axis. 

The coordinate axes of Eq. (1) are fixed in space. 
It is more convenient, however, to choose a new set of 
axes fixed to and rolling with the airfoil, the x-axis 
coinciding with the axis of roll. The relations between 
the components of the two sets of axes are 


, 


x= 

y = y’ cos pt’ + 2’ sin pt’ 9 
z = — y’sin pt’ + 2’ cos pt’ = 
t = 7?’ 


where the primes denote the fixed axes and p is the rate 
of roll of the airfoil. 

The flow over bodies rolling with constant velocities 
is independent of time; thus derivatives with respect 
to time in the new coordinate system are zero for the 
cases we shall study. It follows that the differential 
equation (correct to the second order) for the flow over 
bodies with constant p is given by 


2pM? 


7p (VPs — BB ey) + 


—B®,, “+ Py, + ®., =e 


M* [(y — 1) ®x(Pr2 + Py + Biz) + 20,072 + 
2,9®,, + 2®,9,, | (4) 


where M is the Mach Number and 6? = M? — 1. 
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Similarly, the pressure coefficient is found to be 


p 


it. 
| 


Cy = -2(¢, a ~ ?, | .) a ®,” =. ®,” = B°®,” (5) 


METHOD OF ITERATION 


We shall assume that Eq. (4) can be solved by an 
iteration procedure. The solution to the first-order 
partial differential equation is taken as the first approx- 
imation to the solution of Eq. (4). The first-order 
equation, obtained by neglecting the second order 


terms in Eq. (4), is 
—B*drr -— Pry -|- pz: = 0 (6) 
We assume that the second approximation can be 
found by substituting the first-order solution into the 
right side of Eq. (4) and solving the resulting non- 
homogeneous equation, which from Eqs. (4) and (6) is 
—B Yr + Voy + Ve at —(2pM* V) (Vorz — 2hry) i 
24? ((N a 1 )B*broy, = a Py Pry =F ¢:9::| (7) 
where 
N = (7 + 1)M?*/(26?) 
The solution of Eq. (7) will be referred to as the second- 
order solution. 


EXAMINATION OF POTENTIAL FORM 


It is helpful to investigate the type of solution ob- 
tained from Eq. (7). The first-order solution will be 


of the form 
d = agi(x, v, 2) + pgo(x, y, 2) + Egs(x, y, 2) 


where a is the angle of attack of the airfoil, € is a thick- 
ness parameter, and a, p, and ¢ are small compared to 
unity. It follows from the preceding expression and 
Eq. (7) that the second-order solution will be of the 


form 
Y = avh(x, y, 2) + prhe(x, y, 2) + 
e7h3(x, y, 3) + apha(x, vy, 3) + aehs(x, y, 3) + 


pehs(x, y, z) (8) 


The terms a*/,(x, y, 2), @hs(x, y, 3), and aehs(x, y, & 
are independent of the rate of roll and do not con- 
tribute to the lifting pressure due to rolling motion of 
the airfoil.* Since the present paper is concerned only 
with the lifting pressure due to rolling, all terms of the 
form a®/;, e°h3, and aes can be neglected in the follow- 
ing analysis. 

The following argument shows that p*/.(x, y, z) and 
aphs(x, y, z) do not contribute to the lifting pressure: 
The thickness parameter, ¢, is not present; therefore 
the airfoil can be considered a flat plate insofar as these 


terms are concerned. If the potential on the upper 


* These terms are associated with steady supersonic flow. 
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surface of a flat plate is expressed as 
P(a, Pp) = ag, + pg, + a®hy + p*h, + aph, 
the potential on the lower surface is 
P(—a, —p) = —agi — pe, + a®hy + p*h, + aph, 
The potential difference is 
Ab = 2ag, + 2pgr 


Since the lifting pressure can be found directly from 
A® it follows that p*/, and aphy do not contribute to 
the lifting pressure. 

The remaining term in Eq. (8) will be found by 
evaluating the right-hand side of Eq. (7), neglecting 
all expressions multiplied by a’, ae, p?, e®, and ap. 


BOUNDARY CONDITIONS 


Physical considerations require that the flow be 
tangent to the surface and that all velocity perturba- 
tions vanish upstream of the airfoil. These boundary 
conditions are expressible mathematically as (x, y, 2) 
= 0 upstream of the airfoil and g-‘Vs = 0, where s(x, y, 
z) = Ois the equation of the surface of the body and q 
is the velocity vector. 

The equation of the surface of the body may also be 


expressed as 
z= éf(x, y) 
where ¢€ is small compared with unity and f(x, y); thus, 
Vs = —1te(Of/Ox) — je(Of/Oy) + k 
Since the velocity, g, can be written as 
q = UV + Vo, + Vor) + (Vb, + Vb, + pz) + 
R( io. + Vw. + py 

it follows that the boundary condition on the body sur 
face is 
—(V + Vo, + V¥z)o, — (Vo, + V¥, + p2)o, + 

Vo. + Vy. + py = 0 
where o, = ¢«(Of/Ox) and o, = e(Of/Oy). 


The coordinate axes will be chosen so that the airfoil 


lies approximately in the s = 0 plane. The boundary 
conditions for the first-order solution are $(x, y, 2) 0 


upstream of the airfoil and 


o.\s=0 = (—py/V) + ao; (9) 


Similarly, the boundary conditions for the second-order 
solution are Y = 0 upstream of the airfoil and 


Velenq = OrOr\gu0 + Fy%y\s=0 — f(x, V)zz2\z=0 (10 


For the airfoils considered here, the first-order ve 
locity components are discontinuous across the Mach 
sheet from the leading edge. To evaluate the effect of 
the discontinuities on the second-order solution, we 
shall assume that the leading-edge Mach sheet has a 
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SWEPTBACK ROLLING AIRFOILS 


snall but finite thickness within which the first-order 
velocity components will be made continuous. Thus 
the discontinuities through the Mach sheet are replaced 
by continuous functions. This process is illustrated in 
Fig. | for ¢, for an airfoil with leading edge perpendicu- 
lar to the free stream. The effect of the discontinuities 
on the second-order solution will be found by obtaining 
the second-order solution within the Mach sheet and 
then letting the thickness approach zero.* 

The procedure just outlined will be illustrated by 
considering the case of an infinite rolling rectangular 
airfoil having a constant initial slope, ¢, along its leading 
edge. Only the Mach sheet above the airfoil will be 
considered, since the Mach sheets above and below 
are of the same form. The discontinuities through the 
leading-edge Mach sheet depend only on the initial slope 
of the airfoil; therefore the discontinuities in the first- 
order velocity components through the Mach sheet 
above the airfoil will be the same as those for a flat roll- 
ing airfoil at an angle of attack of minus e. 

The first-order potential function associated with the 
flow over the upper surface of an infinite rolling rec- 
tangular wing with a constant slope e is 


@ = [(py 1") — e] (x — Bz)/8 (11) 


where the airfoil lies approximately in the z = 0 plane.t 
The discontinuities in the first-order velocity com- 
ponents through the leading-edge Mach sheet above 


the airfoil are, from Eq. (11), 


Ad, = [(py/V) — e€]/8@ 
Ad, = 0 
Ad. = |[(—py V) + e] 


It will be assumed that the Mach sheet has a small 
thickness (see Fig. Ib). The velocity components 


within the sheet will be defined as 


or = [(py/V) — e] (x — Bs) (/8) ) 
gd, = 0 (12) 
o. = —([(py V) — e] (x — Bz) iN 


where / is the thickness of the Mach sheet in the x direc- 
tion. Note that within the Mach sheet 


— 26... + by + Ge = 0 


and that the velocity components are continuous func- 
tions in the neighborhood of it. 


* An alternate approach would be to attach a small surface 


to the leading edge of the airfoil which would cause the first 


order velocity components to be continuous. The effect of the 


discontinuities would then be evaluated by a limiting process in 
which the width of the small attached surface would approach 
zero 

t Eq. (11) can be found by the integration of a source distribu 
tion (see references 5 and 6); however, for simple problems such 
as this, a heuristic method generally yields the result with a 


miuimum of effort. 
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a Leading Edge Moch Sheet 


With Zero Thickness b Leading Edge Moch Sheet 


With Assumed Thickness 





pz ” xe pt, ' 








d Plot Of $,Across Leading 
Edge Mach Sheet With 
Assumed Thickness Along 
The Line ZZ, 


c Discontinuity In $, Across 
Leading Edge Moch Sheet 
Along The Line z=z, 

Fic. 1. An illustration of removing the discontinuity in ¢, 
across the leading edge Mach sheet by assuming the Mach sheet 
to have thickness. 


From Eqs. (7) and (12) the differential equation 


within the Mach sheet is 


o 2M*pey 2N(x — 62) 
— Bre + Vy + 22 = = IV l + ] (13) 


It is well known * 7 that the solution of the nonhomo- 


geneous equation 
—Bbr2 + Vw + v2: = F(x, y, 2) 


is given by 


1 . . . F(é, . ¢) ; 
caine De J J J R dé dn dé + 

l “SP 1 

a | / (ve ry ov) da_ (14) 


R = V (x — &)* 


Nn= —1B?», + jro + kos 


where 


— By — n)? — Bz — £)? 


(&, 
are the direction cosines of the outward normal to the 
element of area da, s is an arbitrary surface that cuts 
the forward Mach cone from the point (x, y, 2), v is the 
volume inclosed in the forward Mach cone from the 

° ° ; bd 
point (x, 4, z) and the surface s, and the symbol S'S 


denotes the finite part of an integral as defined by Hada- 


n, ¢) are auxiliary cartesian coordinates, (1, v2, ¥3) 


mard in reference 7. 
For planar problems the surface integral in Eq. (14) 


can be reduced to (reference 6, p. 22) 








702 JOURNAL OF THE AERONAU 


TICAL 











x 


Fic. 2. Coordinates for sweptback wing. 
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It follows that the potential function within the 


Mach sheet is 


V(x, y, 2 


y= Ref ff: af a 4 EA) 5 
l 
1 f[ f1 
neat f [0-)uan (15) 


where 5; is the area of the z = 0 plane which is enclosed 
within the Mach sheet and the forward Mach cone from 
the point (x, y, 2). 

To evaluate the second-order solution on the down- 
stream side of the Mach sheet, we take (x, y, z) to be 
an arbitrary point on the surface x = Bz + /. Since 
the thickness of the Mach sheet is small, it follows that 
for given values of — and ¢ the variation of the inte- 
grands of the vdbane and surface integrals is made up 
almost completely of the variation of 1/R. This al- 
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lows the substitution of y = 7 in the integrands except 
in R, and this in turn allows integration in the n direc. 
tion. Then the potential at a point (6z + /, y, z) can 
be expressed as 


Mpey 
¥(62 + /, y, 2) = 2 fi fe x 
2N py yf 
I t — Be) |dé a E dt 
r + ] ( at) |as +2 € i E dé 


Performing the remaining integrations and taking the 
limit as / approaches zero yields 


Ay|,-2 = M2(N + 1)peyx/ Vp? (16 


The value of the discontinuity in the second-order 
potential is given by Eq. (16) for that part of it which 
contributes to the lifting pressure due to the flow over a 
rolling airfoil with a nonsweptback leading edge of 
constant initial slope. Similarly, when the leading 
edge is swept back, the second-order discontinuity is 
found to be 


, M?°mpe/B?m*N 
AY} mx y-Az=0 = 


1}(2Ay — 2)2 (17) 
2Vr? Nes ° ye F 


where m is the cotangent of the sweepback angle A (see 


Fig. 2) andA = V ‘B? m? — 1. 


SOLUTION OF DIFFERENTIAL EQUATION 


Having determined the differential equation and its 
boundary conditions, we proceed to its solution. 

The coordinate axes are chosen as indicated in Fig. 
2. The equation of the upper surface is represented 
by 


z= e[(mx — y)/m] (18) 
The first-order solution is 


p 
¢ = (mx — y — \s) [—mx + 


(282m? — 1)y—Az] — 
OV a! 8 ) ] 
m (= —-y- 2) 
f : (19) 
A m 


The second-order potential function must satisfy the nonhomogeneous equation (neglecting terms multiplied by 


p* and e?) 


— Bez ar Vyy + V2 


chosen to be 


— pel (Nem 
PERM w 


By inspection the complementary function, 4, 
(17)] and 


+ 1) yf + =f" + 


= ee (“= 
" &¥§ \! 


where the primes denote derivatives with respect to [(mx — y — Xz) /ml. 


“ m 
l ) a 
2 


(20) 


B'm*N .,  NGm* _ 
aie | xf 
nn ~~ 


The particular integral of Eq. (20) was 


apt = (en ‘ 1) f VB?m*4 | >| 
2*f’ — sf — : xz (21) 
: r3 2 . > . 


of Eq. (20) satisfying the boundary conditions [Eqs. (10) and 
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— B?(Ox)rr + (02) yy + (02)22 = O 


is found to be 


pe ae ) , Pee ”( M?*NB‘m‘. _.. Bim + ‘) 
, = + 1)\(mx — y)* + A*3s")/'(0) -{ — 387 + Il - x 
V \ 4\4 A | : | * 3 ? m* 


| m(2 4 ) m? {~ ieee fu) d | my shen (Mere :) (= "| 4 
—_ 2 = u) du i 2 = 5 — 
A \A " A? J0 3 ? \; nN 
Bm? [7 - le m ( M?NB?2m4 , coals Alana \ 
f(u) du | — - x2 +m?+ 1 f(u) du (22) 


\” 
The part of the potential function contributing to the lifting pressure is 
Y= A+ % 


From Eqs. (5), (19), and (23) the lifting pressure for symmetrical airfoils is 


: 4pm — tbe M?m?| 1 /NB?m? ; NBim4 be 
AC, = (ms — p-m-y) — _ — + 1) (mx — y)f’(O) + x2 — liv — 


VAr3 VM At 
NB*m? . 2NB?m? . : 
—— <= | imax — = + 117 (24) 
ea ‘ - 


AC, = —(4py/VB) [1 + (€/8)MPNf'(x) | 


For unswept wings, #77 = © and Eq. (24) reduces to 
| 


Eq. (25) indicates that for the unswept wing the second-order effect is a function only of the local slope of the 


airfoil surface. 
Eqs. (24) and (25) are linear in e. 
to different types of thickness distributions for the same plan form. 


This leads one to inquire whether the effect of thickness is linear with regard 
Suppose the equation of the surface 1s ex- 


pressed as 
= efi(x, vy) + eofo(x, y) 


The first-order solution is of the form 


bd = ag,(x, ¥, 3) + pgolx, V, 2) + e1g3(xX, V, 2) €224(X, V, 3) 
The second-order solution will be of the form 


Y = ath,(x, y, 3) + prhe(x, ¥, 3) + arhs(x, vy, 2) + @*hi(x, y, 2) + 

aphs(x, y, 2) + aehs(x, y, 5) + aehs(x, y, 3) + pehe(x, y, 2) + peho(x, y, 2) + €;€2 (x, VY, 2) 
By the same arguments as were used in the section on the form of the second-order potential function, it can be 
shown that only the terms peg and pests contribute to the lifting pressure due to rolling. It follows that the effects 
of thickness on the lifting pressure can be added linearly up to the second order. The linearity of thickness effects 
on the lifting pressure is also true for angle of attack, steady pitching, and other types of motions. 





THE DAMPING IN ROLL 


The damping in roll coefficient, C;,, can be written 


as 


. . 
C,, = lim (pb/2V = | | yAC, dx dy (26) 
p—0 I0/Z )O0. ./ wing area 


where / is the wing span and S is the wing area. Sub- 
stituting Eq. (24) into Eq. (26) vields 


, 2m el/°NB?m? | ws 
Cy = — Pag - {(d) (27) 
3A ed 


Eq. (27) indicates that the effect of thickness on the 


C,, of an infinite swept wing with an arbitrary sym 
metrical thickness distribution varies directly as the 
thickness of the trailing edge. 

Fig. 3 presents the variation of the C,, with Mach 
Number for a 10 per cent infinite wedge for various 
sweepback angles. This figure indicates that the effect 
of thickness on the C,, increases slightly with the 


sweepback angle. 


CONCLUDING REMARKS 


It is possible to estimate fairly accurately the effect 
of thickness on the C,, of many airfoils not considered 
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Fic. 3. The variation of the damping in roll coefficient, Cip, reference 8. This figure indicates that the agreement 


with Mach number for 10 per cent thick infinite wedges for vari 
ous sweepback angles. 


in the analysis. As an illustration consider the C), of a 
rectangular airfoil with a wedge cross section such as 
considered in reference 8. An estimate of the effect 
of thickness can be made by increasing the C,, of the 
linearized theory by the same percentage as the thick- 


between theory and experiment has been improved 
considerably. 

Since the theory in this paper is closely associated 
with the Busemann second-order theory, it seems 
likely that the present results have limitations similar 
to those of the Busemann theory (which are treated in 


reference 9). 
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A Method of Estimating the Compressive 
Strength of Optimum Sheet-Stiffener Panels 
for Arbitrary Material Properties, Skin 
Thickness, and Stiffener Shapes 


W. R. MICKS* 
The RAND Corporation 


SUMMARY 


Developed in this paper is a simple analytical method of syn 
thesizing allowable stress curves for optimum sheet-stifiener 
This method may be used with any struc- 


compression panels. 
The only material 


tural material and any stiffener shape 
properties required in the analysis may be determined from a 
compression stress-strain diagram for the material being 
used 

Also included in the analysis are the effects of varying the rela- 
tive amount of panel material allocated to the skin. 

Analytical results are compared with the results of tests con 


ducted by the NACA. 


SYMBOLS 


A, = area of one stringer 

A = area of skin between stringer centerlines 

b = stiffener spacing 

b, = distance between points of attachment, measured 


across the stiffener cross section (see Fig. 9) 
b, = distance between points of attachment, measured be- 
tween stiffener cross sections (see Fig. 9) 
= radius of gyration of the stiffener cross section 


p: 

P, = radius of gyration of the sheet-stiffener combination 

t, = skin thickness 

i = average panel thickness (average cross-sectional area 
per unit width) 

(t/b) = thickness-width ratio for a flat plate element, in gen 
eral 

- = distance from the centroid of one stiffener to the skin 
centerline 

k, = t,/t 

k, = x,/p- shape factor measuring stiffener depth 

Re = /,/A. = p,/A, shape factor measuring bending ef 
ficiency of the stiffener cross section 

ky = shape factor measuring stiffener width 

k, = (p.2/A.)(to/bo) shape factor measuring bending ef 
ficiency of the basic stiffener shape 

k = p,/i shape factor measuring bending efficiency of the 


sheet-stiffener combination 


INTRODUCTION 


MM WORK HAS BEEN DONE on design and strength 
testing of sheet-stiffener compression panels. 
Several comprehensive testing programs have been car- 
ried out to determine optimum panel configurations 
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and to develop design charts for panels made of various 
types of material and having various shapes of stiffen- 
ers.' A separate series of tests was conducted for each 
different material and/or stiffener shape. 

However, since new structural materials are de- 
veloped from time to time, there is need for a method 
of estimating the strength properties of sheet-stiffener 
panels made of a new material, without resorting to a 
complete new testing program. The recent introduc- 
tion of elevated temperatures into the design conditions 
for aircraft structures has greatly increased this need, 
since one of the primary effects of elevated temperature 
on the design of structural elements is to change the 
physical properties of the material. 

Strength properties of compression panels are best 
presented and compared on a plot of allowable com- 
pressive stress, f, versus the structural index, q/Lo, 
for wide columns," * where q is the load per inch width 
and Ly is the effective pin-ended column length. Shown 
in Fig. | isa plot of f vs. g/ Lo for aluminum-alloy panels 
having various stiffener shapes. (This figure is re- 
produced from reference 2.) These curves represent 
the highest stresses attained in tests using the material 
and stiffener shape indicated on the figure. Such curves 
have been called ‘“‘optimum”’ curves. 

However, one important thing that such plots of the 
“optimum” stress have not ind:cated is how much of 
the material is in the stiffeners and how much is in the 
skin. This information is necessary in making a pre 
liminary estimate of the torsional stiffness of a wing 
using sheet-stiffener panels. Some of the configura 
tions corresponding to the curves of Fig. | tend to have 
a relatively small amount of material in the skin and 
therefore are not likely to be used in the design of a 
wing panel for a wing required to have high torsional 
stiffness. The basic parameter that indicates the rela- 
tive amount of material allocated to the skin is /,//, 
the ratio of skin thickness to average panel thick- 
ness. 

In the actual design of a high-speed wing, a minimum 
allowable skin thickness is usually specified on the basis 
of aeroelastic effects. Therefore, any optimization of 
the sheet-stringer panel design should contain /,// as an 
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al 
estimating the effects of different stiffener shapes on 
— ls wreioas: stall A ston ee 
{_ the strength properties of sheet-stiffener panels, without D 
1 . i 7 | conducting a complete testing program for each different 
“Z ae im =e shape. 
a ° . at 
! | The present paper develops a simple method for esti- I 
— mating the strength properties of optimum-designed e 
Centroidal axis of sheet-stiffener panels (1) made of any material, (2) fe 
- Stringer centroid a containing any specified fraction of the panel material 
in the skin, and (3) having any shape of stiffeners. The v 
Ej aia as ef estimated strength properties from this method are s 
0-2 ne ere Coe compared with those obtained from tests on a plot of 
: ’ — ; _ the f vs. q/Lo. 
independent variable. For this reason, a method is The only material properties required for the analysis I 
needed for determining the maximum allowable stress 2. he determined from a compression stress-strain 
for any specified value of t/t. _ , diagram for the material being used. 
The third variable that should be accounted for is 
the stiffener shape. It is desirable to have a method for 
Ta . : _ CROSS-SECTIONAL SHAPE FACTOR FOR WIDE COLUMNS 
*In an analysis of compression panels having stiffeners of 
uniform thickness, Zahorski* used the ratio of skin thickness to For wide columns the tangent-modulus column- 
stiffener thickness as a parameter in expressing the relationship endian eameat nial Sow einen ta Cs Sevens ; 
between the relative amount of skin and stiffener material. In ny ee ay ee ree niles 
the case of formed stiffeners this parameter has a definite mean- f= g R tp,t/3 q/Lo */, (1) 
ing. However, for extruded stiffeners composed of elements , 








having different thicknesses, this particular ratio is no longer where 
applicable. It is then desirable to use the parameter ¢,/?7, which 
represents a slightly higher level of abstraction. qd = load per inch width 
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Ly = effective pin-ended column length The total moment of inertia Jy of one bay of the panel 
¢, = tangent modulus about its centroid is 
ks = pr i =\9 as = 
: : 4 ; Ir = I.,, + I, + A(x, — ®)? + Ask (7) 
pr = radius of gyration of the wide-column cross 
section where 
7 = average thickness _ , ' ; , 
/,, = moment of inertia of the stringer about its 
The overall shape factor k; gives a measure of the centroidal axis 
bending efficiency of the wide-column cross section. /,. = moment of inertia of the skin element about 
This factor will be affected by stiffener shape and spac- its centroidal axis 
ing, skin thickness, and other factors that determine the , ' : ae — ? a 
Substitute Eqs. (3), (4), and (6) into Eq. (7) and 
overall panel cross-sectional shape. on ; 
write the moment of inertia as 


For a fixed value of structural index g/L, a measure 
of loading intensity,” an increase in the value of k3 will Jr = J., + 7;, + t0(1 — Rix Pk, + tybx2(1 — ky)* (8) 
increase the long-column buckling stress [see Eq. (1) ]. sa - ’ : 

Esai & : Ml | The moment of inertia of the sheet element about its 
However, as k3 is increased, perhaps by making the 
stiffeners thinner and deeper, the local buckling stress 
is decreased. The theory of optimum design assumes 
that the highest allowable stress for the column will _ . eer : : 
The moment of inertia of the stringer about its cen- 
occur when the long-column stress equals the local- ; , . , 
: : troidal axis may be written in the form 
buckling stress. It is therefore necessary to develop : 
an expression for ks in terms of design parameters that I,, = k,A2” (10) 


centroidal axis is 


I,, = bt,8/12 (9) 


will be known. oe ome eae ale 
; ; E ; : where k, = J,,/A,°, which is the shape factor indicating 
As will be shown later, k3 is a function of stress, ma- Nye ; : . 
. j a aaa the efficiency of the cross-sectional shape in bending. 
terial properties, and panel geometry. Substituting ..,. aie ‘ : 
ag a ae  < This factor k, is identical with the shape factor k; de- 
this function for k3 in Eq. (1) will be equivalent to : : : 
: , veloped in reference 2. 
equating the long-column buckling stress to the local : nS SS ae . ‘ 
rtng ' . Substituting Eq. (3) for A, into Eq. (10), 
buckling stress, while also specifying a given set of 
boundary conditions that depend on material proper- T,, = Rts*b?(1 — k,)?/k,? (11) 
ties, stiffener shape, and relative skin thickness. . . ¥ . ° . 
; : , ; é' . Eqs. (9) and (11) are substituted into Eq. (8) to give 
After the expression for k; is determined, Eq. (1) 
will be used to plot curves of f vs. g/Lo for the optimum 


k,t,*b?(1 — k,)? | bts 
compression panel that represents the specified bound- = 


3 
. - 2 pa 9 
I, be 4 12 + bt,x"(1 k,) (12) 
ary conditions. 
The radius of gyration for the panel is 


DERIVATION OF THE EXPRESSION FOR SHAPE FACTOR k; Ip] 
Pr = 1 (13) 
The shape factor k; has been defined as ks = pr/t. ii 
The following derivation gives the development of an Since Ar = t,b/k, and 7 = t,/k,, the expression for 
equation that expresses p,/? in terms of stiffener shape (= p_/7) is written as 
factors, the skin-thickness parameter ¢,/7, and the stress. j T,t/tp.4/ 4/2 14 
‘ . ~~ . . . . ., = mp “P - ( ) 
Shown in Fig. 2 is an idealized sheet-stiffener panel ro i 
with dimensional notation used. The area of one Eq. (12) is now substituted into Eq. (14) to give the 
stringer can be expressed expression for ks, 
A, = (i — t,)b (2) kk,b(1 — k,)* k, — x.?k (1 — ky) )” . 
kh» = $+ + = ( 15) 
i. 12 [” 


Let k, = t,/Zi and substitute for 7, giving stringer area as 
In order to express x, in nondimensional form, it is 


A, = t,b(1 — k,)/k; (3) : 
= tere necessary to introduce a shape factor for the stiffener. 
rhe area of skin between stringers is Let 

A. = t.b (4) k, = x,/p, (16) 
rhe distance of the panel centroid from the centerline y yore p: = radius of gyration of the stiffener The 
of the skin is given by radius of gyration of a stiffener is given by 

# = A.x,/(A; + A;) (5) p: = (I.,/A,)” (17) 

Substituting for A, and A,, Substituting Eqs. (11) and (3) for /,, and A., respec- 
& = x,(1 — k,) (6) tively, into Eq. (17) 
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Fig. 3— Examples of cross-sections having equal 
values of the shape factors &, and A, 


es [k,t.b(1 — k,) k,]’ . (18) 


The expression for x, is now obtained from Eq. (16) 
by substituting Eq. (18) for p, and solving for x.. 


x, = k,[kt.b(1 — k,)/k]” (19) 


This expression for x, is now put into Eq. (15), and 
simplified, giving 


kyk(l — k,)2°1 + Rk) Re] 
kz = | obei\ = + | (20) 
t,/b 12 


where k, = p,”/A, or I,,/A,’. 

The next step here is to determine the relationship 
between k; and the local buckling stress. 

Since the buckling stress of a flat-plate element is a 
function of the thickness-width ratio (¢/6) for that 
element, it is possible to relate k; to buckling stress 
through the factor ¢,/b of Eq. (20). However, all other 
factors in Eq. (20) should first be made independent of 
thickness-width ratio. 

It can be shown that the shape factor k, is practically 
independent of the ratio ¢/6 of the stiffener elements 
if the following conditions are fulfilled while the size 
or thickness is being changed: 

(1) The skin thickness is small relative to the depth 
of the stiffener. 

(2) The stiffener cross section is composed of rela- 
tively thin-walled elements—i.e., the wall thicknesses 
are small relative to the overall dimensions of the stiff- 
ener. 

(3) The ratio between values of thickness-width 
ratio (¢/b) for all of the stiffener elements is held con- 
stant while the actual value of (¢/b) is being changed. 
(This is also a necessary condition for having all stiff- 
ener elements buckle at the same stress. ) 

These three conditions are usually fulfilled in aircraft 
sheet-stiffener design over most of the range of variables. 
In cases where this is not true a correction factor may 
be derived for any particular type of configuration. 

The nondimensional shape factor, k, = (p.°/A,;), is 
not independent of the (¢/b) ratio of the stiffener ele- 
ments. Therefore, a new shape factor was developed 
for thin-walled sections which includes the bending- 
efficiency factor p.°/A, but does not vary with wall 


thickness. This factor, denoted by &,, is defined as 


k, = (p-” A,)(to bo) (21) 
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or 


k, = k, (to ‘bo) (29 


where f)/bo is the thickness-width ratio for a specified 
element of the stiffener cross section (the choice of this 
element is discussed later). 

This factor is independent of the ratio ¢/b of the stiff 
ener elements if the three conditions outlined above are 
fulfilled. 

The shape factors k, and k, were developed primarily 
for use with thin-walled cross sections composed of 
straight-line elements. However, it will be shown later 
that they can be applied to optimum cross sections con- 
taining a combination of straight and curved elements, 

Fig. 3 shows three variations of the same basic cross- 
sectional shape which have the same values of k, and 
ks. 

One of the basic assumptions made in this analysis 
is that, for the condition of overall optimum design, 
all elements of the panel must buckle at the same 
stress.* Since no element is made overstrength, the 
buckling stress will equal the ultimate stress. This 
assumption will apply to the stiffener design, as well as 
the sheet-stiffener combination. 

Solving Eq. (22) for k, and substituting into Eq. (20), 


kki(l — k,)2(1 + R,?k,) key" 
k; = 
| (t,/b) (to/ bo) + | \ 


23) 
In evaluating k,, it is desirable to use the ratio of 
t/b determined from a stringer element having the 
same edge-restraint as that of the sheet between stiff- 
eners.t For this condition and the assumption of 
optimum design conditions mentioned above, 


to bo = f, b (24 


If the stringer has no element with the same edge- 
restraint coefficient as the sheet elements, the use of 
some other stringer element gives the relation 


to/bo = k-(t,/b) (25) 
where 


edge-restraint coeff.cient of the sheet elements 


edge-restraint coefficient of the stiffener ele- 
ment used 


The factor k; can be expressed in terms of (¢/b)° by 
substituting Eq. (25) into Eq. (23) to give 


kk(l — R21 + 2k) ke]? 
ks = t t , t 4 t (26) 
k,(t,/b)? 12 
Since the buckling stress of a flat-plate element is a 
function of (¢/b)*, the above expression for kj may be 
related to stress through the plate-buckling equation 


* The maintenance of buckle-free conditions up to ultimate 
load was found to be in agreement with optimum design condi 
tions determined by the testing program outlined in reference 4 

+ A discussion of edge-restraint coefficients existing at the op 


timum design condition is presented in Appendix (A). 
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Fig. 4-Buckling stress vs 2 for a simply-supported flat plate 
of aluminum alloy 


fer = [kw?En/12(1 — pw?) ](t/b)? (27) 
where 7 is the Stowell effective modulus factor (see 
reference 5). 

The factor k; could be expressed as a function of stress 
by solving Eq. (27) for (¢/b)* and substituting in Eq. 
(26) for (t,/b)?. However, it is more convenient here 
to use the two equations |Eqs. (26) and (27)] sepa- 
rately, as illustrated in the next section. 

For the condition of optimum design, the local buck- 
ling stress of the panel elements must equal the long- 
column buckling stress. Therefore /,, of Eq. (27) must 
equal f of Eq. (1), and k3 can be determined as a func- 





tion of f for use in Eq. (1). 

Since comparisons of various panel designs will be 
made on the basis of f vs. g/ Lo, as shown in Fig. 1, a 
procedure will be outlined for plotting these curves 
from the relation given by Eq. (1). 

The example curves used to illustrate this procedure 
are for 24S-T aluminum-alloy panels having Z stiffen- 
ers. The shape parameters are based on the Z-section 
used in references 1 and 4 (flange length = 0.4 web 
length), and have the values given in Table | of this 


report. 


PROCEDURE USED IN OBTAINING FINAL CURVES 


For the material chosen, a curve is plotted from Eq. 
Such a curve is shown in Fig. 4 


(27) giving f vs. b/t. 
(Constants used in 


for two types of aluminum alloy. 
plotting Fig. 4 are k = 4, u = 0.3, and E = 10.5 X 
For 24S-T material the values for 


10° Ibs. per sq.in. 
The values of 7 for 


n were taken from reference 5. 
the 75S-T material were computed from stress-strain 
test data of reference 6.) 

Values of the shape factors k, and k, are determined 
for the basic stiffener shape being investigated. Then 
k; can be determined as a function of stress; k,, by use 
of Eqs. (26) and (27). It is assumed here that the 
sheet between stiffeners is simply supported at the 
edges. Therefore, a value of k; = 1.00 was used. 

Fig. 5 shows how the local buckling stress varies with 
k; for several values of skin-thickness parameter ,. 

Fig. 6 shows a plot of k; vs. k, for a constant value of 
stress equal to 40,000 Ibs. per sq.in. When k, equals 
unity, all of the panel material is in the skin. As 
k, is decreased, while holding f constant, material 
must be taken from the skin and put into the stiff- 


eners. 
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Fig.5—Stress vs A, for 24S-T Z-stiffener 
compression panels of optimum proportions 


It is interesting to note that, as the value of k, is de 
creased from unity, ks; reaches a maximum value and 
then drops off again as still more material is transferred 
from the skin to the stiffeners. 

As the value of k, is decreased toward zero, while 
holding the stress constant, the stiffener spacing will 
tend to approach zero. Therefore the actual width of 
the stiffener would place a lower limit on the value of 
k,. However, it would be impractical to design a panel 
for a value of k, below that which gives the maximum 
value of 3. 

Although the end point of the curve at k, 1.0 is 
established, a portion of the curve for relatively large 
values of k, (say, k; > 0.8) may tend to be somewhat 
too high, since there is doubt that the resulting stiff- 
eners in this range will have sufficient stiffness to stabi- 
lize the sheet and produce the simply supported edge 
conditions. This would, in turn, make the 
f vs. g/L» too high for these same values of ;,. 


curves of 
This 
effect could be accounted for by a more elaborate ver- 
sion of the analysis made here. However, a few check 
tests would still have to be made in order to verify the 
results. It is unlikely, however, that any designs would 
be made at either extreme of the range of k,. 

The v 
vary somewhat with the stiffener shape used (see ex- 
ample curves Y-stiffeners). For the material 
being used, the relation between tangent modulus and 
7, which is 


ralue of k, giving maximum allowable stress will 
for 


stress may be shown by a curve such as Fig. 
) 


reproduced from reference 3. 


Solving Eq. (1) for g/Ly and using E, = rE gives the 


expression 


1 


q/Ly = f'"/eE rks (28) 









































Fig.6—Variation of 4, with k, for 24S-T 
Z- stiffener panels (£=40,000 psi) 


For any value of f chosen, we can determine corre 


sponding values of k; and 7 from Figs. 5 and 7, re 
Substitution of these values into Eq. (28 
Then plotting 


spectively. 
gives a corresponding value of g/Lo. 
f vs. g/Lo results in curves as shown in Fig. 8. The 
highest curve shown is that corresponding to the maxi 
mum value of ks (see Fig. 6). 


COMPARISON OF RESULTS 


The dashed line of Fig. 8 represents the “‘optimum’ 


tests as described in reference | 


curve obtained by 
It can be seen that this curve lies fairly close to the 
highest curve obtained by the present method.  Pos- 
sible sources of difference are: 

(a) Difference in material properties: The curves 
obtained by the present method are based on ‘“‘nominal’ 














aluminum alloy taken from Figs 
24S-T 


properties of 24S-T 
4 and 7. If the tangent-modulus curve of the 
material used in the tests differs from the values used 
here, particularly in the yield region, the f vs. q L 
curves will be affected accordingly. 

(b) Shape factors: In deriving the f vs. g/Lo curves 
by the present method, the shape factors were assumed 
This 
condition was probably not exactly fulfilled in _ tests. 





to be constant over the entire range of loading. 
) Differences in end-restraint condttions: iS 
iiaaty that the edge-restraint conditions for Pot 
of the test specimens were somewhat different from the 
condition assumed here, particularly if some of the 
panel elements buckled before the maximum load was 
reached. However, it has been assumed here that, for 
optimum design, all of the panel elements buckle simul- 
When is fulfilled, edge-re- 


taneously. this condition 
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Fig. 7 — Variation of tangent modulus with stress—typical values 
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Fig. 8 — Allowable stress vs structural index for 24S-T Z- stiffener 
compression panels having optimum proportions 
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Fig.9— Additional symbols used for 
two points of attachment 


straint coefficients seem to vary only slightly from the 
theoretical values for simple edge support [see Ap- 


pendix (A)]. 


STIFFENERS HAVING Two PoINTs oF ATTACHMENT 


The first portion of this paper dealt with the case 
having the stringer attached at only one point. How 
ever, sirice many stiffener shapes are attached at two 
points, the analysis is extended to cover this case. 

The sketch shown in Fig. 9 gives the nomenclature 
used for two points of attachment. 

It can be seen that if the spacing 0), is not equal to 
b., the two elements of sheet will not buckle at the same 
stress. Because of practical considerations in fabrica 
tion, we will deal only with the case where b, 2 /.. 

The stiffener spacing can be written 


(29) 


b=b,4+ b, 


Substituting Eq. (29) into Eq. (20) gives the expres 
sion for k3 as 
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kkA(l — ke2(1 + R,7R,) 
k; ” | ; . ua -+- 


‘,/0, 
k,kA1 — k,2)(1 + k,7hk,) 4 
- (30 
t,/b- 12 
Expressing /,/b, in terms of /¢,/b, gives the relation 
ship 
b, l b. tse 
= ( [R + (R? + R)”] (31 
i, 2 Xt, 
where 
R = k,ko?(1 — Ry) /R: 
= 0/0, 


The factor ky is a stiffener shape factor that describes 
the attachment. This 
factor is independent of the ratio ¢/b for the stiffener 


distance between points of 
elements, subject to the three conditions listed previ 
ously. 

Substituting Eq. (3!) into Eq. (30) and simplifying 
gives the expression 
k3 = 

kk — R,)?(1 + R7R,) (1/2) [2 + R + (R? + R) 

t,/b, 


Since it has been specified that the sheet between 
stiffeners must buckle at the same stress as the stiffener 
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Fig, |! —Allowable stress vs structural index for 75 S-T Y-stiffener 
compression panels having optimum proportions 


elements, Eqs. (22) and (25) are substituted into Eq. 


(32) to give 


>). + k,2k,)(1/2)[2 + R + (R? + R) 
Bee b.)” 


This expression is used in the same way as before in 
obtaining the final curves of f vs. g/Lo. 

Shown in Fig. 10 are the final curves obtained by 
using Eq. (33) and Fig. 4, together with the shape 
Table 1. The 
“optimum” test curve was taken from reference 1. 

The ratio of sheet width 5, under the stiffener (see 


factors for the Y-stiffener as listed in 


Fig. 9) to sheet width 6, between stiffeners can be ob- 
tained by solving Eq. (31) for b./d,. 


b./b, = (1/2)[R + (R? + R)”] (34) 
For any specified value of ¢,/7, the parameter ks will 


have its greatest value when b,/b, equals unity. For 


TABLE | 
Stiffener Shape Factors 
k k, ky 
Straight-Y 0.0935 1.33 1.48 
Curved-Y 0.109 1.34 1.48 
Z 0.089 1.36 


any value of b,/b, less than unity, the sheet element 
under the stiffener will not be stressed to its buckling 
The 
condition where b, is greater than }, has not been con- 


load, thus introducing a source of inefficiency. 


sidered here because of practical considerations in de- 
sign and fabrication. 

For the Y-stiffener panel, the highest curve of f vs. 
q/ Lo is obtained when ¢t,/? = 0.3. 

In the case of the 75S-T Y-stiffener panel, the sheet 
material and the extruded stiffeners have somewhat 
different material properties. No attempt has been 
made here to account for this difference in plotting 
the curves of Fig. 11 for 75S-T Y-stiffener panels. The 
relationship between f and }/t given in Fig. 4 for “‘nom1- 
nal’’ 75S-T material was used in the computations. 
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a. Straight web b. Curved web 


Fig. l12— Y-stiffener sections designed by the NACA 


EFFECT OF IMPROVING THE EFFICIENCY OF 
STIFFENER CROSS SECTION 


THE 


Following the test program conducted by the NACA 
on Y-stiffeners having straight webs, some work was 
done on a more efficient cross section having curved 
webs,’ which is shown in Fig. 12(b). 

In using the present method with a stiffener section 
having a combination of straight and curved elements, 
it is necessary only to use the thickness-width ratio of 
the straight element in evaluating k, and assume that, 
if the stiffener is efficiently designed, the curved por- 
tion will buckle at the same stress as the straight por- 
tion. 

Comparison of the values of k, in Table | show the 
curved web is somewhat more efficient than the straight 
webs. Therefore the curve of f vs. g/L» for the curved- 
web Y-stiffener will lie above that of the straight-web 
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Fig. '3- Comparison of allowable stress for 
75S-T sheet-stiffener compression panels 
having straight-web and curved-web Y- stiffeners 


SCIENCES OCTropspeEeR, 1953 
Y-stiffener. Fig. 13, comparing the highest allowabk 
stress curves for these two stiffener shapes, illustrates 
the effect of improving the stiffener efficiency by the 
amount shown in Table 1. The dashed lines represen; 
curves derived by the NACA from test results on thes 
two sections.’ It can be seen that the amount of jm 
provement obtained from the curved web is about the 
same for the analytical curves as for the curves based 
on tests, although the absolute values of stress do not 
agree exactly. 


CONCLUSIONS 


The present paper outlines a simple method for syn 
thesizing allowable stress curves for sheet-stiffener 
compression panels of optimum design. As shown in 
Figs. 8, 10, and 11, these curves closely approximate 
corresponding curves that have been obtained from 
tests. 

The basic method outlined here can be used to evalu- 
ate quickly, simply, and inexpensively the merits of 
any new structural material considered for use in sheet- 
stiffener wing panels. The only material properties 
required are those that can be derived from a compres- 
sion stress-strain diagram for the material. 

The synthesized curves are more general than the 
“optimum” curves based on tests, because they can be 
made to show the effect of varying the relative amount 
of skin material. This leads to a method of eva'uating 
structural weight penalties associated with increased 
skin thickness. 

It appears that with the aid of a few carefully de- 
signed tests the present method could replace much of 
the testing work ordinarily done in the investigation ol 
This 


is particularly true when design curves have already 


the strength properties of sheet-stiffener panels. 


been determined for a given material and _ stiffener 
shape and it is desired to evaluate the effect of changing 
the material or the shape of the stiffener. 

The stiffener shape factors developed for this analysis 
may be of value in the design of improved stiffener 
sections, since they furnish a simple criterion by which 
the efficiency of new shapes can be quickly and easily 
checked. 


(A)—EDGE-RESTRAINT COEFFICIENTS 
OPTIMUM DESIGN CONDITIONS 


APPENDIX FOR 


It may be desirable to include a brief discussion on 
how optimum design principles apply to the question 
of edge-restraint of the different flat-plate elements of 
a panel cross section and to the interaction between 
adjacent elements at the buckling load. This discus- 
sion deals with sheet-stiffener panels composed of flat- 
plate elements attached along one or both of the un- 
loaded edges. 

Reference 8 presents charts showing computed values 
of edge-restraint coefficient k of the plate-buckling equa- 
tion [see Eq. (27)] for the individual flat-plate elements 
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namely, the I-section, 


of several thin-walled shapes 
These 


channel, Z-section, and a rectangular section. 
charts illustrate how the effective value of k changes 
when the proportions of the cross section are changed. 
The relative values of ¢/6 of the elements are varied 
from the condition where the outstanding flanges buckle 
first (and are restrained by the web) to the condition 
where the web tends to buckle first and is restrained 
The parameters plotted are k vs. 


by the flanges. 
The subscript F 


by /bp for constant values of tp/tw. 
denotes an outstanding flange (supported along one 
edge) and W denotes the web (supported along two 
edges). 

These charts are of value in determining the condi- 
tions of edge-restraint which exist at the optimum de 
sign condition, wherein all elements buckle at the same 
stress. It might be thought that, because the web and 
flange elements tend to buckle in waves of different 
length, mutual restraint between adjacent elements 
might appreciably affect the value of k when the rela- 
tive size of the webs and flanges is changed (while keep- 
ing the relative values of ¢/b the same). However, it 
appears that there is little interaction between the webs 
and the flanges at the optimum design condition. 

Fig. 5 of reference 8 gives values of k for the individ- 
ual elements of a hollow rectangular cross section. 
For the condition of equal values of ¢/b for the sides 
and ends, the effective edge-restraint coefficient k was 
nearly equal to 4.0 for all of the individual elements 
over the entire applicable range of variables. This 
clearly indicates an effective end-restraint of simple 
support at the optimum design condition with little or 
no interaction between the ends and the sides, even 
though they have different wall thicknesses and/or 
lengths. 

To check the value of k for the other shapes, the fol- 
lowing method was used: First, values of k = 4.0 
and k = 0.525* were assumed for the webs and flanges, 


0.525 for outstanding flanges was deter 
2 By assuming that 


* The value of k = 
mined indirectly from Fig. 2 of reference 8. 
k = 4.0 for the web and then working backwards, it was possible 
to determine what the value of ke should be in order to remain 
constant while values of the plotted parameters bw/by and tw /ty 
were changed [and while keeping (tr/br)/(tw/bw) constant, as 
explained above]. That such a value existed is considered suf 
ficient justification for the assumptions made in the body of the 


paper regarding optimum design configurations. 
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respectively. This defined the relative values of t/b 
for the flat-plate elements as 


2.76 


(tr bp) (tw/bw) = V kw kr = 
for optimum proportions. Plotting lines representing 
this geometry on the charts of Figs. 1 through 4 of 
reference 8 showed that the effective values of kp and 
kw were practically equal to the assumed values over 
the entire applicable range of variables, thus indicating 
that at the optimum design condition the effect of inter- 
action between the individual elements was negligible, 
at least for purposes of this simplified analysis. 

The value of k used for the webs is the theoretical 
value for simply supported edges,’ while the value 
used for the flanges is somewhat higher than the theoret- 
ical value.* 

Since the sheet is supported along two edges with 
closely spaced rivets, the condition of simple edge sup- 
port has been assumed for the sheet in the example 
calculations. This gives a value of k, = 1.00. 
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On Similarity of Thermal Stresses 
in Elastic Bodies 


Martin Lessen 
Professor of Aeronautical Engineering, The Pennsylvania State 
College, State College, Pa. 


June 3, 1953 


F THERE EXISTED geometrically similar model and prototype, 
] what criterion would exist for a similitude of thermal stresses 
in model and prototype? 
consideration of the relevant dimensionless equations. 


This question can be answered by a 


The dimensional equations for an isotropic thermally conduct- 
ing elastic solid may be given in Cartesian form as follows: 
Equilibrium 

07;;/Ox; = 0 
Stress-Strain 


e;, = ((1 + v)/E) rij — (v/E) rex6i; + af T — Ty)6;; 


Energy 
p(OU/dt) = K(0?T/0x,0X») 
State 
v= % = AT — ty) + Me, ~ «,,) 

where 

Tij = stress tensor 

%j = space coordinate 

«ij = strain tensor 

v = Poisson’s ratio (isothermal ) 

E = Young’s modulus (isothermal ) 

a = thermal coefficient of expansion 

i i = temperature 

T) = reference (original) body temperature 

p = density 

U = internal energy per unit mass 

t = time coordinate 

K = thermal conductivity 

Uo = reference internal energy at 7) and e,;, 

Y = thermodynamic property corresponding to specific 

heat at constant strain 
ny = thermodynamic property defined later 
€kko = reference (original) strain taken as equal to zero 


The equilibrium equation neglects the effect of body forces; it 
is assumed that the inertia effects are small. The stress-strain 
equation contains the effects of temperature change. The energy 
equation is a statement of the first law of thermodynamics in 
which the work term is eliminated by linearization. The state 
equation is obtained from the linear terms of a Taylor series ex- 


pansion about the reference state as follows: 
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Contributions should not exceed 800 
ro 
U = U(T, «;) 
= one aU\T =h 
U = U+(T — TM) “aa + (€j; — &;,) 
a. = €ijo sii 3 O¢€;; €1)] =Cijy 
+ higher order terms 
Let 

QU/OT = y, OU /deij = ij 

If the material is isotropic, \;; = 6;;. Therefore, U — Ul, = 


y(T — To) + Mek — exo) 


The following scaling transformation is then made: 


Tij = Ea(T, — T)7;;' 
€ij = a(T, — 7 )e;;’ 
T— 7 = (1% — %%)T’ 

Xs = lx,’ 

t = ot’ 


where 7; = surrounding temperature that the body is exposed 


to, / = characteristic length scale, and o = characteristic time 
scale. 

If the state equation is substituted in the energy equation, 
the resulting dimensionless equations (with all primes dropped 


follow: Equilibrium 


O7,;/Ox; = OV 
Stress-Strain 
é:5 = (1 + v)riy — vreeds3 + 75; 
Energy 
oF... a Ko oT 
r + ex} = x 
ol ¥ pyl? OXm OXm 


The time scale, ¢, can be adjusted to keep the coefficient 
Ka/pyl? constant. The remaining conditions to be satisfied for 
similarity are that vy and A\a@/y be constant. 

For most metals in the elastic range, v is almost constant and 


* 


equal to about 0.3. It can be shown,* using the Maxwell rela- 


tions, that 


a E,-E l B 
¥ E, (1 — 2y) ¥ 
where E, = Young’s modulus (adiabatic) and 8 = specific heat 


at constant stress 
For aluminum (and most metals) 


ha/y = 0(0.01) 


It, therefore, seems reasonable to neglect the effect of \a/y in 
the energy equation and to arrive at the conclusion that simi- 

* Duke, C. E., Analysis of Damping in Elastic Solids, Doctorate Disser- 
tation, Aeronautical Engineering Department, The Pennsylvania State Col- 
lege, December, 1952. 
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ity of thermal stresses exists for elastic metal bodies of similar 


lar’ 


geomet ry. 
A conclusion from the foregoing considerations is that, if a 


model and prototype of the same material are submitted to the 
same thermal excitation, the corresponding thermal stresses will 
be the same but the time scale of the thermal stress history will 
yary as the square of the linear dimension scale. 


Discrete Vortex Systems in the Transition 
Range of Fully Developed Flow in a Pipe 


JR. Weske and A. H. Planthoilt 
The Johns Hopkins University and Consolidated Gas Electric Light 
and Power Company of Baltimore, Respectively 


June 1, 1953 


| pram RESULTS OF STUDIES of the transitions from fully 
developed laminar to turbulent flow in straight pipes (see, 
eg., reference 1, frontispiece) do not give a satisfactory insight 
into the details of the process by which the original unidirec- 
tional (peripheral) vorticity of laminar flow is augmented and 
transformed to create the omnidirectional vorticity distribution 
encountered in turbulent flow. Further study of this process, 
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undertaken by the junior author as a subject for an M.S. thesis, 
was influenced by the study of Theodorsen’s paper? and by con- 
versations with its author. Observations were made of the mo- 
tion of discrete vortex systems generated in the fluid which re- 
sults from interaction of these vortices and the surrounding rota- 
tional stream. It was recognized that the ‘‘Magnus effect" 
plays an important role in this interaction. When applied to an 
elementary vortex filament in rotational flow this effect may be 
defined by the Kutta-Joukowski theorem as follows for the two- 
dimensional case (Fig. 1): 


df = p(&r — &)Vrei. dA (1) 
where 
te = vorticity in the elementary cross-sectional area dA 
of the vortex filament 
fy = vorticity of the ambient field at the filament 
Vret. = (Ver — Vo), velocity of the vortex filament relative to 


the ambient fluid 
df = force of interaction between vortex filament and 
ambient fluid, normal to the vorticity and the ve- 


locity vector 


Ifte > &, df is positive. 
If tv < &, df is negative. 
If tr = &, the force is zero. 


The case Er = 0, & + O which merits special treatment will not 
be discussed here, since it does not pertain to a discrete vortex 
system. 

In a real fluid, a vortex moving relative to the ambient fluid 
experiences a drag in addition to the above lift force. Acceler 
ation of the vortex with respect to the ambient fluid tends to re- 
duce the difference (Vr — Vo) = Vrei. to zero. This tendency is 
counteracted if, by displacement of the vortex along a positive 
velocity gradient, | Vo| is increased. For conditions of a real fluid 
according to Fig. 1, displacement in the +-y direction should pro 
duce a larger Vrei. than a corresponding displacement in the op 
posite direction. 

Likewise, in a real fluid, the circulation dl! = trdA varies with 
time, and for constant dA, the vorticity &* tends to approach the 
value of the ambient vorticity &. If, for Er > &, the vortex fila- 
ment is displaced laterally in direction of a negative gradient of 
vorticity &, which is the case in Fig. 1 for displacement in the 
+y direction, then & also decreases with + and (fr — £) may re- 
main positive over a longer period. The rate of decay of vor- 
ticity of a discrete vortex filament relative to the rate of its dis 
placement along a gradient of vorticity of the ambient fluid 
clearly is a function of the Reynolds Number of the vortex fila 
ment. This relationship, however, will not be defined further 
here. Decrease of the cross-sectional area dA at constant circu 
lation d¢ of the vortex filament through stretching (see below) 
will also increase the product (te — £)dA, since the distributed 
vorticity £ of the ambient fluid is constant with variation of dA. 

To maintain relative motion between vortex filament and 
ambient fluid in a real fluid there must be a reaction force, in addi 
tion to D’Alembert’s inertia force, to counteract the aforesaid 
lift and drag forces. This reaction is supplied by the component 
in the x-y plane, Fig. 1, of the tension force in the vortex filament. 
Consequently, the vortex filament must be three-dimensional. 
For an ideal elementary vortex filament of radius r and of circu- 


lation ' = tprr?, Fig. 2, the tension force created by static pres- 
sure T = dA (Apmar./2) is 
T = (p/16xr)T? (2) 


which, for constant circulation along the filament remains the 
same irrespective of the amount of stretching imposed upon the 
filament. 

The work of stretching a vortex filament, disregarding fluid 
friction, equals the increase of the energy of the vortex filament 


2 


dT-Al = A |(pIp/2) (tr/2)*1) = (p/167) T2Al (3 
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where Jp = ar‘/2 is the polar moment of inertia, / is the length 
of the filament, and A/ is the elongation through stretching. 
This equation describes the important process of conversion of 
energy of mean motion of flow into energy of motion of discrete 
vortex filaments. 

In order to show that the above relationships are pertinent to 
the process of converting the unidirectional vorticity distribution 
that exists in iaminar flow into the omnidirectional vorticity 
distribution of turbulent flow, the motion of discrete vortex fila- 
ments was observed and recorded experimentally in the following 
manner. 

The experiments were made with water in a straight pipe of 
1°/s in. inside diameter in a section made of a clear plastic tube, 
150 diameters downstream of a well-rounded inlet, this length 
being considered sufficient to produce fully developed flow 
Disturbances were minimized by steadying the flow in an equal 
izing chamber ahead of the inlet. 
from 1,600 to 2,300. 
while at 


Reynolds Number was varied 
Below this range vortices decayed quickly, 
higher Reynolds 
interfered with observations. 


Numbers extraneous disturbances 
The formation of vortex filaments was induced by a trip 
namely, a dise placed adjacent to the wall of the tube and normal 


to the stream. Three trip designs were used: 


(a) A low aspect ratio trip, '/, by '/, in., the outer edge 
rounded to '/s in. radius. 
(b) A higher aspect ratio trip, '/2 by */i,in., outer edge gently 


curved. 
(c) A “scimitar” trip, '/is in. uniform width, curved camber 
,in., °/s in. from wall to wall, such that fluid could pass under 


neath the arch formed by it. 


The vortices were made visible by allowing ink to ooze from the 
sharp outer edge of the trip at a slow but steady rate or by bond- 
ing tellurium powder electroconductively to the upstream surface 
of the trip and rendering the fluid motion visible through elec 
trochemical action, by the new method developed by Wortmann.* 

It was found that a succession of horseshoe vortices were 
formed in the velocity discontinuity surface downstream of the 
trip. Generally these vortices appeared at a definite frequency 
that varied with Reynolds Number and perhaps with other factors 
not as yet identified. These vortices were stretched in down 
stream direction at a rate slightly less than the speed of the am 
bient fluid, and at the same time they were displaced across the 
stream toward the center of the pipe. Fig. 3 shows a series of 
such vortices photographed in a light flash of 30 microsec. dura- 
tion, at right angles to the plane of the trip. Only on exceptional 
occasions did the tips of horseshoe vortices penetrate beyond the 
centerline of the tube, this line being the limit up to which an 
excess of velocity of the ambient fluid above the velocity of the 
vortex filament could be maintained. 

Fig. 4, taken of a similar train of vortices from the same trip 
at 45 deg. with respect to the plane of the trip, shows more 
clearly the loop shape of the vortex filament. 


setup did not lend itself to the observation of the motion of ran- 


The experimental 
domly generated vortex filaments. Therefore the present evi 
dence merely indicates that artificially induced vortex filaments, 
originating from a finite disturbance at Reynolds Numbers in the 
transition range, by the action of the “Magnus” effect, aug- 
menting existing vorticity, generate vorticity components in both 
radial and axial directions. 

An interesting corollary of these observations is the fact that 
vortex chains of the horseshoe type in a shear flow no longer coin 
cide with stream surfaces but are actually pervious to flow. It 
was of interest to ascertain whether penetration of ambient flow 
into a chain of vortices like that of a von K4rm4n vortex street 
takes place. 

Fig. 5, taken at R = 2,200, shows the two rows of vortices of 
opposite vorticity produced by the 
only the center portion was covered with tellurium. 


“scimitar” trip, of which 
The 
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Fic. 3. Horseshoe vortices at R = 2,000 photographed at right 


angles to the plane of the trip. 





Fic. 4. Horseshoe vortices at R = 2,000 photographed at 45 


deg. to the plane of the trip. 





‘scimitar” 


Divergent vortex street downstream of the 


Fic. 5 
trip, at R = 2,200 
tices of the outer row appear to behave like the ones shown in 
Fig. 3, the vortices of the inner row move at lesser speed down 
stream and toward the wall and decay more rapidly. There ts 
however, clear evidence that the two vortex rows diverge in shear 
flow and that fluid from outside the original discontinuity stream 
Both the 


divergence of vortex loops and penetration of external fluid should 


surfaces enters the region between the vortex rows 


contribute to the rapid decay of the wake behind the obstacle 
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Variation of Transition Reynolds Number with 
Mach Number 


A. H. Lange, L. P. Gieseler, and R. E. Lee 
U.S. Naval Ordnance Laboratory, White Oak, Silver Spring, Md 
June 15, 1953 


ad | Maz REYNOLDS NUMBER OF TRANSITION of the boundary layer 
on a slender cone and on a hollow cylinder has been deter 
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Fic. 1. Boundary-layer transition in the NOL Aeroballistics 


Wind Tunnel No. 2 at various Mach Numbers. 


mined from spark-schlieren photographs at Mach Numbers be 
tween 1.86 and 5.00 in the intermittent NOL 40- by 40-cm 
Aeroballistics Wind Tunnel No. 2. The cone had a total apex 
angle of 5 deg., a length of 20 in., and a base diameter of 1.75 
in. The hollow cylinder had an outside diameter of 4 in., a 
length of 19 in., and an internal bevel of 10 deg. Supersonic 
flow was entering the hollow cylinder (swallowed normal shock) 
Both models have constant pressure along the surface in an ideal 
flow. They were nearly at equilibrium temperature, 1 < Twa 
Tequil- 1.07, during the test. The Reynolds Number was 
computed from the free-stream density, velocity, and viscosity, 
the latter being determined at the free-stream temperature 
The length used was the distance from the tip of the models to the 
place where the boundary layer became turbulent 

Results are shown in Fig. 1. Each point represents an average 
A decrease in Reynolds Number 
This 
trend agrees with stability theory! but is opposite to that indi 
Further results and 


of at least twelve observations 
of transition with increasing Mach Number was observed. 


cated by Gazley? in Fig. 9 of his article 
additional experimental details are given in NAVORD Reports 


2752’ and 2823,‘ which have been published recently. 
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Measurement of Boundary-Layer 


Boundary Layers with Variable Heat Capacity 
on Nonisothermal Surfaces 


Martin Bloom 

Assistant Professor of Aeronautical Engineering, Polytechnic 
Institute of Brooklyn, Brooklyn, N.Y 

July 29,1959 


HE ANALYSIS OF CHAPMAN AND RuBesIN! for laminar bound 
ary-layer flow over a nonisothermal flat surface can be modi 
fied in a simple way to take into account the temperature de 
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pendence of cy (see reference | for notation). The energy Eq. 6' 


for constant Pr can be written as follows: 


oi oi L @ ( *) (*)’ +) 
uu ~ v = rm + (1) 
e Ox “ oy (Pr) oy Ov 4 oy 


where 


the enthalpy fuuction 
The transformed form of this equation, analogous to Eq. 20,' is 


O*s* oi* o7* (Pr) u,? 
+ (Pr)f APr)i'<* = (f’)? (2) 
On? On Ox* 4 do 
where 
g* =m 4/9 


The solution of Eq. (2) and the boundary conditions satisfied 
by the solution are given as follows: 


a*(x*, nn) = 1 + (42/2 @)r(n) 4 ) a,x*”" Y,(n) (3) 


n=0 

iy* = 1 + (u,2/272)r(0) 4 ) a,x™ (4) 
n=0 

e*(x*, ©) = ] (5) 


The functions r and Y, are given in Figs. 4 and 5.!. The trans 


formation of Eq. (41)! remains unchanged —viz., 
. en 7 
(y/2) (u./v¥.xC) i T* dn (6) 
J0 


However, the integral on the right-hand side of Eq. (6) must 
be evaluated by means of tables or formulas relating i and T 
Eq. (42)! is not valid when c, depends upon T 

The rate of heat transfer is given as follows: 


(27) = Ss (=7=)/ Ua ) 
g= —k, = — 
OV Ju 2Cp« ee VorC 


n=0 


enn"Y.00) (7) 


where coefficients a, are determined by the prescribed variation 
of i,* in Eq. (4). Eq. (7) is almost identical to Eq. (50).! For 
low values of 7x, it is seen that 
to = Cpo Ta 
Uat/ia = (yo — 1)Mo? 

and Eq. (7) is identical in form to Eq. (50),' although the a,’s 
are defined differently. Eq. (57)' for the shear remains un 
changed 
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Concluding Discussion on the Hydraulic 
Analogy 


E. V. Laitone 
Associate Professor, University of California at Berkeley 


June 25, 1953 


ig COMMENTING ON MY PAPER,' Bruman? is correct in stating 
that the surface wave group velocity has no relation with the 
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strong hydraulic jump. However, he is completely wrong, and 
has entirely missed the main point of the paper, when he says 
that ‘‘experiments with the hydraulic analogy are usually de- 
signed to represent shock waves by hydraulic jumps.’ In the 
first place, practically all the analogy experiments are more con- 
cerned with the potential or isentropic flow which is related only 
to the group velocity. As a matter of fact, even the papers he 
refers to, by Orlin? and Preiswerk,* are chiefly on potential flow 
with no hydraulic jumps present. The Orlin’ paper deals with 
the analogy of subsonic flow about a circular cylinder; the 
Preiswerk‘ paper gives the isentropic supersonic expansion in a 
Laval nozzle. Preiswerk‘ states (Part II, page 49) that water 
depths above 1 in. gave increasing errors as compared to theory 
because the vertical acceleration was no longer negligible for the 
larger water depths. (This corresponds to the group velocity 
varying with wave length as I pointed out in reference 1.) The 
agreement shown in Fig. 78 of Preiswerk’s‘ paper would have been 
much better if he could have used '/;-in. depths. 
the water was flowing, the boundary-layer growth on the tank 
bottom made this ideal situation impossible. 


However, since 


Fig. 14 of Orlin’s* paper clearly shows the deleterious effect of 
water depths greater than '/,-in. The deeper water depths show 
a strong detached shock wave at supposedly subsonic speeds. 
This is not a capillarity effect, since it disappears for shallow 
water. It is the effect of a varying group velocity. Part of the 
flow pattern corresponds to supersonic flow, even though most of 
the flow field corresponds to subsonic flow. 

The supersonic flow pattern is smaller and rides over the top 
of the subsonic flow pattern and destroys the analogy. It is pro- 
duced by the deeper water having a much lower group velocity 
than given by V gh for the shorter wave lengths (see reference 1); 
consequently, the shorter wave lengths all produce a common 
supersonic flow pattern. Again, the agreement would be best 
with '/,-in. water depth, but the boundary layer below the flow- 
ing water caused a more serious error. In all cases the experi- 
mental data are not in good agreement with wind-tunnel tests 
over the front portion of a cylinder because the water depths 


were greater than '/, in. - 


Another point that still seems to need further emphasis is that 
the group velocity, not the wave velocity, provides the only 
analogy to the speed of sound, since it is the actual rate of prop- 
agation of energy and disturbance (see reference 1). For ex- 
ample, the Mach or characteristic lines in supersonic flow are 
the envelopes of the infinitesimal disturbances traveling at the 
local speed of sound, which is not a constant but depends on the 
or enve- 


‘ ” 


local temperature. These same ‘‘characteristic lines 
lopes are formed in free-water surface flow by the group velocity 
that must vary with the local water depth as V gh but should not 
vary with wave length in order to maintain any analogy. The 
only possibility of maintaining a nearly constant group velocity 
for the entire spectrum of wave lengths produced is to use as 
large a model as possible and try to keep all the water depths be- 
tween the limits of 0.10 to 0.25 in. for standard surface tension. 
This even requires that, if strong compression regions are present, 
even though isentropic, the free-stream water depth must be less 
than '/, in. so it will not rise above this value in the strong com- 
pression region. If the above conditions are safis‘ied, then all 
potential disturbances, except the negligible so-called capillary 
ripples, are propagated at the local group velocity of \ ‘oh, and the 
flow corresponds to the potential flow of a perfect gas with a spe- 
cific heat ratio of 2. Weak hydraulic jumps may be present, 
since they have a negligible increase in entropy and therefore do 
not invalidate the analogy. Strong hydraulic jumps cannot be 
present, since they do not correspond to a shock wave in a perfect 
gas (see Preiswerk,* Part II, page 3). The energy equation 
cannot even be used in this case, since there is a loss of total head 
corresponding to a change in stagnation enthalpy. However, 
certain comparisons can be made with strong hydraulic jumps if 
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there are no regions corresponding to isentropic expansion or com 


pression. These are mainly one-dimensional problems. 
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Open-Nose Body of Revolution of Minimum 
Drag at Small Supersonic Speeds* 


Elliot T. Benedikt 
Consolidated Vultee Aircraft Corporation, San Diego, Calif 
June 11, 1953 


ABSTRACT 


The drag coefficient of an open-nose body of revolution is approximated 
by means of a formula, valid at small supersonic speeds, due to Lighthill 
The variational conditions that the equation of the profile of the body must 
satisfy in order to minimize the drag are converted into a singular integral 
equation of the Hilbert type. This equation is solved by means of a new 
procedure of G. Tricomi The extremizing profile thus obtained is shown 
to correspond to an actual minimum for the drag 


SYMBOLS 


Cp = drag coefficient 

x (or y); R cylindrical coordinates (see Fig. 1) 
Ro = R(0) 

l length of profile 

E = x/l 

Ri = Ril) 

Ss S(x) cross-sectional area of body 
So = S(O) 

S; = Sil 


* 7b \ fc-e *b / 
f f(x) dx lim -« f(x) dx + J f(x) dx? 
a p->0 (Ja cTp \ 


principal value of integral of function with singularity at 


M = Mach Number 
a V M?- 1 
€ e(x) arbitrary function [variation of S(x) 
a, b, 
C, co 
C1, C2, ¢ constants of integration 


INTRODUCTION 
A MERIDIONAL CROSS SECTION of the open-nose body of revo- 
lution considered here is shown in Fig. 1. The body will 


* This paper is based on work carried out in January, 1951, for the Fort 


Worth Division of Consolidated Vultee Aircraft Corporation 
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READERS’ FORUM 


to x = © being cylindrical with radius 
0 and having, from the last point on, 


Making use of an approximate 


= o, and to be cylindrical tending from x = — « 


be ia 
with radius Ry from the cross section x =/ on. The tangent Ro from x = —& tox = 
ee . . the shape described above. 
planes at the orifice and at the point x = / will be assumed to - ’ ; : : 
form given by Lighthill' valid for slender bodies of revolution 
having a discontinuous profile and moving at small supersonic 


assumed to extend from x = O tox 


be parallel to the axis of symmetry of the body. The flow will 
R*S Ro. speed, the external drag coefficient can be expressed in our case 


by the formula (applicable as long as the change in diameter 


of the body is sufficiently small), 


be further assumed to be uniform in the region 0 ‘ 
The last condition implies that the flow, and consequently the 


drag, will be the same as the one that obtains with a body ex- 


+) > ’ *) 
Cp = ("/27) log S"(x)S"(y) dx dy — 2S'(0) S"(x) log - dx + [S’(0)]? log (2/aRo) { 1) 
0 /0 eae Jo x 
The function S(x) corresponding to the profile of equation 
, Fe (2) 
R(x) = VV S(x)/2 « 
and minimizing the above expression will now be sought. 
FORMULATION OF THE VARIATIONAL PROBLEM IN INTEGRAL FORM 
Consider a variation ACp in Cp caused by a bounded variation ¢«(x) of S(x) satisfying the boundary conditions 
(3) 


ex) = €'(0) = ei) = e(l) = 0 
We shall have 


"s ™ 
27ACpy = / | log [S”’(x) + €"(x)}[S’(y) + €"(y) ldx dy 
/ 0 0 x-y 
+) 


2S8'(0) [S"(x) + €"(x)] log —dx + [.S’(0)]? log (2/aRo) -—2rCp 
/% x 
"1 ~ 
> log e"(x)S"(y) + €"(y)S"(x) |dx dy — 
Ft oF x ¥y 
* 7 fl 
28'(0) e"(x) log — dx 4+ log e"(x)e"(y)dx dy 
0 x J 0 0 x y 


or, because of the symmetry of the first integral, 


7] *) */ “1 i 
. ‘ " l ” ” l I ” ” 
27ACp = 2 / e”"(x)S"(y) log dx dy — 25'(0) e’(x) log - dx 4 | log e"(x)e"(y) dx dy 
J? 0 el 0 x 0 0 x y 


This expression can be transformed by repeated integration by parts, and taking into account conditions (3) into 


? ne . ™] aa | 
d? : 1 so) | l 
rACp = (x) ) S”(y) log dy dx + 3 / e(x)e(y) dx dy (4) 
; y x? f J 0 (x y)4 


J 0 { dx? 0 x 
The first integral converges since, because of the second of conditions (3), 
e(x) = O(x?) (5) 
forx—>0O. If € is small (!e!? ~ 0), expression (3) becomes the first variation 
-) ? ) ; 
sC : J} a , so) 
rbCpy) = e(x) - S”(y) log dy { dx 
UL dx? So «—y x? 


0 
Its vanishing is a necessary condition for the value of Cp to be a minimum. Due to the arbitrariness [apart from conditions 


of Cp. 
(3) (4)] of the function e(x), the above requirements will be satisfied only if 

1? > l S’(0 

a* ” . } 

S"(y) log dy = () 

dx? Jo [*t-y = 

Integrating twice, we obtain the singular integral equation 
7) 
S"(y) log dy = — S'(0) logx +ax+) (6) 
/0 ~~ =F 


which the function S’(x) must satisfy in order to minimize Cp. It is convenient to further modify the integral equation by an inte 


gration by parts. We obtain, thus, 
’ I }| 1 | \ {' si 
/ S"(y) log dy = lim 4] S’(y) log . + | S’(y) log { “ — dy 
/J0 a ~~ Fi |X — Yi Jjx+o , e~'?2 


om 











Imposing upon S(x) the not too restrictive condition 


lim [S’(x + p) — S’(x — p)| log p = 0 


p—0 


and because of the assumption 


S(0) = S'(l) = 0 (7) 


we have 


1 
S’(y) 
dy = (3 ax 
3—-x 


S”(y) log 
0 Ix at | 0 


Replacing this result into (6) we obtain the integral equation 
*/*] 
{/ S’(y) 


F dy =ax+b 


(8) 
aS 7 =F 


for S’. 


INVERSION OF THE INTEGRAL EQUATION 


The solution of the integral equation (8), if it exists, can be 
expressed by means of an inversion formula established recently 


1 i Val — y 


l ) 
C- — (ay + b)dy 
rV x(1 — x) a 0 y-x a . 


by Tricomi? in the form, 


S"(x) = 
(9) 
By means of the substitutions 
y = (1/2)(1 — cos 8) 
x = (1/2)(1 — cos ¢) 
the integral 


* l J ? - 
‘2-2 (ay + b) 
y-x 


J = 
T 0 
appearing in Eq. (9) can be written in the form 
* ® P 
1 sin? 6 
[ = = [4A cos 6 — B] dé 
7 Jo cosé?—cos¢ 
Using the identities, 
sin? @ = (1 — cos 26)/2 
sin? 6 cos 8 = (cos 6 — cos 30)/4 


and 


cos n6 sin ng 
= FT 


de 


0 cos@?@—cos¢ sin ¢ 


I becomes a quadratic form in cos ¢ and consequently in x. 
Replacing this result into Eq. (9), we obtain 
Co + Ox + Cox? 


/ 


x(p — x) 


S(x) = (10) 
Because of Eq. (7), we must have co = 0, c,/1 = —co=c. Re 
placing this determination of the constants of integration into 
Eq. (10) and integrating once more, we get 


i 
¢ 


S(x) = So +c Vv “yl lL —y)dy = 
0} ‘5 
So + (¢ 4)[(2x —/1)Vxil — x) —Fsin7'! Vx 1] (11) 
From the boundary condition S(/) = S|, we obtain the deter- 
mination 
c/4 = 2(S,; — So)/l?x = 2(R,? — Ro?)/l? (i¥*) 


for the remaining integration constant. 


THE PROFILE OF MINIMUM DRAG 
The expression 
’ ie = 
—-1)V x1 -2¢)] 


(12) 


S(x) = Sy) + (2/r)(S, — Sy) [sin-' Ve + (22 





° e 3 . 1.0 
X/L 
Fic. 2. Minimum drag profiles for Ro/R; = 0, 0.25, 0.50, and 
0.75 
iB 


(L/2R2) Cp 





0 2 4 6 8 1.0 
Ry R, 


Fic. 8. Minimum drag coeflicient as a function of R)/R;, 


obtained replacing Eq. (11’) into Eq. (11) represents the cross- 
sectional area of a surface of external drag. This result has to 
be further qualified by an analysis of the second variation which, 


as shown by Eq. (4), is 


1 fl 
iain (x )ey) 
(3/3)8Cp = aeatl 
0 o (x — y)* 


Because of the assumed boundedness of the function ¢«(x), the 


dx dy 


largest contribution to this integral will arise for values in the 
neighborhood of x = y. The numerator of the integrand can 
then be replaced by e*(x), from which it follows that 6°Cp = 
o>. 
fore minimize the drag of the body of revolution considered here 
within the range of validity of formula (1). Because of Eq. (2), 
the profile corresponding to this cross section will be given by 


R(x) = Ro 1 + (2/7)[(Ri/Ro)? — 1] 


[sin-?§ Vt + (2 —1) Vet — £)]} 


The cross-sectional area defined by Eq. (12) will there- 
] 


Graphical representations of this function are shown for various 
values of the ratio R:/Rp in Fig. 2. The dependence of the mini- 
mum drag upon the ratio between the rear-end and front-end 
base diameters of the profile is represented graphically in Fig. 3 
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SCALE OF TURBULENCE, s= “/, 


An Expression for the Gust Alleviation Factor 


M.L. Williams 
Assistant Professor of Aeronautics, 
Technology, Pasadena, Calif 


June 22, 1953 


California Institute of 


N A RECENT PAPER,' Fung has investigated the application of 
statistical analysis to structures loaded dynamically. As one 
example of the procedure, he has worked out the root mean square 


vertical acceleration of a rigid two-dimensional wing in turbulent 


flow as 
52 = w? [4U2/rb1 + «)?]I(a, s) 
where 
mean square vertical acceleration, ft. per sec. per sec 
w? = mean square turbulent vertical velocity, ft. per sec 
l = flight speed, ft. per sec 
b = semichord, ft 
k = mass density ratio for the airplane (=2W/mpShe 


W’/S = wing loading, Ibs. per sq.ft 


a = 2/[s(1 + x)] 

S = b/L 

L = characteristic length of turbulence, ft 
p = air density, slugs per cu.ft. 


= gravity acceleration, ft. per sec. per sec 

I(a, s) = function given graphically (Fig. 5) and also analyti 
cally [Eqs. (56)-(60)] in the reference 

By making use of the approximation that x >> 1 in most 


practical cases, the relation can be manipulated into the form 


(W/g)N 2 (24)p Vw? US I(a, s) 
W 2W 


which, upon defining (W/g) (s?) W as the mean acceleration, 
An, (w?)/? by V, a mean vertical velocity, and replacing 27 by 


m, the lift curve slope, there results the form 
pUVSm _{I(a, s) 
An = . 
2u 7 


is equated 


which is similar to the usual gust formula if (//7) 
? 


to A, the gust alleviation factor. Fig. 1 shows a graph of (//7) 
against the chord to turbulence ratio with the mass ratio as a 
parameter. It may be noted that in American literature it has 
been customary to use K = K (W/S), although other countries 
have used the mass ratio dependence. The two differ essentially 
by the length of the chord, but it is thought that the mass ratio 
is the more definitive parameter. 

Until experimental evidence for the scale of atmospheric tur 
bulence is complete—i.e., the proper number of chord lengths to 
use for L, the scale or characteristic length of the gust, in obtaining 
s and thus A 
which vary as a function of the mass ratio; secondly, reference to 


it is suggested that one use the maximum values 


the figure gives a good idea as to the sensitivity of the incremental 
load factor as affected by the mass ratio and turbulence scale. 

It is believed therefore that the Fung expression may be used 
as a basis for separating the gust velocity ( V), which is a function 
of the atmosphere, from the gust alleviation factor (K), which 
characterizes the airplane, thus permitting a quantitative pre 
diction of the gust load factor. Improvements will be effected by 
considering (1) the probability distribution of An, (2) a nonrigid 
wing, (3) finite span effects, and (4) a more precise representation 


of the gust power spectrum. 
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| REFERENCE 1, A SOURCE-SUPERPOSITION METHOD of solution 
of a periodically oscillating wing at supersonic speed was 
proposed by the writers. Using this method, the problem of a 
rectangular wing with supersonic periodic motions was treated.* * 
Miles,‘ by a Fourier transform method, Stewartson,® by the La- 
place transform method, and recently Goodman,® by a method 
proposed by Gardner,’ have also treated this rectangular oscil- 
lating wing problem. While our solutions of the rectangular 
wing yield numerical results in close agreement with the results 
of Miles, et al., over a wide range of frequencies and Mach Num 
bers, it has been discovered that they are not identical. In solv- 
ing the basic integral equation [Eq. (25), reference 1] which 
appears in our method of approach, we have made an implicit 
assumption that the magnitude of the upwash velocity in the re 
gion outside the wing tip is, in a certain sense, independent of the 
oscillatory motion. This assumption imposes severe limitations 
on the solutions such that they are, as will be shown later, valid 
in general only for terms of zero and first order in frequency 
As the general method of approach outlined in reference 1 can be 
used in treating problems of a class of wings of general plan forms, 
it seems worth while to make some attempt to improve the ac- 
In the present note, the basic integral 
A procedure has been found whereby 


curacy of our method. 
equation is re-examined. 
it is possible to improve the accuracy of the method to yield solu 
tions valid to terms of any order in frequency. 

t Since the completion of the present note, two interesting papers!” !! 
have appeared which deal with problems related to the present subject by 





The basic integral equation under consideration is 


arn 
“D du "= Ap — Ar 
; ; - F(u, 0, v, Up, Up) dv = 
J 0 Vup—tJn 2Vrp—-v 
.. 
“uD du vp Ap 
J F(u, Vv, v, UD, Vp) dv (] 
0 Vup—u YW tp —-? 
where 
J ip | 
F = exp 4 — [((up — u) + (vp — v)] 7 X 
| Ba f 
2v , 
cos V (up u) (vp — v)} (2) 
BaM 


where f? 
a represents the speed of sound in the free stream, and » represents 
The quan- 


= MM? — 1, M represents the free-stream Mach Number, 


the frequency of the periodic motions of the wing. 
tities 1), v2, which represent the leading edges of the wing (Fig. 1), 
are given functions of u; Ag and Az, which represent the effective 
slopes of the stream lines on the bottom and the top sides of the 
wing, are given functions of the characteristic coordinates (u, v): 
up and vp define a point on the diaphragm region outside the wing 
tip; and Ap represents the unknown upwash function in the 
diaphragm region. 
the simple procedure in reference 1 under the assumption that 
In general, Ap is an unknown 


This integral equation for Ap was solved by 


Ap is a function of (u, v) only. 
function of (u, v, v, 7), since the physical parameters v and 
enter the kernel of the integral equation F as in Eq. (2). Sup 


pose Ap is expanded in an infinite series in v, thus, 


+ (tv/Ba)" A, (3 


n 0 


++ 


Ap = 
where » are integers. Substituting Eq. (3) into Eq. (1) and equat 
ing like power terms in v yields a system of integral equations 


with A,(u, v) as the unknowns. Thus, 


t Terms involving fractional powers and logarithms of v do not enter 


different approaches the solution because of the nature of Eq. (1) 
1 {> du ™ Aw — Ar { $2 i 
/ 4| (up u) + (vp = V (Up — tu) (vp — 0) T 
' / 
nN: 0 Vup—Uudn 2V vp v | M 
' 2 / "| “ (—)m “up du "Up A, 
(up — u) + (vn — 0) + —V (up — 2) (vp — ¥) { dy = / Xx 
vi m=0(n — m)! Jo Y/ tp —- &, Vvup-—1 
j 2 iis 2, ies. 
| (up — u) + (wo — 2) — V (up — u) (vp — 2) + | (up — u) + (vo — 2) 4 V (tp — u) (tp — v) ( dv (4 
| M M 
From Eq. (4), the following equations defining Ao, Ai, A2, and A;, respectively, can be derived: 
“Ay — Aq ""p No 
; dv = dv (5 
e 2V vp v : Vp? 
An — Az - oa Ai . 
V vp v di Ao V up ydv = - dv (6 
; 2 J Jn Vo -' 
‘Ap — Ar . “D 
Up v)/?* dv — Ao(vp — v)/*?*dv +2 
9 
dv (7) 
1 Ap — Aq ; 1 "D 
< (vp — v) ’? dv + Agivp — v) /? dv — 
6 J, 2 6 J» 
9 me “?'n he | 7 
diy dv + du; 
M2 Jo : u) Vito —v M? Jo ‘ 
dv (8) 
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similar equations define the A,’s for m = 4, 5,... The scheme of solution of this set of integral equations is: (1) Solve Eq. (5) for 
and solve for Ao; 


Ao; (2) substitute Ao thus found into Eq. (6) and solve for Ai; (3) substitute Ao, Ai thus found into Eq. (7), 
It is clear that we can carry out this scheme indefinitely to obtain A, for arbitrary n. It is further op that the A,’s are all 


etc. 
For a flat plate wing, (1/2) (Ans — Ar) = a, we 


defined by integral equations of the Abelian type whose solutions are well known 


obtain the following results for Ao, Ai, and Ag: 


2a lyo(u) — v(u) v(u) — 2; #) 
Ao = ( PAO) pat 5 [AED — A? (9)t 
rLN v — vu) } vy — v(u) 
A, = —(2a/37) } [v2() — v(u)] 7? V0 a ve(u) f (10) 
3/2 5/. 
= ‘ 2 3 
= < v(u) — v4(2) u(u) — —v(u)} — v(u) — v4(2) + 
QrV iv — v(t) | 3 5 
1 8/s 8/o 
4 s & 1 
32 vo(u) — 2(%) du + ve(u) — m(u) 20 — ve(u) (S$ (11) 
7 0 3 { 


At the subsonic leading edge v = v2(u), Ao, A1, and Ay become infinite; thus, they possess the correct singularity property. 
The contribution of the diaphragm to the velocity potential at a point («,, 7, +0) on the top or bottom surface of the wing is given 


1 





as 
a . 
Sunlts, Ve, 0,8) = F exp(irt) aie V, V, Uw, Vy) dv (12) 
T J 0 
where u’ is defined by v2(u’) = v, (Fig. 1). Expanding ®,p in power series in v yields 
x F n Un (-—) u’ 
a v 
Pp = F exp(ivt) 7. _ : y, 
Qq n=o\ Bad ma=o™ —™)! Jo Van 
» : n—m . 9 ; n—m) 
(te. — &) + (%e — 9) — uy © (He — u)(V» — v) + 1 (Ue — Ut) + (Vx v) + uy ¥ SM - u) (Vm — v) co (13 


We can rewrite our integral equations defining the A,’s, for instance, Eqs. (5)-(8), with vp replaced by v,, and these equations remain 


valid. Using these results, we can eliminate the A,’s in Eq. (13). The result can be written as 


a du ” Ap — Az — iy \" { 
Pip = F —exp(irt) 7 F(u, 0, v, Uw, Vw) dv + ? ~~ APn | (14) 
7 ' 0 \/ Ue —~UJn 2V0w —-? n=O Ba 
On the basis of the assumption Ap = Ap(u, v) we previously obtained the result that A®, = 0 for all » [Eq. (33), reference 1]. The 
corrected computations show, however, aes A®, = O is true only for » = 0 and » = 1 and that A®,, ~ 0 forn = 2, 3, This 


result substantiates our statement in the opening paragraph that our results in references 1-3 are valid in general for terms of zero and 


first order in frequency. For higher order of accuracy, we have to introduce the correction potentials A®,, which can be evaluated 


by our present method; for instance, 





= - ve Aj 
At, = Vu, — un’ du 7 v (15) 
Ma Jo oJ v2(u) V%-? 
2 u“ Un Ag u 
A, = Vly — u’'| —2 du dy + du 
VW? 0 is niu) Ve — v 0 
” Ai 
du dy 7 ; dv (16] 
Vv J(u) V Vw — v 
+ First obtained by Evvard and Turner in reference 
As an example of the application of the present method, we can __ rection term is 
C. C, 2C. 
ty\ _ ou a Mo (19) 


deduce some correction terms for the lift coefficient expression Cio — 
given in reference 2 for a rectangular wing with vertical oscilla- Li(ke) = BAR. ( Cro* 
tions. For this simple geometrical configuration, A’; and A®; 
can be simply evaluated from Eqs. (15) and (16). Thus, the lift 
coefficient in the mixed supersonic region, Cz,, as given in Eq 


C1,* 
where A.R. denotes the aspect ratio, and Cz, and Cy, denote the 
lift and moment coefficients in the purely supersonic region and 


(32) of reference 2, has to be subjected to the following correc Cr,* = (4a/B) exp(ivt) 
tions: 
The expansion of L;(k.) in reference 2 as a power series of k,, up to 
~~ e a k,? 7 the &.3-term, is 
As = correction in v? term = — exp (ivt) (17) 
12 8M? 1 i 1 1 7 
Lie) = — — +  . ke + l + k,* - 
1 1 k.3 a M? 24 2M? 
A; = correction in v* term = ia( +- exp (ivt) (18) = ; 
20 60M?/8M? t 15 3 F 
1 + + ke+... (20) 
120 2M? 2M‘ 


where k. = Myvx./82a is the reduced frequency, x, being the wing 


chord. From reference 2, [Eqs. (68) and (76), ~ wing-tip cor With the correction terms from Eqs. (18) and (17), 
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which is identical with the corresponding series from Miles’ re- 
sult.4 ; 
can be improved to an arbitrary degree of exactness. 
parent that the present method does not yield solutions in a 
closed form, so the particular problems of the rectangular wing 
are best treated by the methods of Miles and Stewartson. More 
complex plan forms can be treated by the present method or by 
that of Goodman. 


Therefore, by the present method, the previous results?» 
It is ap- 
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HE DETERMINATION OF THE STRESSES in a highly redundant 
structure inevitably requires a large number of numerical 
calculations, and in the design of such a structure these calcula 
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tions must be repeated many times to determine the optimuy 


combination of structural components. Some of these com 
binations may be identical except for modifications at a fer 
points in the structure, so that it would be advantageous to evaly 
ate the effect of these modifications without the necessity of re where [6 
and colu 
bAmm] a 
the elem¢ 


rewritten 


peatedly solving the entire system of equations. In a recent 
article! attention has been drawn to a simple method, originally 
due to Best,? which in some instances greatly decreases the 
amount of numerical work required. The purpose of this artic le 
is to generalize the method and to express it in a mathematica 


language that may be applied to any problem requiring the soly 


tion of a system of linear equations. The vec! 
In the method it is assumed that a solution has been obtained — 

for one combination of structural components and that a struc symbol, 

tural modification gives rise to changes in only a few of the co 

efficients of the unknown quantities. The equations of th The si 

modified structure are then manipulated so as to restore the co terminat 


efficients to their original values leaving the effect of the struc imu! 

on ¥ ‘ n, § 

tural modification in the constant terms. An elementary appli 
¥ app rows anc 


of [5Um| 


cation of this idea has been employed in electrical engineering for 
many years where it has been called the ‘‘Compensation Theo 
rem’’;? electrical engineers, however, have not taken full ad 


vantage of its possibilities because they are seldom interested i Hence, t 





accurate, detailed solutions of large networks 

Any method of analyzing redundant structures eventually 
Eq. (9) 
once [6 


duces to the solution of a system of simultaneous linear algebrai 
equations, and these equations will be taken as the starting point 


of the method only un! 


The system of 7 equations for the modified structure, writter [FP] + 
in matrix form are fy 
[A] [F] = [U 
and for the original unmodified structure | oe 
are reqt 
[A°] [F°] = [U] (2 be the | 
I, | be 
The n-rowed column matrices [F] and [F®] are the unknown “a 
quantities in the modified and unmodified structures, respectively 
and the m-rowed column matrix [U’] is the matrix of constant 
terms, assumed to be the same for both the modified and unmodi The : 
fied structures. determ 


column 


F, F,° U 
F F.° 7 Eqs. (1 
fi of [A°] 
[F] = , Pj=| | U] = will be 
n and 
Fy Fy 0 ie easily 
: Pabe involvi 
The » by ™ square matrix of coefficients [4] is assumed tobe} jg eyig 


the same as the matrix [A°], except for differences occurring 1 re 


the first m rows and the first m columns. It is therefore con 


: ae , “ ‘ Som 
venient to partition [A] and [A°] into four submatrices, respec case W 
al calles TR PO 

tively, of which three are pairwise identical, may b 
LA A mm'A mn ; , a) ee In this 

[A] = red Fa 1°] = y ty FS \o tion c: 

where 
ha.~- Aes ee a The 
S csie't = ° q | a 1 = 
Fro 
Am.1 r. faates 1 1 | ¥ 
As, mts di, | In < 
indeps 
Amn] = - ; , etc If 
v 
As at: + « Ae a 
The matrices [F] and [A] will be written as the sum of the un If ¢ 
inodified matrices plus correction terms portic 
iF] F F ‘ 15) s 
7) om [PO] (aF) (4 
eae | } [ol | \ 1/0 
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A] = [A®] + [64mm] (5) 


where [6F] has » nonzero rows and [64 ,..,] has m nonzero rows 
and columns. It should be emphasized that the elements of 
54mm] are not necessarily small and may be of the same order as 
theelements of [4]. With the aid of Eq. (5), Eq. (1) can now be 
rewritten as 

A] [F] = —[6Amn] [F] + [U] (6) 


t 


The vector — [6A mm] [F] is regarded as an inhomogeneous, o1 
‘constant,’ term in this equation, and it will be given a new 
symbol, 

[5Um] = —[5Anmm] [F] (7 


The significance of the method rests upon the fact that the de 
termination of [6U,,] requires the solution of only m, rather than 
n, simultaneous equations. Since [64,»,] has only m nonzero 
rows and columns, it is evident that only the first m components 


of [6U,]are nonzero. Eq. (6) can be rewritten as 


A°] [F°] + [A] [6F] = (8Um] + [U] (8) 
Hence, using Eq. (2) 

[A°] [6F] = [6Un] (9) 
Eq. (9) can be used to evaluate the increments in the unknowns 
An equation in which [6U,,] is the 
only unknown may be obtained by expanding [F] in Eq. (7) into 

[F°] + [6F] and evaluating [6F] from Eq. (9), 

[5Um] = —[5Amm] {[F°] + [6F]} = 

— [5A mm] [F°] — [64mm] [4°]! [6U nm) 


once [5U] is calculated. 


(10) 





own 
ively 
istant 


modi 


to be | 
ng in | 
con 


spec 


Only the first » rows and columns of the inverse matrix [A °] 
Let (Baw! 
be the matrix of the first m rows and columns of [A°]~! and let 
I,] be the identity matrix of the first m rows and columns 
Then, 


are required, since [6U,,] has only m nonzero values. 


oh m| [Bmm]} [5Um] = —[5Amm]°([F°] (11) 


The solution of the original problem, Eq. (1), thus requires the 
determination of the influence numbers of the first m rows and 
columns of the unmodified structure, [Bm], and the solutions of 
Eqs. (11) and (9). The inversion of Eq. (9) requires all the rows 
of [4°]! but only the first m columns. It is evident that little 
will be gained by the method if m is not a great deal smaller than 
n and that methods of solving the unmodified structure which 
If a program 


(6A, 


easily yield the influence numbers are desirable.‘ 
involving many structural modifications is to be undertaken, it 
is evidently desirable to write the equations in such a way that 
the modifications appear in as few of the coefficients as possible. 

Some interesting observations can be made concerning the 
case where there is only a single structural modification, which 
may be assumed to occur in the first row and the first column 
In this case Eq. (11) becomes a scalar equation for which the solu 
tion can be immediately written as 


6U; = —b6AyF,°/(11 + 6An By) (12) 
The increments in the unknowns are from Eq. (9), 
6F, = —By 6AyF,°/1 + 6An* Bu) (13) 
From this equation we observe that 
6Fy/6F; = Bu / Bi (14) 


In other words the ratio of the changes in any two unknowns is 
independent of the size of the structural modification. 
If we write R = B,,+6A,,; and let k = 1 in Eq. (13), then 


5F,/F,° = —R/(1 + R) (15) 


lf 6Ay, is extremely small (R< 1) then there is a direct pro 
portionality between the per unit change in F; and 6Ay. Eq 
15) states that this proportionality is to be altered by the factor 


1/11 + R) where R is the product of the change in the coefficient 
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and the corresponding element of the inverse matrix. Hence, 
once By is determined, the effect of a coefficient change of any 
size is dctermined by a simple curve. From Eq. (13) it is evident 
that the only way in which the external loads enter into the 
problem of determining the effect of a single structural modifica 
tion is through the quantity F\°. Hence, results for one loading 
condition can be immediately extended to all other loading condi 
tions if the solutions for the unmodified structure under the other 
loading conditions are known. 

The above results can sometimes be used in experimental or 
analog computer investigations to reduce the number of cases to 
be studied. For example, the effect of variations in the rigidity 
of a structural element can be determined for all values of its 
stiffness and all loading conditions by making measurements for 
the standard value of the rigidity under all loading conditions 
and for one other value under one loading condition. The non- 
standard value should preferably represent a large change 
Then, from Eqs. (14) and (15), 


oF; 1 . 
By = = — ro we, (16) 
6Ay F,° + bF, PAn 
Br = By *(6Fy/5F;) (17) 
Eq. (13) then gives 6F; for any values of 6A, and F,°. For ex- 


ample, if the structural element is a member of length land area 
Ao carrying axial forces only and F, is the force in that member, 
then the coefficient of F,° is 

Ay°® = (l/EAo) + C 
where C may be regarded as the flexibility of the path by which F\° 
If the structural modi 


(18) 


traverses the remainder of the structure. 
fication consists of changing the area of the member to a new 


san = we [(7.) - (Ga) 


From the generality of the method the other unknowns, 
5Fy, may be axial forces, moments, displacements, ete. The 
only restriction is the trivial one that the quantity /, must be 
statically indeterminate. 


value, A;, then 


(19) 
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pees EXACT SOLUTION of the energy equation for a rotating plate’ 
may be extended without any limitation of generality to all 
rotationally symmetrical flows for which the exact or approximate 
solution of the Navier-Stokes equations may be represented in the 
same form as that of the von Karmn solution for the flow over a 


rotating plate--i.e., in the form, 
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rw F(¢) 
rwG($) 
w= Vw H(¢) 


For some of these flows, tables have been computed of the func 
tions F, G, and H so that they may be used to obtain function 
tables for the solution of the energy equation. 
to tabulate the functions S and Q! also for the following flows: 

(1) Rotational flow over solid ground.? 

(2) Flow over a rotating plate exposed to axial airflow,* where 
various heat-transfer effects of the parameter ratio of the inflow 
velocity and the angular velocity of the flow may be considered 
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wu A RIGID AIRPLANE MODEL flying in a closed wind tunnel 


is given freedom in pitch and heave, it possesses two oscil- 
latory modes. The shorter period mode corresponds to the usual 
short-period longitudinal oscillation in free flight and is commonl+ 
heavily damped. The longer mode is caused by tunnel wall 
interference and is commonly divergent.' This unstable motion 
is called ‘‘galloping,’’ and its existence seriously limits wind- 
tunnel testing of models provided with these degrees of freedom. 

Recently, a means of stabilizing this mode has been developed 
by students at Massachusetts Institute of Technology working 
under the direction of the writer. The device in question is 
shown in Fig. 1. It consists of an auxiliary airfoil attached to 
the model mount in such a way that a nose-up attitude of the 
model causes a vertically downward load to be applied to the 
mount. Since the auxiliary airfoil is pivoted about its aerody 
namic center, no moment is applied to the model by the device. 

This device has been built and tested and preliminary results 
are presented in an M.I.T. Bachelor’s Thesis by Messrs. Gion 
friddo and Stewart.? The auxiliary airfoil used had 15 per cent 
of the wing area of the airplane model, had an aspect ratio of 8, 
and was geared to the model so that a 1l-deg. nose-up change of 
attitude of the model caused a 1-deg. nose-down change of atti 
tude of the auxiliary wing. 

The stabilizing device operates in the following manner: Sup 
pose that the model, initially in equilibrium with no vertical ve 
locity, falls with a constant vertical velocity, w. If the model 
now rotates through an angle —8, equal to the ratio of the verti 
cal velocity to the tunnel velocity, it is seen that the model with 
out stabilizing device preserves its equilibrium angle of attack and 
therefore its condition of equilibrium. In the absence of tunnel 
wall effect, the model without stabilizing device is seen to be in- 
different not only to vertical position but also to vertical ve- 
locity. In the presence of tunnel wall interference the model ex- 
periences a nose-up change of trim as it approaches the floor anda 
nose-down change of trim as it approaches the ceiling. This tend- 
ency, coupled with the indifference to a certain combination of 


Thus it is possible 
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vertical velocity and pitch attitude, together with inertia effects, 
gives rise to the unstable galloping motion. With the stabilizing 
device, however, the nose-down attitude associated with a verti- 
cal velocity gives rise to a nose-up attitude of the auxiliary wing, 
which produces a lift force opposing the vertical velocity—i.e., a 
damping force—and this effect in combination with the effect of 
tunnel wall interference causes the longer period oscillation of the 
model to be a damped oscillation. 

The model flown in the tests reported in references 2 and 3 was 
a conventional high-wing monoplane of reasonable proportions. 
Without the auxiliary wing its longer period oscillation had a 
negative damping ratio of about 10 per cent critical, and with the 
auxiliary wing, a positive damping ratio of about 30 per cent 
critical. With the auxiliary wing in use, the model would center 
itself in the tunnel at its trim speed, and its behavior was docile in 
every way. 

The auxiliary wing also affects the shorter period oscillation, of 
course, but here its effect is essentially the same as a reduction of 
the lift curve slope of the model by a factor proportional to the 
ratio of auxiliary wing lift curve slope to model lift curve slope 
and the auxiliary wing gearing ratio, since the angle of attack is 
nearly in phase with the angle of pitch attitude in this mode. 
This change theoretically reduces the damping of the shorter 
mode from about 60 to 59 per cent critical damping for the model 
tested and causes a slight increase in the period of the motion. 

Since the stabilizing device described appears to permit testing 
of models with simultaneous freedoms in heave and pitch with 
little distortion of the short-period longitudinal mode, work is in 
hand to investigate its dynamics in detail, and publication of the 
results of this investigation is intended at an early date. 
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